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Abstract 

We investigate color superconducting phases of cold quark matter at densities rele- 
vant for the interiors of compact stars. At these densities, electrically neutral and 
weak-equilibrated quark matter can have unequal numbers of up, down, and strange 
quarks. The QCD interaction favors Cooper pairs that are antisymmetric in color 
and in flavor, and a crystalline color superconducting phase can occur which accom- 
modates pairing between flavors with unequal number densities. In the crystalline 
color superconductor, quarks of different flavor form Cooper pairs with nonzero to- 
tal momentum, yielding a condensate that varies in space like a sum of plane waves. 
Rotational and translational symmetry are spontaneously broken. We use a Ginzburg- 
Landau method to evaluate candidate crystal structures and predict that the favored 
structure is face-centered-cubic. We predict a robust crystalline phase with gaps 
comparable in magnitude to those of the color-flavor-locked phase that occurs when 
the flavor number densities are equal. Crystalline color superconductivity will be a 
generic feature of the QCD phase diagram, occurring wherever quark matter that is 
not color-flavor locked is to be found. If a very large flavor asymmetry forbids even 
the crystalline state, single-flavor pairing will occur; we investigate this and other 
spin-one color superconductors in a survey of generic color, flavor, and spin pairing 
channels. Our predictions for the crystalline phase may be tested in an ultracold gas 
of fermionic atoms, where a similar crystalline superfluid state can occur. If a layer 
of crystalline quark matter occurs inside of a compact star, it could pin rotational 
vortices, leading to observable pulsar glitches. 
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Chapter 1 
Introduction 



1.1 Overview 

In this thesis we shall discuss the behavior of cold quark matter at densities that 
are relevant for the interiors of compact stars. It is well known that cold dense 
quark matter is unstable to the formation of a condensate of quark Cooper pairs, 
making it a color superconductor. Various phases of color superconductivity have 
been proposed, and in section 1.2 we review the phase diagram of QCD to provide a 
larger context for a discussion of the color superconducting phases. In section 1.3 we 
discuss the color-flavor-locked (CFL) color superconductor, the ground state of cold 
quark matter at very high densities. In section 1.4 we describe how the CFL phase 
can be disrupted at intermediate densities that are relevant for compact stars. At 
these intermediate densities, neutral quark matter has unequal numbers of up, down, 
and strange quarks, and a crystalline color superconducting phase is favored. The 
crystalline color superconductor has the remarkable virtue of allowing pairing between 
quarks with unequal Fermi surfaces. Cooper pairs with nonzero total momentum are 
favored; the condensate spontaneously breaks translational and rotational invariance, 
leading to gaps which vary periodically in a crystalline pattern as a superposition 
of plane waves. In section 1.5 we discuss the crystalline phase and look ahead to 
the detailed calculations of chapters 2 and 3 where we investigate single-plane-wave 
and multiple-plane- wave crystalline phases, respectively. In section 1.6 we discuss 
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single-flavor color superconductivity, which might occur when there is a very large 
flavor asymmetry that forbids the crystalline state. We also preview chapter 4, which 
presents a larger survey of various single-color and single-flavor color superconductors. 
Many of these are spin-one phases with unusual spectra of elementary excitations. 
Finally, in section 1.7 we discuss physical contexts in which the crystalline phase may 
occur with observable consequences. In 1.7.1 we review the astrophysical implications 
of color superconductivity for compact stars. If a layer of crystalline quark matter 
occurs inside of a compact star, it could pin rotational vortices, leading to observable 
pulsar glitches. In 1.7.2, we describe how a crystalline superfluid might be created 
and detected in a trapped gas of ultracold fermions. These are previews of more 
detailed discussions of glitches and cold atoms that appear in chapter 5. 



1.2 The phase diagram of QCD 

In recent years much theoretical and experimental effort has been devoted to under- 
standing the behavior of quantum chromodynamics (QCD) in extreme conditions of 
very high temperature or density. Because QCD is asymptotically free [1], its high 
temperature and high density phases are readily described in terms of quark and 
gluon degrees of freedom [2]. At high temperature, the familiar hadronic phase of 
QCD gives way to a deconfined plasma of quarks and gluons, in which all the sym- 
metries of the QCD Lagrangian are unbroken [3] . This phase preceded hadronization 
in the early universe, and efforts are underway to produce and probe this phase in 
thermalized collisions of relativistic heavy ions at Brookhaven and CERN labora- 
tories [4]. At low temperature and high density, on the other hand, the hadronic 
phase gives way to a degenerate Fermi system of quarks. Under the influence of 
the QCD interaction, quarks near the Fermi surface can bind together as Cooper 
pairs, which condense by the BCS mechanism [5] to form a color superconduc- 
tor [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. While not accessible in the lab- 
oratory, this cold dense quark matter might occur inside compact stars, with a host 
of potential astrophysical implications [19]. 
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Figure 1-1: Schematic QCD phase diagram with hadronic, plasma, and superconduct- 
ing phases, as a function of temperature T and chemical potential /i qua rk = /-tbaryon/3. 
The gray-shaded areas are regions of phenomenological interest. 

These various states of QCD fit together in a phase diagram as a function of 
temperature and chemical potential, as shown in figure 1-1. Let us first consider the 
temperature axis of this phase diagram. The early universe moved down the vertical 
axis during the first tens of microseconds after the big bang [20] . The temperature axis 
can also be studied with lattice simulations, which reveal a transition from a quark- 
gluon plasma (QGP) to a gas of hadrons at a temperature of about 170 MeV [21]. 
With realistic quark masses (m s ^> m U ( i ^ 0), chiral symmetry is explicitly broken 
everywhere and the transition is a rapid but smooth crossover. Chiral symmetry is 
(approximately) restored in the QGP phase and broken in the hadronic phase. 

Now let us consider the chemical potential axis. Nuclei are droplets of the hadron 
liquid phase; they sit on the horizontal axis just to the right of the small "curlicue" 
which is a line of first-order phase transitions between the hadron gas and the hadron 
liquid (this line terminates at a second-order critical point at a temperature of about 
10 MeV) [22]. Moving further to the right along the chemical potential axis, the 
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density increases beyond nuclear density and eventually the nucleons overlap to such 
an extent that a description in terms of quark degrees of freedom is more suitable. In 
fact various model calculations suggest that on the horizontal axis there is a first-order 
transition from the hadronic phase to deconfined quark matter [9, 10, 23, 24, 25, 26, 
27]. Starting from this first-order transition at fi q h, a first-order line separating quark 
and hadron phases extends upwards and to the left, terminating at a second-order 
critical point before it reaches the temperature axis [23, 24]. In relativistic heavy ion 
experiments at Brookhaven and CERN, the collided nuclei might follow trajectories in 
T-fi space like those shown in figure 1-1: the nuclei start at zero temperature, depart 
from equilibrium during the first moments of collision, then perhaps reappear on the 
phase diagram at high temperature, thermalizing to create a brief fireball of quark- 
gluon plasma before expanding and cooling through the quark-hadron transition to 
produce hadrons. 

To the right of jjt, q h is the regime of color superconductivity. In this regime decon- 
fined quarks fill large Fermi seas and the interesting physics is that of quarks near 
their Fermi surfaces (quarks that are deep within Fermi seas are Pauli-blocked and 
therefore essentially behave as free particles). Pairs of quarks near the Fermi surface 
that are antisymmetric in color feel an attractive QCD interaction; this is not sur- 
prising because at low density and temperature the quarks bind strongly together to 
form baryons. This attractive interaction makes the system unstable to the formation 
of a BCS condensate of quark-quark Cooper pairs. Because a pair of quarks cannot 
be a color singlet, the BCS condensate breaks color gauge symmetry, and the system 
is a color superconductor. 

The phase boundary separating the QGP and color superconducting phases in fig- 
ure 1-1 occurs at the critical temperature T c (/x) above which the di quark condensates 
vanish. This temperature is expected to be a few tens of MeV [18]. Hot protoneutron 
stars may cool through the color superconducting phase a few seconds after the birth 
of the star, with interesting implications for neutrino transport in supernovae [28]. 
Compact stars rapidly cool to temperatures of a few keV and then reside on the chem- 
ical potential axis of the QCD phase diagram (at essentially zero temperature) [29]. 



12 



The chemical potential will vary with depth and the star could be a "hybrid star" 
with a hadronic mantle and a quark matter core [30, 31, 32, 33]. Because the star 
temperature is small compared to the critical temperature for color superconductivity, 
if any quark matter is present then it will be a color superconductor. 

At very high densities (far to the right on the phase diagram) the ground state of 
QCD is the color-flavor-locked (CFL) color superconductor [11, 34, 15, 18]. This is 
a robust and symmetric phase that is invariant under simultaneous SU(3) rotations 
in color and flavor. At very high densities the strange quark mass can be neglected 
and flavor SU (3) is a good symmetry of the QCD Lagrangian. The nonzero strange 
quark mass explicitly breaks this symmetry and the effect is more pronounced at lower 
densities. Moving down the chemical potential axis from very high density, this flavor 
asymmetry can eventually disrupt the CFL phase at an "unlocking" transition [35, 
36, 37, 38, 39, 40, 26]. It is uncertain whether unlocking occurs before hadronization. 
If, as in figure 1-1, we assume that unlocking does occur first, then there is a window 
in the QCD phase diagram at intermediate density (fi q h < < /U un i oc k) for which the 
ground state of QCD is deconfined quark matter that is not color-flavor-locked. 

In unlocked quark matter, pairing should still occur because there is still an attrac- 
tive interaction between quarks. Any pairing breaks color symmetry so the system 
remains a color superconductor. In the present work we investigate the novel pair- 
ing patterns of quark matter that can occur within this window of non-CFL color 
superconductivity. We propose that the most favorable candidate is crystalline color 
superconductivity [41, 42, 43, 44, 45, 46, 47, 48]. Unlocked quark matter has dif- 
ferent number densities of u, d, and s quarks, and we will show that the crystalline 
phase can accommodate this by forming Cooper pairs with nonzero total momentum. 
Condensates of this sort spontaneously break translational and rotational invariance, 
leading to gaps which vary periodically in a crystalline pattern. 
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1.3 High density and color-flavor- locking 



At extremely high densities, the quarks at the Fermi surface have very large mo- 
menta and their interactions are asymptotically weak. In this limit, a rigorous the- 
oretical analysis of the QCD ground state is possible using weak-coupling, but non- 
perturbative, methods of BCS theory [6, 12, 16, 13, 49, 50, 51, 52, 53, 54, 17, 55]. This 
analysis predicts that at asymptotically-high density, the preferred ground state of 
QCD is the color-flavor-locked (CFL) color superconductor involving three massless 
flavors of quarks (at asymptotic densities, one can reasonably neglect the up, down, 
and strange quark masses). All nine quarks (three flavors, three colors) together form 
a simple and elegant diquark condensate of the form [11] 

(?/w(x)VW x )) ^ Aoe Q/L 4e iiA (C , 75)aft (1-1) 

with indices for color (a, (3), flavor and spin (a, b). The e a pA tensor indicates 

that the condensate is an SU(3) color antitriplet of rg, rb, and gb Cooper pairs. 
The color 3 channel is favored because this is the attractive channel for perturbative 
single-gluon exchange. The e^A tensor indicates that the condensate is simultaneously 
an SU(3) flavor antitriplet of ud, us, and ds Cooper pairs. The common index A is 
summed and therefore "locks" color to flavor. The Dirac matrix C75 indicates that 
the condensate is both rotationally invariant (J = 0) and parity even. Quarks are 
paired with the same helicity and opposite momentum; therefore they have opposite 
spin and form a spin singlet. 

The value of the CFL gap Ao can be calculated in the asymptotic limit, where 
the leading interaction is just single gluon exchange. Unfortunately the result can be 
trusted only at unphysically large chemical potentials, of order 10 8 MeV or higher [56]. 
Extrapolation to densities of interest for compact stars (fi ~ 400 MeV) is unjustified 
but yields a gap of about 10-100 MeV. Alternatively, the gap can be calculated by 
using phenomenological toy models whose free parameters are chosen to give rea- 
sonable vacuum physics [9, 10, 23, 11, 14, 15, 57, 58]. For example, one might use 
a Nambu-Jona-Lasinio (NJL) model in which the interaction between quarks is re- 
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placed by a pointlike four-fermion interaction (with the quantum numbers of single- 
gluon exchange or the instanton interaction), and choose the coupling constant to 
fit the magnitude of the vacuum chiral condensate at \x = 0. It is gratifying that 
these model calculations also yield gaps of about 10-100 MeV, in agreement with the 
extrapolation from the asymptotic limit. 

The CFL phase has a remarkable pattern of symmetry breaking and elementary 
excitations [11]. The exact microscopic symmetry group of QCD with three massless 
flavors is 

G micr osco P ic = SU(3) color x SU(3) L x SU(3) R x U(1) B (1.2) 

where the first factor is the color gauge symmetry, the next two factors are the chiral 
flavor symmetries, and the last factor is baryon number. Electromagnetism is included 
by gauging a vector U (1) subgroup of the flavor groups. In the CFL phase, G m icrosco P ic 
is broken to the subgroup 

Gcfl = SU(3) co \ ot+ l+b. x Z 2 . (1.3) 

The color and chiral flavor symmetries are broken to a "diagonal" global symmetry 
group of equal SU(3) transformations in all three sectors (color, left-handed flavor, 
and right-handed flavor). All nine quarks are gapped and the system has nine massive 
spin- 1 quasiquark excitations that decompose into a octet and a singlet under the 
diagonal SU(3). The gluon fields produce an SU(3) octet of spin-1 bosonic excitations 
that acquire a mass by the Meissner-Higgs mechanism. Chiral symmetry is broken 
and there is an associated octet of massless psuedoscalar bosons. Baryon number 
is broken to a discrete Z 2 symmetry in which all the quark fields are multiplied by 
— 1; there is an associated massless scalar boson which is the superfluid mode of the 
CFL phase. The gauge symmetries are completely broken except for a residual U(l) 
generated by Q = Q + -^T 8 where Q is the electromagnetic charge operator and T 8 
is a generator in the Lie algebra of the SU(3) color group. This U(1)q symmetry 
is a gauged subgroup of the unbroken SU(3) co \ ot+ l+r- There is a corresponding Q 
photon which is a linear combination of the usual photon and the T§ gluon (in fact 
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there is just a small admixture of the gluon, so an ordinary magnetic field externally 
applied to a chunk of CFL matter is mostly admitted as a Q magnetic field and only 
a small fraction of the flux is expelled by the Meissner effect [59]). All the elementary 
excitations have integer Q charges. 

1.4 Intermediate density and unlocking 

The CFL state pairs quarks with different flavors, forming a flavor antitriplet of ud, 
us, and ds pairs. As in any BCS state, a quark with momentum p pairs with a quark 
of momentum — p, and the condensate is dominated by those pairs for which each 
quark is in the vicinity of its Fermi surface: | |p| — pp\ < A , where A is the BCS gap 
parameter. Since different flavors are paired together, the different flavors must have 
equal Fermi surfaces for the BCS pairing scheme to work. If the Fermi momenta are 
different, then it is no longer possible to guarantee that the formation of pairs lowers 
the free energy: the two fermions in a pair have equal and opposite momentum, so at 
most one member of each pair can be created at its Fermi surface. This implies that 
the CFL phase can be disrupted by any flavor asymmetry that would, in the absence 
of pairing, separate the Fermi surfaces. The nonzero mass of the strange quark has 
precisely this effect. To understand this, consider noninteracting quark matter with 
m u = md — and m s ^ 0. Chemical equilibrium under weak decay reactions requires 

2 1 

Vu= A* - g^e, (J, d = Us = V + gAte (1-4) 

where (i = |/ibaryon is the average quark chemical potential, and (i e is a chemical 
potential for electrons. With these chemical potentials the noninteracting quarks and 
electrons fill Fermi seas up to Fermi momenta given by 

P U / = (iu,dl Pf = Vv 2 s - m2 s, PF = (1.5) 
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with corresponding number densities 



^M=A<* N s = \^ 2 s -m 2 s f'\ N e =^\. (1.6) 

Electrical neutrality then requires 

2 -N u - l -N d - l -N s - N e = 0. (1.7) 

These equations are readily solved and figure 1-2 shows how the Fermi momenta of 
the quarks and electrons respond as you vary m s at fixed /i. Notice that the effect 
of the strange quark mass, combined with the requirement of electric neutrality, is 
to push pf, up and pp down relative to p F , and at the same time induce a nonzero 
electron density. To leading order in m s the electron chemical potential is 

A*e = ^ (1.8) 

and the Fermi momenta are 



FF P 12/i FF 4/i 

u m l 
Pp = U i 

5m|_ m| 

12/1 4/1 (1J) 

The magnitude of the splitting between Fermi surfaces is 5pF = m^/A/i. Notice 
that decreasing /i enhances this flavor disparity so the effect is more important at 
intermediate densities. 

The system responds differently when pairing interactions at the Fermi surface 
are taken into account [37]. At asymptotically high densities, the system is in the 
CFL state with equal numbers of u, d, and s quarks. As /i decreases, the CFL state 
remains "rigid" with equal numbers of u, d, and s quarks, despite the presence of a 
stress which seeks to separate the Fermi surfaces. This rigidity maximizes the binding 
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Figure 1-2: Fermi momenta for noninteracting electrons and quarks, in a system 
that is electrically neutral and in weak equilibrium. All three quark Fermi surfaces 
separate when the strange quark mass is nonzero, with 5pp ~ m^/4/i. 

energy for the Cooper pairs in the CFL phase. The CFL state is the stable ground 
state of the system only when its negative interaction energy offsets the large positive 
free energy cost associated with forcing the Fermi seas to deviate from their normal 
state distributions. The free energy of the CFL state must be compared to that of 
the unpaired or "normal" state in which the quarks simply distribute themselves in 
Fermi seas as in figure 1-2. The result is [40] 

3 3 

^CFL - ^normal ~ + T^ m t ( L1 °) 

where the negative first term represents the pairing energy gain of the CFL phase, 
and the positive second term is the cost associated with enforcing equal numbers of 
u, d, and s quarks. We have neglected terms in the free energy that are of high order 
in A or m s . We find that the CFL phase is favored over unpaired quark matter only 
for [35, 36, 37, 39, 40, 26] 

2 

Tfl 

H > /Amiock ~ -rr- (1-11) 

and the CFL pairing vanishes at a first-order unlocking transition. Notice that the 
unlocking occurs when the Fermi surface splitting SpF in the unpaired phase becomes 
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equal to the gap in the CFL phase: this is consistent with the physical intuition 
that the CFL phase can "force" the Fermi seas to deviate from their normal state 
distributions only for 5pF < Aq. 

In interpreting equation (1.11), recall that the value of the CFL gap is not precisely 
known: it is of order 10-100 MeV and is also density-dependent. The strange quark 
mass parameter m s includes the contribution from any (ss) condensate induced by the 
nonzero current strange quark mass, making it a density-dependent effective mass. 
At densities that may occur at the center of compact stars, corresponding to fi ~ 
400 — 500 MeV, m s is certainly significantly larger than the current quark mass, and 
its value is not well known. In fact, m s decreases discontinuously at the unlocking 
transition [26]. Thus, the criterion (1.11) can only be used as a rough guide to 
the location of the unlocking transition in nature. As in figure 1-1 we assume that 
unlocking occurs before hadronization, so there is a window of intermediate densities 
where non-CFL quark matter will occur. 

It has been proposed [60] that the CFL state may not be completely rigid above 
the unlocking transition, but may instead respond to the imposed stress by forming 
a condensate of CFL Goldstone bosons. To respond to the stress the system wants 
to reduce the number of strange quarks and increase the number of up quarks, as 
in figure 1-2. Introducing an extra up quark and a strange hole lowers the energy 
by ~ m 2 s /2fi, but appears to require the breaking of a pair and therefore involves an 
energy cost which is of the order of the gap Ao- However, a down-particle/strange- 
hole pair has the quantum numbers of a kaon, so the energy cost is actually just 
the mass of the kaon-like elementary excitation in the CFL phase. The CFL kaon is 
a pseudo-Goldstone boson which acquires a small mass rriK ~ Aoy/m s md/ \i <C Ao 
when nonzero quark masses are introduced which explicitly break chiral symmetry. 
So the CFL vacuum may lower its free energy by decaying into K° collective modes 
by the process — > (us)(du). A K condensate corresponds to a relative rotation of 
left- and right-handed condensates in flavor space. If kaon condensation occurs, it 
lowers the free energy of the CFL phase with another term of order m\ in equation 
(1.10). As a result the 4 in the denominator of (1.11) increases but remains smaller 
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than 5, so this is a relatively minor effect as it concerns unlocking [40]. 

The unlocking transition was originally studied without the requirement of charge 
neutrality [35, 36]. It was assumed that fi u = fid = fi s = fi and therefore in the 
unpaired phase the up and down quarks would have equal Fermi momenta p F = p F = 
fi but the strange quark would have a smaller Fermi momentum p s F = a//i 2 — m^, so 
the system would have a net positive charge density. In this context, us and ds pairing 
are disrupted at the unlocking transition, but ud pairing should persist because there 
is no "stress" trying to separate the u and d Fermi surfaces. What is left over is the 
simple and well-known "2SC" state [9, 10, 18] with a condensate of the form 

(ip iaa (x.)ip jf3b (x)) oc Aoe Q/g3 ei i 3(C , 7 5 ) ab . (1.12) 

This is quite similar to equation (1.1) for the CFL condensate, except that there is 
no summation over a common color and flavor antitriplet index A. In the 2SC phase 
the condensate involves only two colors and two flavors. Four of the nine quarks are 
paired: rd pairs with gu, and ru pairs with gd. Five quarks are left unpaired (the 
blue u and d quarks, and all three colors of the strange quark). 

In a charge neutral system, however, the 2SC phase is unlikely to occur. As we 
have seen, imposing charge neutrality in the unpaired phase requires the introduction 
of an electron chemical potential /z e . This is an isospin chemical potential which 
separates the up and down Fermi surfaces, so the stress that disrupts us and ds pairing 
should equally disrupt ud pairing. A detailed calculation confirms that this physical 
intuition is correct [40, 61]: a calculation of the free energies of charge-neutral CFL, 
2SC, and normal (unpaired) states reveals that the 2SC phase is nowhere favored. 

Below the unlocking transition, the system is still unstable to the formation of 
a condensate of Cooper pairs, but the pairing pattern must be something other 
than CFL or 2SC. Perhaps the most obvious possibility is that the system simply 
abandons inter-species pairing and forms single-flavor (uu), (dd), and (ss) conden- 
sates [8, 62, 63, 64]. Unfortunately these condensates are rather feeble: the gaps 
are no larger than a few MeV [63, 62], and they could even be much smaller than 
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this [62, 9]. The single-flavor condensates are weak because they must be either 
symmetric in spin (and therefore J = 1) or symmetric in color, whereas the QCD 
interaction favors condensates that are both antisymmetric in spin ( J = 0) and anti- 
symmetric in color. This is easy to understand. The color 3 channel is the dominant 
attractive channel, perturbatively (with Coulombic gluons), nonperturbatively (via 
the instanton interaction), and phenomenologically (from the fact that pairs of va- 
lence quarks in a baryon are in a color 3 state). The J = channel is enhanced by 
its rotational symmetry (larger symmetry generally implies more robust pairing). 

Condensates that are antisymmetric in color and spin are also antisymmetric 
in flavor (by the Pauli principle), so the QCD interaction naturally favors inter- 
species (inter-flavor) pairing. So, inter-species pairing is likely to be favored, but 
a novel pairing arrangement must be proposed that can accommodate inter-species 
pairing even at intermediate densities when the different species have different Fermi 
momenta. The pairing arrangement of the crystalline color superconductor [41, 42, 
43, 44, 45, 46, 47, 48] is the most well-studied option, and it is the main emphasis of 
the current work. 

The crystalline phase was originally described by Larkin, Ovchinnikov, Fulde, and 
Ferrell (LOFF) [65, 66] as a novel pairing mechanism for an electron superconductor 
with a Zeeman splitting between spin-up and spin-down Fermi surfaces, neglecting 
all orbital effects of the magnetic perturbation. Quark matter is a more natural 
setting for the LOFF phase, as it features a "flavor Zeeman effect" with no orbital 
complications, as a consequence of the large strange quark mass. Cooper pairs in the 
LOFF phase have nonzero total momentum: a quark with momentum p is paired 
with a quark with momentum — p + 2q such that each quark is near its respective 
Fermi surface, even though the two Fermi surfaces are disjoint. The magnitude |q| is 
determined by the separation between Fermi surfaces (we expect |q| ~ Spp) while the 
direction q is chosen spontaneously. This generalization of the pairing ansatz (beyond 
BCS ansatze in which only quarks with momenta which add to zero pair) is favored 
because it gives rise to a region of phase space where each quark in a pair resides 
near its respective Fermi surface; as a result, pairs can be created at a low cost in 
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free energy and a condensate can form. In contrast to the BCS phase, where pairing 
occurs over the entire Fermi surface, LOFF pairing has a restricted phase space. For 
a given q, the quarks that pair are only those in "pairing rings" , one on each Fermi 
surface; as explained below and in figure 1-3, these circular bands are antipodal to 
each other and perpendicular to q. 

As we shall see in chapter 2, the phase space is restricted to the ring-like pair- 
ing regions by the formation of "blocking regions" (see figure 2-2) in which pairing 
is forbidden. Momentum modes inside these blocking regions are occupied by one 
species of quark but not the other, so pairing cannot occur. Quasiquark excitations 
are gapless for momenta at the boundaries of these blocking regions, but within the 
pairing rings the quasiquarks are gapped. Because it has both gapped and gapless 
quasiparticle excitations, the crystalline state is simultaneously superconducting and 
metallic. 

If each Cooper pair in the condensate carries the same total momentum 2q, then 
in position space the condensate varies like a plane wave: 



meaning that translational and rotational symmetry are spontaneously broken. This 
justifies calling it a crystalline color superconductor. Of course, if the system is 
unstable to the formation of a single plane-wave condensate, then we expect that 
a superposition of multiple plane waves is still more favorable, leading to a more 
complicated spatial variation of the condensate: 



Each A q corresponds to condensation of Cooper pairs with momentum 2q, i.e. another 
pairing ring on each Fermi surface. As we add more plane waves, we utilize more of the 
Fermi surface for pairing, with a corresponding gain in condensation energy. On the 
other hand, the rings can "interact" with each other: condensation in one mode can 



(^(x)^(x)) ~ Ae 2iqx 
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enhance or deter condensation in another mode. The true ground state of the system 
is obtained by exploring the infinite-dimensional parameter space of crystalline order 
parameters {A q } to find that particular crystal structure which is a global minimum 
of the free energy functional f2[A(x)] = f2({A q }). 

As an aside, it is worth noting that crystalline phases have appeared in other QCD 
contexts. In their analysis of quark matter with a very large isospin density (with large 
Fermi momenta for down and anti-up quarks) Son and Stephanov have noted that if 
the d and u Fermi momenta differ suitably, a LOFF crystalline phase will arise [67]. 
A LOFF crystalline phase can also occur for neutron-proton pairing in asymmetric 
nuclear matter with a splitting between the neutron and proton Fermi surfaces [68]. 
Moreover, the LOFF state is not the only crystalline phase that has been investigated: 
at large baryon number density, pairing between quarks and holes with nonzero total 
momentum has also been discussed [69, 70, 71, 72]. This results in a chiral condensate 
which varies in space with a wave number equal to 2/i; in contrast, the LOFF phase 
describes a diquark condensate which varies with a wave number 2|q| comparable to 
5pF- Several possible crystal structures have been analyzed for the crystalline chiral 
condensate [72], but this phase is favored at asymptotically high densities only if the 
number of colors is very large [69], greater than about iV c = 1000 [70, 71]. It may 
arise at lower densities in QCD with fewer colors, but apparently not in QCD with 
N c = 3 [72]. 

At least two alternatives to the crystalline color superconducting phase have been 
proposed that also allow pairing between quarks with unequal Fermi surfaces. The 
first alternative is the deformed Fermi sphere (DFS) superconductor [73]. In this 
phase, the unequal Fermi surfaces of the two paired species are deformed so that 
they can intersect, and then pairing can occur in the vicinity of this intersection. 
The deformations are volume-conserving so that they do not change particle numbers 
for the two species. The larger Fermi surface has a prolate deformation, while the 
smaller Fermi surface has an oblate deformation, and pairing occurs along two bands 
just above and below the equator of each spheroid. The DFS phase breaks rotational 
symmetry, but unlike the crystalline (LOFF) phase it does not break translational 
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symmetry because the Cooper pairs still have zero total momentum. In both the DFS 
and LOFF phases the pairing is along circular bands on the Fermi surfaces and there- 
fore the two phases would seem to have similar condensation energies. Meanwhile, in 
the DFS phase there is also a large kinetic energy cost associated with deforming the 
Fermi spheres away from their preferred spherical shapes. There is no such cost for 
the crystalline phase, so we expect the crystalline phase to have a lower free energy. 

The second alternative is the "breached-pair" color superconductor [74, 75, 76]. 
The breached-pair superconductor is translationally and rotationally invariant, in 
contrast to the crystalline and DFS phases. This state was first encountered by 
Sarma [74] in the context of an electron superconductor with a Zeeman splitting 
between the spin-up and spin-down Fermi surfaces, the same context in which the 
crystalline (LOFF) phase was first proposed. In this historical context it was found 
that the breached-pair state was never a minimum of the free energy, but recent de- 
velopments [75, 76] suggest that the breached-pair state might be a stable ground 
state for pairing between a light species and a heavy species with different Fermi mo- 
menta, and therefore might accommodate us and ds pairing in intermediate-density 
quark matter. In the breached-pair superconductor, the Fermi sea of heavy quarks 
(Fermi momentum ph) is redistributed to accommodate pairing at the light quark 
Fermi surface (Fermi momentum p{). The redistribution has a small energy cost be- 
cause the heavy quark has a very flat single-particle dispersion relation. For pi < pj L , 
heavy quarks are promoted from |p| ~ pi to |p| — Ph, creating a "trench" of unoc- 
cupied heavy quark states in the interior of the heavy quark Fermi sea. This trench 
is coincident with the light quark Fermi surface and allows the formation a conden- 
sate of Cooper pairs at this surface, a so-called "interior gap" [75]. For ph < pi (the 
scenario of interest for quark matter), heavy quarks are promoted from |p| ~ ph to 
|p| ~ pi, creating a "berm" of occupied heavy quark states far above the heavy quark 
Fermi sea. This berm is coincident with the light quark Fermi surface and allows 
the formation of a condensate of Cooper pairs at this surface, a so-called "exterior 
gap" [76]. 

The common mechanism in both interior and exterior gap phases is the promo- 
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tion of a shell of heavy quarks across a momentum "breach" of magnitude \pi — Ph\, 
thereby creating a second edge in the momentum distribution of heavy quarks. This 
edge behaves like a new Fermi surface coincident with the light quark Fermi surface, 
accommodating the formation of Cooper pairs. The breach is a "blocking region" 
analogous to the aforementioned blocking regions in the crystalline phase. In the 
crystalline phase, the blocking regions restrict the pairing to occur on rings and for- 
bid pairing away from these rings. In the breached-pair phase, the blocking region 
is spherically symmetric and forbids pairing in the breach between pi and ph- Just 
as in the blocking regions of the crystal, momentum modes within the breach are 
occupied by one species of quark but not the other. In either context, quasiparticle 
excitations are gapless for momenta at the boundaries of the blocking regions. In the 
breached-pair state this means that the light quark Fermi surface is gapped while the 
heavy quark Fermi surface remains ungapped, and the system is simultaneously su- 
perconducting and metallic, like the crystalline state. The condensation energy must 
be weighed against the cost of promotion across the breach: the cost of promotion 
is small only if the heavy quark dispersion relation is sufficiently flat, so exterior gap 
(us) and (ds) condensates can occur only when the strange quark is nonrelativistic. 
A breached-pair state has been proposed for ud pairing [77] , but the preceding argu- 
ment suggests that this is not a stable ground state. A gapless CFL state, also with 
breached pairing, has been investigated and was found to be metastable [78]. The 
evaluation of the energy of breached-pair phases is subtle, and the stability is sen- 
sitive to whether a microcanonical (fixed number density) or grand canonical (fixed 
chemical potential) approach is used [79, 76]. 



1.5 Crystalline color superconductivity 

Crystalline color superconductivity has only been studied in simplified models with 
pairing between two quark species whose Fermi momenta are pushed apart by a 
chemical potential difference [41, 42, 48, 46, 47] or a mass difference [45]. We suspect 
that in reality, in three-flavor quark matter whose unpaired Fermi momenta are split 
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as in (1.9), the pattern of pairing in the crystalline phase will involve ud, us and ds 
pairs, with color and flavor quantum numbers just as in the CFL phase. However, 
studying the simpler two-flavor problem should elucidate the nature of the crystalline 
ground state, including its crystal structure. We therefore simplify the color-flavor 
pattern to one involving massless u and d quarks only, with Fermi momenta split by 
introducing chemical potentials 

fid = P> + SfJ, 

fi u = p-8fj, . (1-15) 

In this toy model, we vary 8fi by hand. In three-flavor quark matter, the analogue of 
Sfi is controlled by the nonzero strange quark mass and the requirement of electrical 
neutrality and would be of order m 2 s /Afi as in (1.9). 

The LOFF crystalline phase was originally studied in the context of an electron 
superconductor with a Zeeman splitting between the spin-up and spin-down Fermi 
surfaces [65, 66]. The authors considered a magnetic perturbation AH = hif)^a z ijj but 
disregarded any orbital effects of the magnetic field. They were seeking to model the 
physics of magnetic impurities in a superconductor. Magnetic effects on the motion of 
the electrons [80] and the scattering of electrons off non-magnetic impurities [81, 82] 
disfavor the LOFF state. Although signs of the BCS to LOFF transition in the 
heavy fermion superconductor UPd2Al 3 have been reported [83], the interpretation 
of these experiments is not unambiguous [84]. It has also been suggested that the 
LOFF phase may be more easily realized in condensed matter systems which are 
two-dimensional [85, 86] or one-dimensional [87], both because the LOFF state is 
expected to occur over a wider range of the Zeeman field h than in three-dimensional 
systems and because the magnetic field applied precisely parallel to a one- or two- 
dimensional system does not affect the motion of electrons therein. Evidence for a 
LOFF phase in a quasi-two-dimensional layered organic superconductor has recently 
been reported [88]. 

None of the difficulties which have beset attempts to realize the LOFF phase in 
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a system of electrons in a magnetic field arise in the QCD context of interest to us. 
Differences between quark chemical potentials are generic and the physics which leads 
to these differences has nothing to do with the motion of the quarks. We therefore 
expect the original analysis of LOFF (without the later complications added in order 
to treat the difficulties in the condensed matter physics context) to be a good starting 
point. 

In our two-flavor model, we shall take the interaction between quarks to be point- 
like, with the quantum numbers of single-gluon exchange. This s-wave interaction 
is a reasonable starting point at accessible densities but is certainly not appropriate 
at asymptotically high density, where the interaction between quarks (by gluon ex- 
change) is dominated by forward scattering. The crystalline color superconducting 
state has been analyzed at asymptotically high densities in Refs. [46, 47]. We expect 
a qualitatively different crystalline phase in this asymptotic regime, but this may not 
be relevant for densities of interest for compact star physics. 

With this astrophysical context in mind, it is also appropriate for us to work at 
zero temperature. Compact stars that are more than a few minutes old are sev- 
eral orders of magnitude colder than the critical temperature (of order tens of MeV) 
for CFL or crystalline color superconductivity. The crystalline color superconduc- 
tor has been studied at finite temperature, and the critical temperature is given by 
T c /A(t=o) — 0.44 [42] (a result previously known in the historical LOFF context [89]). 
It is interesting that this differs from the usual BCS relation T C /A( T=0 ) — 0.57 [5]. 
The critical temperatures for the CFL and single-flavor color superconductors also 
differ from the BCS result [64]. 

Now let us consider how the crystalline phase occurs in our two-flavor model. 
Starting at 5fi = 0, the system forms a BCS superconductor with gap Ao- In fact this 
BCS superconductor is precisely the 2SC phase of equation (1 .12): the Cooper pairs 
are color antisymmetric (red pairs with green) and flavor antisymmetric (up pairs 
with down). The blue quarks are left unpaired. The up and down Fermi surfaces are 
coincident. As we begin to increase 5fi, the system exhibits a "rigidity" analogous 
to that of the CFL phase: despite the imposed stress the gap stays constant and 
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the Fermi surfaces remain coincident. The BCS state is the stable ground state of 
the system only when its negative interaction energy offsets the large positive free 
energy cost associated with forcing the Fermi seas to deviate from their normal state 
distributions. The free energy of the BCS state, relative to that of of the normal 
state in which the quarks simply distribute themselves in Fermi seas with p F = fi u 
and pf, = fid, is approximately 

1 2 

^BCS - ^normal ~ o^ltf + ^^J^ (1-16) 

where the first term is the negative pairing energy of the BCS state, and the second 
term is the cost associated with enforcing equal numbers of up and down quarks in the 
presence of the imposed stress 5fi. This result is exact only in the weak-coupling limit 
in which the gap Ao <C ft. This expression (1.16) should be compared to equation 
(1.10) for the free energy of the neutral CFL phase with the imposed stress of a 
nonzero m s . When 5fi reaches a critical value 

^!«-^A = (0.7071 ...)A , (1.17) 

the BCS phase "breaks" and the Fermi surfaces separate (again, the expression is 
exact only in the weak coupling limit in which A <C p). This is the two-flavor 
analogue of the CFL unlocking transition. (Bedaque [38] has investigated the mixed 
phase associated with this first-order transition, where the unpaired blue quarks also 
play a role.) This result was first derived by Clogston and Chandrasekhar [90] in the 
context of an electron superconductor with a Zeeman splitting. 

For 5fi > 5fi 1 , the up and down quarks have unequal Fermi surfaces and a crys- 
talline state is possible. In the simplest LOFF state, up quarks with momentum p 
are paired with down quarks with momentum — p + 2q. Each Cooper pair carries 
the same total momentum 2q. The allowed phase space for p is determined by the 
requirement that each quark in the Cooper pair should sit near its Fermi surface, i.e. 

||p|-^|<A and || -p + 2q| - fi d \ < A (1.18) 
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Figure 1-3: The LOFF pairing geometry for Cooper pairs with total momentum 2q. 
The dashed and solid spheres are the up and down quark Fermi surfaces, respectively. 
An up quark with momentum p (or p') near its Fermi surface pairs with a down quark 
with momentum — p+2q (or — p' + 2q) near its Fermi surface. The pairing is strongest 
for up quarks in a band centered on the dashed ring shown on the up Fermi surface, 
and down quarks in a band centered on the solid ring shown on the down Fermi 
surface. 

where A is the LOFF gap parameter. This phase space corresponds to a circular 
band on each Fermi surface as shown figure 1-3. As indicated in the figure, the bands 
are perpendicular to the spontaneously chosen direction q for the total momentum 
of each Cooper pair. The magnitude |q| = qo is determined energetically from the 
separation between Fermi surfaces. We shall find that the relation is go ~ 1.205/x. 

It is useful to discuss the various scales involved in the problem. The BCS gap Ao 
can be thought of as the fundamental energy scale for physics at the Fermi surface. If 
we consider the weak coupling limit in which A -C /2, then the two scales are cleanly 
separated and all the other Fermi-surface energy scales in the problem (i.e. |q|, 5fi, and 
the LOFF gap A) should be proportional to the fundamental scale A . We achieve 
this by taking a "double scaling" limit in which we choose to hold <5/i/A fixed while 
taking the Ao//i — > limit. In this double scaling limit, |q|/Ao and A/Aq also stay 
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fixed. In fact, every dimensional quantity for the Fermi-surface physics stays fixed as 
"measured" in units of A . If we fail to take the double scaling limit, instead keeping 
Sn/jl fixed as Aq/h — > 0, we would not find crystalline color superconductivity at 
weak coupling [91]. We will not always work in the double scaling limit (see, for 
example, chapter 2), but we will often quote analytic results that are exact in this 
limit (as we did in equation 1.17). In the double scaling limit, 5fi <C /2 and the two 
Fermi surfaces in figure 1-3 are very close together. The opening angles ip u and ipd of 
the two pairing bands become degenerate and take on the value 

« 2 cos" 1 f j^-J « 2 cos' 1 ^ = 67.1°. (1.19) 

The radial thickness of each pairing band is of order A, while the angular width is 
Sip ~ A/ ' \J |q| 2 — 5ji 2 ~ 1.5A/5/1. If we use double scaling then both i/jq and 8if) 
are constant because |q|/Ao, 5/x/Ao, and A/Ao are all held constant while Ao — > 0. 
Hereafter, when we speak of the "weak coupling limit" we shall always mean the 
double scaling limit. 

If all the Cooper pairs in the condensate have the same nonzero total momentum 
2q, then the condensate varies like a single plane wave in position space, as in equation 
(1.13). In chapter 2 we present a careful study of this single plane wave condensate, 
the simplest example of a LOFF phase. We show that there is a range of 5fi in which 
quark matter is unstable to the spontaneous breaking of translational invariance by 
the formation of a plane wave condensate. Of course, once one has demonstrated an 
instability to the formation of a plane wave, it is natural to expect that the state which 
actually develops has a crystalline structure consisting of multiple plane waves, as in 
equation (1.14). In chapter 3 we investigate this possibility. Larkin and Ovchinnikov 
in fact argue that the favored configuration is a crystalline condensate which varies in 
space like a one-dimensional standing wave, cos(2q • r). Such a condensate vanishes 
along nodal planes [65]. Subsequent analyses suggest that the crystal structure may be 
more complicated. Shimahara [85] has shown that in two dimensions, the LOFF state 
favors different crystal structures at different temperatures: a hexagonal crystal at low 
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temperatures, square at higher temperatures, then a triangular crystal and finally a 
one-dimensional standing wave as Larkin and Ovchinnikov suggested at temperatures 
that are higher still. In three dimensions, the question of which crystal structure is 
favored was unresolved [92]. Our analysis, shown in chapter 3, suggests that the 
favored crystal structure in three dimensions is face-centered-cubic. 

The crystalline states appear for 8fi > 5fi l . In chapter 2 we will show that 
the simplest LOFF state, a single plane wave condensate, can occur in an interval 
5fi 1 < 5fi < 5fi 2 - At Sfj, 2 there is a second-order transition from LOFF to the normal 
state (unpaired quarks). The second-order point occurs at 

Sfi 2 « (0.7544 ...)A (1.20) 

where this relation is exact in the weak coupling limit (the numerical coefficient is 
known exactly; it is the solution of a simple transcendental equation). At the second- 
order phase transition, A/Aq ^ and |q|/Ao tends to a nonzero limit, which we shall 
denote qo/A , where q ~ 0.90A ~ 1.20<5/i 2 . Near the second-order phase transition, 
the quarks that participate in the crystalline pairing lie on thin circular rings on their 
Fermi surfaces that are characterized by an opening angle ipo — 2 cos~ 1 (5/i 2 /go) — 
67.1° and an angular width 8ip that is of order A/8fi and therefore tends to zero 
as A/A — > 0. At 8 Hi there is a first-order phase transition at which the LOFF 
solution with gap A is superseded by the BCS solution with gap A . (The analogue 
in three-flavor QCD would be a LOFF window in rn^/fi, with CFL at lower rn^/fi 
(higher density) and unpaired quark matter at higher m^/fi (lower density).) These 
results are summarized in figure 1-4, where we have shown the free energies and gaps 
for the competing BCS, plane-wave LOFF, and unpaired quark matter phases (the 
figure also shows the gap and free energy for the multiple-plane-wave LOFF state, 
which we discuss below). Keep in mind that this figure is just a qualitative sketch 
which exaggerates the size of the plane- wave LOFF window [8^1,8^2]- Quantitative 
plots are shown in figure 2-4 in chapter 2. In the vicinity of the second-order critical 
point 8^21 our mean- field analysis yields a gap and free energy for the plane- wave 
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Figure 1-4: Free energies and gaps for competing states: normal, BCS, and single 
and multiple plane wave crystalline phases ("pw" and "xtal", respectively). For the 
plane wave state, the LOFF interval is [<5/x 1 ,5/x 2 ] ~ [0.707A , 0.754A ] (A is the 
BCS gap) and the transitions to BCS and normal states are first and second order, 
respectively. For the multiple plane wave (crystalline) state, the LOFF interval is 
[S/i^, Sfij] and both transitions are first order. Note that the plane wave state is 
exaggerated (compare figure 2-4). In reality (5/i 2 — <C (5^ — S/j,^) and the 
crystalline state is much more favorable than the plane wave state. 

state that obey simple power-law relations 

A pw oc y/ 5fi 2 - 5fi, Q pw oc -(5fi 2 - 5fi) 2 (1-21) 

with the expected mean field theory critical exponents for a second-order transition. 

Proceeding beyond just a single plane wave might seem to be a daunting task. We 
have to explore the infinite-dimensional parameter space of crystalline order param- 
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eters {A q } to find the unique crystal structure that is a global minimum of the free 
energy. However, we can exploit the fact that there is a second-order point 8 fa which 
indicates the onset of plane-wave instability in the system. In the vicinity of this 
second-order point, we can express the free energy as a Ginzburg-Landau potential; 
the potential is written as a series expansion of the exact free energy in powers of the 
order parameters {A q }. 

In chapter 3 we explicitly construct the Ginzburg-Landau potential and apply it 
to a large survey of candidate crystal structures. The Ginzburg-Landau calculation 
finds many crystal structures that are much more favorable than the single plane 
wave (1.13). For many crystal structures, the calculation actually predicts a strong 
first-order phase transition, at some 5 fa 3> 5 fa, between unpaired quark matter 
and a crystalline phase with a A that is comparable in magnitude to A . Once Sfi is 
reduced to 8 fa, where the single plane wave would just be beginning to develop, these 
more favorable solutions already have very robust condensation energies, perhaps even 
comparable to that of the BCS phase. Therefore they can even compete with the BCS 
phase and move the position of the first-order transition between BCS and LOFF to 
a new point 8 fa < 5 fa. All of this is shown schematically in figure 1-4. These results 
are exciting, because they suggest that the crystalline phase is much more robust 
than previously thought. However, they cannot be trusted quantitatively because 
the Ginzburg-Landau analysis is only controlled in the limit A — > 0, and we find a 
first-order phase transition to a state with A ^ 0. 

Even though it is quite a different problem, we can look for inspiration to the 
Ginzburg-Landau analysis of the crystallization of a solid from a liquid [93]. There 
too, a Ginzburg-Landau analysis predicts a first-order phase transition, and thus 
predicts its own quantitative downfall. But, qualitatively it is correct: it predicts 
the formation of a body-centered-cubic crystal and experiment shows that most el- 
ementary solids are body-centered-cubic (BCC) near their first-order crystallization 
transition. 

Thus inspired, let us look at how the Ginzburg-Landau calculation will proceed. 
We can start by writing down the most general expression for the Ginzburg-Landau 
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potential that is consistent with translational and rotational symmetry. We will 
include only the modes on the sphere |q| = qo = 1.25/i since these are the modes that 
become unstable at 5fi 2 . To order A 6 , the expression looks like 



Odd powers are not allowed because the potential is invariant under U(l) baryon 
number (which multiples every A q by a common phase). The symbol □ represents a 
set of four equal-length vectors (qi, q2, q3, q4), |qj| = go, with qi — q 2 + q 3 — q 4 = 0, 
i.e. the four vectors are joined together to form a closed (not necessarily planar) 
figure. Similarly, the symbol O represents a set of six equal-length vectors (qi, . . . , q 6 ), 
l0.il = Qo, with qi — q 2 + q3 — q4 — qs + q6 = 0, i.e. the six vectors form a closed 
"hexagon" . We sum only over closed sets of q- vectors because otherwise the Ginzburg- 
Landau potential, expressed in position space as a functional f2[A(x)], would not be 
translationally invariant. Rotational invariance implies that the coefficients </(□) and 
K(o) are the same for any two shapes related by a rigid rotation. 

The quadratic coefficient a changes sign at Sfi 2 showing the onset of the LOFF 
plane- wave instability: a ~ (<5/i — Sfi 2 )/S^i 2 . If there was no interaction between the 
different modes, they would just simply all condense at once, because they would all 
become unstable at the second-order point. The answer is more complicated than 
this because condensation in one mode can enhance or deter condensation in another 
mode. This interaction between modes is implemented in our Ginzburg-Landau po- 
tential by the higher order terms involving multiple modes; thus the coefficients J, 
K, . . . characterize the interactions between modes and thereby determine the crystal 
structure. 

As we shall see in chapter 3, these coefficients can actually be calculated from the 
microscopic theory, as loop integrals in a Nambu-Gorkov formalism. So for a candi- 
date crystal structures with all A q 's equal in magnitude, we can evaluate aggregate 




q,|q|=<?o 
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Ginzburg-Landau quartic and sextic coefficients (3 and 7 as sums over all rhombic 
and hexagonal combinations of the q's: 



/? = $>(□), 7 = 5>(0). 



(1.23) 



□ O 



Then for a crystal consisting of P plane waves we obtain 




(1.24) 



and we can then compare crystals by calculating (3 and 7 to find the structure with 
the lowest fl. 

Evaluating the quadratic coefficient a determines the location of the plane-wave 
instability point, i.e. the value of 5fi 2 - It also tells us that |q| ~ 1.205/i, which 
means that each pairing ring has an opening angle of 67.1°, as in equation (1.19). As 
mentioned above, on its own the quadratic term indicates that adding more plane 
waves (i.e. adding more pairing rings to the Fermi surface) always lowers the free 
energy. But this conclusion is modified by the higher order terms in two important 
ways: 

1. Crystal structures with intersecting pairing rings are strongly disfavored. Recall 
that each q is associated with pairing among quarks that lie on one ring of 
opening angle ipo ~ 67.1° on each Fermi surface. We find that any crystal 
structure in which such rings intersect pays a large free energy price. Therefore 
the favored crystal structures are those that feature a maximal number of rings 
"packed" onto the Fermi surface without intersections. No more than nine rings 
of opening angle 67.1° can be packed on a sphere without intersections [94, 95]. 

2. Crystal structures are favored if they have a set of q's that allow many closed 
combinations of four or six vectors, leading to many terms in the Y2n an( ^ 
J2o summations in equation (1.22). Speaking loosely, "regular" structures are 
favored over "irregular" structures. All configurations of nine nonintersecting 
rings are rather irregular, whereas if we limit ourselves to eight rings, there is 
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a regular choice which is favored by this criterion: choose eight q's pointing 
towards the corners of a cube. In fact, a deformed cube which is slightly taller 
or shorter than it is wide (a cuboid) is just as good. 

These qualitative arguments are supported by the quantitative results of our 
Ginzburg-Landau analysis, which do indeed indicate that the most favored crystal 
structure is a cuboid that is very close to a cube. This crystal structure is so fa- 
vorable that the coefficients (3 and 7 in the Ginzburg-Landau potential (equation 
(1.24)) are large and negative. In fact, we find several crystal structures with nega- 
tive coefficients, but the cube has by far the most negative (3 and 7. In other words, 
starting at the origin in the space of crystalline order parameters {A q }, the "steepest 
descent" in free energy is achieved by moving in the direction of the cube structure. 
Our Ginzburg-Landau potential is unbounded from below, so our analysis is unable 
to discover the actual free energy minimum or the value of the gap at which this 
minimum occurs. But we can reasonably presume that the lowest free energy and 
largest gap are achieved by moving in the direction of steepest descent. We could go 
on, to order A 8 or higher, until we found a Ginzburg-Landau free energy for the cube 
which is bounded from below. However, we know that this free energy would give 
a strongly first-order phase transition, meaning that the Ginzburg-Landau analysis 
would anyway not be under quantitative control. A better strategy, then, is to use 
the Ginzburg-Landau analysis to understand the physics at a qualitative level, as we 
have done. We understand qualitatively what features of the eight-plane-wave solu- 
tion make it most favorable, so the next step is to take this as an ansatz, solve the gap 
equation, and thus obtain a bounded free energy without making a Ginzburg-Landau 
approximation. This calculation is still in progress, but in figure 1-5 we show what 
the bounded free energy might look like (solid curve), compared with the unbounded 
Ginzburg-Landau free energy (dotted curve). The series of plots shows how the un- 
bounded free energy indicates a first order transition: for 5fi > 5fi 2 the quadratic 
coefficient a is positive: increasing A at fixed Sfi, the free energy should first turn up- 
wards, then downwards under the influence of the negative quartic and sextic terms, 
then eventually it will turn upwards again because it must be bounded from below. 
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Figure 1-5: Schematic free energies for the FCC crystal (solid line) and plane wave 
(dashed line), showing how they support first and second order phase transitions, 
respectively. The dotted line is the unbounded Ginzburg-Landau expression for the 
free energy of the FCC crystal. 

The resulting curve can thus generate a first-order transition as 5fi is varied, as shown 
in the figure. For comparison the dashed line shows the free energy of the plane wave 
crystal (with (3, 7 > 0), which demonstrates a typical second order transition. 

The eight q's of our most-favored crystal structure are the eight shortest vectors 
in the reciprocal lattice of a face-centered-cubic crystal. Therefore, we find that 
A(x) ~ (^(x) , 0(x)) exhibits face-centered-cubic symmetry. Explicitly, 



A(x) = 2A 



2tt. . 2tt. . 

cos — [X + y + z) + cos — [x — y + z) 
a a 



2tt, . 2tt. 

cos — (x + y — z) + cos — (— x + y + z) 
a a 



;i.25) 
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Figure 1-6: A unit cell of the LOFF face-centered-cubic crystal. The gray planes are 
surfaces where A(x) = 0. The darker surfaces are contours where A(x) = +4 A, and 
the lighter surfaces are contours where A(x) = —4 A. 



where the lattice constant (i.e. the edge length of the unit cube) is 

V3n 4.536 6.012 



q Sfj, A 



;i.26) 







where the last equality is valid at 5fi = 5/x 2 and where A is the gap of the BCS 
phase that would occur at 8/1 = 0. A unit cell of the crystal is shown in Fig. 1-6. The 
figure clearly reveals a face-centered-cubic structure. Like any crystal, the FCC crys- 
talline color superconductor should have phonon modes which are Goldstone bosons 
of spontaneously broken translation symmetry. Casalbuoni et al have formulated an 
effective theory for the LOFF phonons [48]. 
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1.6 Single-flavor color superconductivity 



If the Fermi momenta of the u, d, and s quarks are very far apart then the system has 
no choice but to abandon inter-species pairing and form single- flavor (uu), (dd), and 
(ss) condensates. In the two-flavor context of section 1.5, with a splitting 5fi between 
the u and d Fermi surfaces, single- flavor pairing will occur for 5p > 5fi^. At 
there is a first-order "crystallization" transition between the single-flavor state and 
the crystal state with ud pairing. The value of the first-order point Sfi^ is not known; 
if it were known we could estimate that an analogous crystallization transition will 
occur in three-flavor neutral quark matter when the Fermi momentum splitting Spp 
approaches 25^, i.e. when 

2 

777 

S PF = 5p* F ^-±~ 25/v (1.27) 

(The factor of two occurs because in our notation p d F — p F = 5pp = 25 jj). We expect 
that 5fi^ is appreciably larger than Ao, which implies that /i* is appreciably smaller 
than /i un iock- Crystalline quark matter occurs in the interval between /i* and /i un iock- If 
/i* is below the hadronization point fi q h then single-flavor quark matter is unlikely to 
occur and the crystalline phase will occupy the entire interval between hadronization 
and unlocking in the QCD phase diagram (figure 1-1). Otherwise, there may be a 
window just above the hadronization point in which single-flavor pairing is possible. 
The structure of a single-flavor (uu) condensate (flavor index 1) is [62] 

(Vw(x)V>#i&(x)) oc Ae a/ 335ii5}i(C7 3 ) afe (1-28) 

with indices for color (a,f3), flavor (i, j), and spin (a,b). The condensate is antisym- 
metric in color (as usual, the color 3 channel is favored because this is the attractive 
channel for the QCD interaction). Only two of the three colors pair; the choice of 
the index 3 for the color tensor e a p^ is arbitrary and excludes the blue quarks from 
pairing. The condensate is obviously symmetric in flavor. By the Pauli principle it 
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is also symmetric in Dirac indices 1 . It is a J = 1 (vector) condensate that breaks 
rotational symmetry: the Dirac matrix C73 in equation (1.28) indicates that the con- 
densate has spontaneously chosen the 3 direction in position space. The condensate is 
parity even. The quarks are paired with opposite helicity (LR pairing) and opposite 
momentum; therefore they have parallel spin and form a (s = 1, m — ±1) spin triplet 
state. 

The gap parameter A can be calculated using an NJL model with a four-fermion 
interaction vertex [62]. These calculations are shown in chapter 4. Unfortunately the 
value of the gap is drastically sensitive to the details of the effective interaction and to 
the chemical potential. It could be as large as 1 MeV, or orders of magnitude smaller 
(see the dash-dotted line in figure 4-2). For \x of 400 to 500 MeV, the NJL calculation 
predicts a gap that ranges from 0.1 to 10 keV (this illustrates the sensitivity to the 
chemical potential). The calculation is also very model dependent. In figure 4-2 the 
gap is calculated using an NJL model with pointlike magnetic gluons. Calibrated 
to give the same CFL gap, a different NJL model that includes pointlike electric 
and magnetic gluons predicts a much larger gap (by about a factor of 10); with an 
instanton interaction, no gap is predicted at all (the channel is flavor symmetric and 
there is no interaction with an instanton vertex). 

At asymptotically high density a model-independent calculation of the gap is pos- 
sible, with a perturbative gluon interaction [63, 64, 51]. Extrapolating to reasonable 
densities (/x ~ 400 MeV), the perturbative calculation predicts spin-one gaps of order 
20 keV - 1 MeV, assuming that the gap in the CFL phase is of order 10-100 MeV. 
The perturbative calculation also predicts that the condensate will have an additional 
C<Jo3 component: this channel, which is repulsive in an NJL model with pointlike glu- 
ons, becomes attractive at asymptotic density when the gluon propagator provides 
a form factor that strongly emphasizes small-angle scattering. In the Ccjqz channel, 
quarks are paired with the same helicity (LL or RR pairing) and opposite momentum; 
their spins are antiparallel and they form a (s = l,m = 0) spin triplet state. 

lr The symmetric Dirac matrices that could appear in equation 1.28 are C70, Cctoi, Ct7oi75, and 
Cji. The first is ruled out because it has no particle-particle component. The second and third are 
disfavored by our NJL model. 
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The elementary excitations of the single-flavor color superconductor are quite dif- 
ferent than those of a spin- zero phase like CFL. The quasi quark excitations have 
anisotropic dispersion relations and they are not fully gapped. For the C73 conden- 
sate, the energy gap to create a quasiquark goes to zero at the poles of the Fermi 
surface; for the CV03 condensate, the energy gap vanishes on the equator of the 
Fermi surface. When a nonzero quark mass is considered, these modes acquire a 
small gap of order mA//i [62]. The gapless or nearly-gapless excitations are likely 
to dominate transport properties of the material (viscosities, conductivities, etc.). 
The system should also have massless spin-waves which are Goldstone bosons of the 
spontaneously-broken rotational symmetry. It is interesting to note that, unlike the 
CFL phase, the single-flavor color superconductor does not have a massless "rotated 
photon" that can admit magnetic flux (there is no leftover U(1)q gauge symmetry). 
Therefore the phase exhibits an electromagnetic Meissner effect [96]. 

The condensate of equation (1.28) spontaneously chooses a spatial direction and 
a color direction. A lower free energy may obtained by making the replacement 
eo/33C73 ~~ ► ^olPaC^ia with a summation over the common color-spin index A. This is 
a "color-spin locking" (CSL) phase in which color structure is correlated with spatial 
direction [63, 64]. The symmetry breaking pattern in the CSL phase is 

SU(3U OI x SO(3)j -> SO(3) color+J (1.29) 

i.e. the color and rotation groups are broken to a diagonal subgroup of simultaneous 
global rotations in color and space. In the CSL phase the quasiquark dispersion 
relations are isotropic and fully gapped. This phase is analogous to the B phase of 
helium-3 [97] (whereas the condensate of equation (1.28) is analogous to the A phase). 
As in helium-3, there are yet other possible phases for spin-one condensates. Some of 
these other phases have been explored in refs. [63, 64]. 

In concert with the NJL model calculations for the single-flavor color superconduc- 
tor, a large catalog of color-flavor-spin channels for diquark condensation has been 
surveyed [62]. This survey is shown in chapter 4. For the survey we use an NJL 
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< 1 MeV 
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1 
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1 


< 1 MeV 


5 


1 


1 






C7075 





< 0.01 MeV 



Table 1.1: Summary of attractive channels from the NJL survey of chapter 4. The 
pairing pattern is shown in equation (1.30); N c and Nf are the numbers of colors and 
flavors that participate in the pairing. 

model that includes four-fermion interactions with the quantum numbers of electric 
gluon exchange, magnetic gluon exchange, and the two-flavor instanton, with Fermi 
couplings Ge, Gm, and Gi, respectively (see equation (4.2)). We investigate diquark 
condensates that factorize into separate color, flavor, and Dirac tensors, i.e. 

(4>iaa?Pj(3b) OC At^ij^ab- (1-30) 

Chapter 4 shows an exhaustive survey of 24 different condensates which have this 
generic form. Many of these condensates are spin-one, with interesting quasiquark 
dispersions like those discussed above. For all of the attractive channels, the NJL 
mean field theory gap equations are solved to estimate values of the gaps. Unfortu- 
nately, most of these NJL gap calculations suffer from the same drastic model de- 
pendence that inflicts the single-flavor color superconductor calculation as described 
above. The results for the five attractive channels are shown in table 1.1. With the 
notable exception of the first channel, the gap estimates in the right column should 
be interpreted very cautiously. Not only are the values of the gaps quite model- 
dependent, they are also extremely sensitive to the chemical potential: as we will see 
in chapter 4, they can vary by more than two orders of magnitude when the chemical 
potential is changed from 400 MeV to 500 MeV. The numerical estimates in the table 
should be interpreted as optimistic upper bounds for gaps which could be orders of 
magnitude smaller. 

The first channel (2 colors, 2 flavors) in table 1.1 is the familiar 2SC phase of 
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equation (1.12). This phase, and its 3-flavor, 3-color cousin (the CFL phase), have 
the largest gaps of any of the color superconducting phases. They are also the only 
phases for which the NJL gap calculations are robust. As we have discussed previously, 
other color superconducting phases are only likely to prevail when the CFL and 2SC 
phases are disrupted by a flavor asymmetry, as occurs in intermediate-density neutral 
quark matter. Channels 2 and 3 (2 flavors, 1 or 2 colors) are unlikely to be of interest: 
they require inter-species pairing, but have gaps that are smaller than 1 MeV, so the 
same stress that disrupts the 2SC and CFL phases will even more readily disrupt 
these phases (channel 3 has been proposed to accompany the 2SC phase and allow 
pairing between blue up and down quarks [9, 98], but we have seen that the 2SC 
phase is unlikely to occur in neutral quark matter). The fourth channel (2 colors, 1 
flavor) and its 3-color cousin (the color-spin-locking phase) are the single-flavor color 
superconducting phases discussed earlier in this section. Channel 5 (1 color, 1 flavor) 
vanishes for light quarks, but it may allow pairing for an "orphaned" color of strange 
quark (i.e. a strange quark that is neglected by all other pairing processes). 

1.7 Applications 
1.7.1 Compact stars 

Our current understanding of the color superconducting state of quark matter leads us 
to believe that it may occur naturally within compact stars. The critical temperature 
below which quark matter is a color superconductor is estimated to be about 10 to 50 
MeV. In compact stars that are more than a few seconds old, the star temperature 
is less than this critical temperature and any quark matter that is present will be in 
a color superconducting state. It is therefore important to explore the astrophysical 
consequences of color superconductivity [19]. 

Much of the work on the consequences of quark matter within a compact star 
has focussed on the effects of quark matter on the equation of state, and hence on 
the mass-radius relationship [30]. The largest contributions to the pressure of quark 
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matter are a positive /i 4 contribution from the Fermi sea, and a negative bag constant 
B. As a Fermi surface effect, the effect of pairing is a contribution of order /i 2 A 2 . 
This is small compared to the two leading terms, so the conventional wisdom has 
been that superconductivity has a minor effect on the equation of state. Recently, 
however, it has been observed that if the bag constant is large enough so that nuclear 
matter and quark matter have comparable pressures at some density that occurs in 
compact stars, then there may be a large cancellation between the two leading terms 
and the Fermi surface term can have a large effect [33]. Therefore mass-radius curves 
can be sensitive to the presence of color superconductivity. 

A gravitational wave detector could yield insight into compact star interiors from 
observations of binary inspirals/mergers. In a hybrid star with a sharp interface 
between a nuclear mantle and a CFL color superconducting core, there is a large 
density discontinuity at the interface [39]. The two sharp density edges (at the core 
radius and at the star radius) could create features at two distinct time scales in the 
gravitational wave profile emitted during the inspiral and merger of a pair of compact 
stars of this type. The first feature would occur when the less dense nuclear mantles of 
the stars begin to deform each other; the second feature would occur only somewhat 
later when the denser cores begin to deform. 

The phase transition at which color superconductivity sets in as a hot proto- 
neutron star cools may yield a detectable signature in the neutrinos received from 
a supernova [28]. At the onset of color superconductivity the quark quasiparticles 
acquire gaps and the density of these quasiparticles is then suppressed by a Boltzmann 
factor exp(— A/T). As a result the mean free path for neutrino transport suddenly 
becomes very long. All of the neutrinos previously trapped in the core of the star 
are able to escape in a sudden burst that may be detectable as a bump in the time 
distribution of neutrinos arriving at an earth detector. 

Color superconductivity has a large effect on cooling and transport processes in 
quark matter [29, 99]. In quark matter, the neutrino emissivity is dominated by 
quasiquark modes that have energies smaller than the temperature T. These modes 
can rapidly radiate neutrinos by direct URCA reactions (d — > u + e + v, u — > d+e + + u, 
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etc.) which then dominate the cooling history of the star as a whole. In the CFL 
phase, all of the quarks have a gap A 3> T; the neutrino emissivity is suppressed by 
a Boltzmann factor exp(— A/T) and the CFL state does not contribute to cooling. 
In a compact star with a CFL core and a nuclear mantle, the cooling will occur only 
by neutrino emission from the mantle. 

This conclusion is revised for non-CFL phases of quark matter. Both the crys- 
talline color superconductor and the breached-pair color superconductor have gapless 
quasiquarks for momenta at the edges of "blocking regions" , as discussed in section 
1.4. These gapless modes could accommodate direct URCA reactions and conceivably 
dominate the entire cooling of the star [29, 99]. A similar effect could occur in the 
single-flavor spin-one color superconductor, which can have gapless quasiquarks at the 
poles or at the equator of the Fermi surface (in its non-color-spin- locked versions) [62]. 
Just how these special gapless modes could affect emissivity rates is unknown and 
is worthy of investigation. The crystalline and single-flavor phases also have collec- 
tive modes that will contribute to the heat capacity (phonons [48] and spin waves, 
respectively) . 

Recent work suggests that the observation of long-period (of order one year) pre- 
cession in isolated pulsars might constrain the possible behavior of magnetic fields 
in the core of a compact star [100]. Rotating compact stars with superfluid interiors 
will be threaded with a regular array of rotational vortices that are aligned along the 
axis of rotation. At the same time, if the core is a type II superconductor then it 
will also be threaded with an array of magnetic flux tubes that are aligned along the 
magnetic axis of the star. If the vortex and flux tube arrays coexist, they prevent 
any rotational precession because a precession would entangle the interwoven arrays. 

Remarkably, the observed precession therefore might rule out the standard model 
of a nuclear core containing coexisting neutron and proton superfluids, with the proton 
component forming a type II superconductor. But color superconducting interiors can 
accommodate the observed precession: magnetic flux tubes do not occur in either the 
CFL phase (which is not an electromagnetic superconductor [59]) or the single-flavor 
phase (which is a type I superconductor [96]). 
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Finally, in this thesis we wish to investigate the possibility that crystalline quark 
matter could be a locus for glitch phenomena in pulsars [41]. As the rotation of a 
pulsar gradually slows, the array of rotational vortices that fills the interior of the 
star should gradually spread apart; thus the star sheds its vortices and loses angular 
momentum. But if a crystalline phase occurs within the star, the rotational vortices 
may be pinned in place by features of the crystal structure. This impedes the smooth 
outward motion of the vortices. The vortex array could remain rigid until an accu- 
mulated stress exceeds the pinning force. Then, a macroscopic movement of vortices 
will occur, leading to an observed glitch in the rotational frequency of the pulsar. 
In chapter 5 we address the feasibility of this proposed glitch mechanism. With the 
crystal structure known, a calculation of the vortex pinning force can proceed. The 
first step is the explicit construction of a vortex state in the crystalline phase, and we 
discuss efforts in this direction. The task is a challenge by virtue of the interesting 
fact that the LOFF state is simultaneously a superfluid and a crystal. 

1.7.2 Atomic physics 

In section 1.5 we investigated crystalline color superconductivity with a two-flavor 
NJL model, i.e. a toy model with two species of fermion and a pointlike four-fermi 
interaction. This toy model may turn out to be a better model for the analysis of 
LOFF pairing in atomic systems. (There, the phenomenon could be called "crys- 
talline superfluidity".) Recently, ultracold gases of fermionic atoms such as lithium-6 
have been cooled down to the degenerate regime, with temperatures less than the 
Fermi temperature [101, 102, 103]. In these atomic systems, a magnetically-tunable 
Feshbach resonance can provide an attractive s-wave interaction between two differ- 
ent atomic hyperfine states [104]. This interaction is short-range but the scattering 
length can be quite long, so the system may be strongly-interacting. The attractive 
interaction renders the system unstable to BCS superfluidity below some critical tem- 
perature, and it seems possible to reach this temperature (perhaps by increasing the 
atom-atom interaction, thereby increasing T c , rather than by further reducing the 
temperature) [105, 103]. In these systems there really are only two species of fermion 
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(two different atomic hyperfine states) that pair with each other, whereas in QCD our 
model is a toy model for a system with nine quarks. The atomic interaction will be 
short-range and s-wave dominated, whereas in QCD it remains to be seen if this is a 
good approximation at accessible densities. In the atomic systems, experimentalists 
can control the densities of the two different atoms that pair, and in particular can 
tune their density difference. This means that experimentalists wishing to search for 
crystalline superfluidity have the ability to dial the most relevant control parame- 
ter [106, 107]. In QCD, in contrast, 5 fx is controlled by m^f/j,, meaning that it is 
up to nature whether, and if so at what depth in a compact star, crystalline color 
superconductivity occurs. 

Indirect observations of crystalline color superconductivity in the interior of a 
distant compact star are formidably difficult. But the atomic system provides a 
terrestrial setting in which the predictions of this thesis can be directly tested. In 
chapter 5 we will further discuss the prospect of creating and observing the crystalline 
state in an atomic gas. Because the spatial variation of the gap parameter can create a 
density modulation in the gas, it may be possible to literally see the crystal structure. 
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Chapter 2 

Crystalline Superconductivity: 
Single Plane Wave 

2.1 Overview 

In this chapter we study the simplest example of a crystalline color superconductor: 
a condensate that varies like a single plane wave in position space [41]. Equivalently, 
each Cooper pair in the condensate has the same total momentum 2q. We will use a 
variational method similar to that originally employed by Fulde and Ferrell [66] and 
described in more detail by Takada and Izuyama [89]. In section 2.2 we will describe 
the variational ansatz for the plane wave LOFF state. We note in particular that, 
unlike in the original LOFF context, there is pairing both in J = and J = 1 channels. 
In section 2.3, we derive the gap equation for the LOFF state for a model Hamiltonian 
in which the full QCD interaction is replaced by a four-fermion interaction with the 
quantum numbers of single gluon exchange. In section 2.4, we use the gap equation 
to evaluate the range of 5fi within which the LOFF state arises. We will see that 
at low Sfj, the translationally invariant BCS state, with gap Ao, is favored. At 5fi 1 
there is a first order transition to the LOFF paired state, which breaks translational 
symmetry. At 5fi 2 all pairing disappears, and translational symmetry is restored at 
a phase transition which is second order in mean field theory. In the weak-coupling 
limit, in which A <C //, we find values of 5^ and Sfi 2 which are in quantitative 
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agreement with those obtained by LOFF. This agreement occurs only because we 
have chosen an interaction which is neither attractive nor repulsive in the J = 1 
channel, making the J = 1 component of our LOFF condensate irrelevant in the 
gap equation. In section 2.5, we consider a more general Hamiltonian in which the 
couplings corresponding to electric and magnetic gluon exchange can be separately 
tuned. This leads to interactions in both J = and J = 1 channels, and we show 
how it affects the range of 8/jL within which the LOFF state arises. In section 2.6, we 
summarize the results for the plane wave crystal. 

2.2 The LOFF plane wave ansatz 

We begin our analysis of a LOFF state for quark matter by constructing a variational 
ansatz for the LOFF wavefunction. We consider Cooper pairs which consist of an up 
quark and a down quark with respective momenta 

k u = q + p, k d = q-p, (2.1) 

so that p identifies a particular quark pair, and every quark pair in the condensate 
has the same nonzero total momentum 2q. This arrangement is shown in Figure 2-1. 
The helicity and color structure are obtained by analogy with the "2SC" state as 
described in previous work [9, 10]: the quark pairs will be color 3 antitriplets, and in 
our ansatz we consider only pairing between quarks of the same helicity. 

With this in mind, here is a suitable trial wavefunction for the LOFF state with 
wavevector q [65, 66, 89]: 

l* q > = 4410), 

4 = II ( cos ^(P) + e a ^(p) sin^(p) oL,(q+p) a^(q-p)) 

peV,a,f3 

(2.2) 

x II a L a (q+p) x II a L/3(q-p)> 

4 = as above, L — > R, 
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Figure 2-1: The momenta k u and of the two members of a LOFF-state Cooper 
pair. We choose the vector q, common to all Cooper pairs, to coincide with the z-axis. 
The angles a u (p) and Qd(p) indicate the polar angles of k u and k^, respectively. 

where a, (3 are color indices, u, d and L, R are the usual flavor and helicity labels, 
and a) is the particle creation operator (for example, a} Lda creates a left-handed down 
quark with color a). The 0's and £'s are the variational parameters of our ansatz: 
they are to be chosen to minimize the free energy of the LOFF state, as described 
in the next section. The first product in equation (2.2) creates quark pairs within 
a restricted region V of the total phase space. This allowed "pairing region" will 
be discussed below. The next product fills a "blocking region" E u with unpaired up 
quarks: these are up quarks with momenta q+p for which there are no corresponding 
down quarks with momenta q — p. The final product fills the blocking region Bd 
with unpaired down quarks. The ansatz does not contain a term that would create 
antiparticle pairs: we have checked the effect of such a term and found that it has no 
qualitative effect on our results. 

To complete the specification of our ansatz we need to describe the allowed pairing 
and blocking regions in phase space. These regions are largely determined by Pauli 
blocking as a result of populated Fermi seas. In the absence of pairing interactions, 
the system is in the "normal" state and up and down quarks are distributed in Fermi 
seas with Fermi momenta = fi u and p d F = fid, respectively (recall that we consider 
massless quarks only, so the single particle energy of a quark with momentum k is 
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Figure 2-2: The LOFF phase space, as a function of p (Eq. (2.1)). We show the 
p y = plane, (a) The phase space in the limit of arbitrarily weak interactions. In 
the shaded blocking regions B u and Ba, no pairing is possible. In the inner unshaded 
region, an interaction can induce hole-hole pairs. In the outer unshaded region, an 
interaction can induce particle-particle pairs. The region V (Eq. (2.2)) is the whole 
unshaded area, (b) When the effects of interactions and the formation of the LOFF 
state are taken into account, the blocking regions shrink. The BCS singularity occurs 
on the dashed ellipse, defined by e u + = \i u + fid, where making a Cooper pair costs 
no free energy in the free case. 

e(k) = |k|). An up quark carries momentum k u = p + q; in p-space, therefore, the 
Fermi sea of up quarks corresponds to a sphere of radius fi u = fl — 5fi centered at — q. 
Similarly, a down quark carries momentum k^ = — p + q, giving a sphere in p-space 
of radius fid — fi+5fi centered at +q. The two offset spheres are shown in Figure 2-2a 
(we have drawn the case |q| > 5fi so that the two Fermi surfaces intersect in p-space). 
In the limit of arbitrarily weak interactions, the blocking region B u corresponds to the 
lower shaded area in the figure: pairing does not occur here since the region is inside 
the Fermi sea of up quarks, but outside the Fermi sea of down quarks. Similarly 
the upper shaded area is the blocking region Bd- The entire unshaded area is the 
pairing region V: it includes the region inside both spheres, where hole-hole pairing 
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can occur, and the region outside both spheres, where particle-particle pairing can 
occur. 

We can now explain how the LOFF wavefunction ansatz can describe the normal 
state with no condensate: we choose #l(p) = Or(p) = tt/2 for p inside both Fermi 
spheres, and otherwise all the #'s are zero. With this choice the first term in Eq. (2.2) 
fills that part of each Fermi sea corresponding to the inner unshaded region of Fig- 
ure 2-2a. The B u and Bd terms fill out the remainder of each Fermi sea to obtain the 
normal state. Note that in the absence of pairing, the normal state can be described 
with any choice of q. The most convenient choice is q = 0, in which case k u = = p, 
B u vanishes, and Bd is a spherical shell. Other choices of q correspond to choosing 
different origins of k n -space and k^-space, but in the absence of any interactions this 
has no consequence. Once we turn on interactions and allow pairing, we expect a 
particular |q| to be favored. 

The phase space picture changes slightly when pairing interactions are included: 
the blocking regions are smaller when a LOFF condensate is present, as indicated in 
Figure 2-2b. We will account for this effect in the next section. With smaller blocking 
regions, a larger portion of the phase space becomes available for LOFF pairing. Such 
pairing is guaranteed to be energetically favorable when it costs zero free energy to 
create an up quark and a down quark, since these quarks can then pair to obtain a 
negative interaction energy. The zero free energy condition is 



where e(k) is the single particle energy of a quark with momentum k. For massless 
quarks, we obtain |q + p| + |q — p| = 2/2, which describes an ellipsoidal surface in 
p-space. This surface is indicated by the ellipse shown in Figure 2-2b; notice that the 
ellipsoid and the two Fermi surfaces all intersect at a circle. 

If the interaction is weak, we expect LOFF pairing to be favored in a thin layer of 
phase space around this ellipsoid. This is manifest in the gap equation derived in the 
next section (Eq. (2.14)) in which, as in BCS theory, we find a divergent integrand 




(2.3) 
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on this ellipsoid in the absence of pairing. Pairing smooths the divergence. As the 
interaction gets stronger, the layer of favored pairing gets thicker. If there were no 
blocking regions, we could use the entire ellipsoid, just as BCS pairs condense over 
the entire spherical surface |p| = \i in the symmetric, 5fi = |q| = case. However, as 
shown in Figure 2-2b, the blocking regions exclude pairing over most of the ellipsoid, 
leaving a ribbon of unsuppressed LOFF pairing in the vicinity of the circle where the 
Fermi surfaces intersect. This agrees with our expectation for the particle distribution 
in the LOFF state: it is as in the normal state, except that there is a restricted region 
(around the aforementioned ribbon) where each quark in a pair can be near its Fermi 
surface and where pairing therefore occurs. 

Although the constant single-particle energy contours for noninteracting up and 
down quarks cross in p-space (see Figure 2-2a), we emphasize that the Fermi surfaces 
of up and down quarks do not cross in momentum (k u - and k^-) space. The p-space 
ribbon of unsuppressed pairing corresponds to unsuppressed pairing between up and 
down quarks with momenta around k-space ribbons near their respective (disjoint) 
Fermi surfaces. 

In the limit of arbitrarily weak interactions, the ribbon in momentum space along 
which pairing is unsuppressed shrinks, as the blocking regions grow to exclude all of 
the ellipsoid except the one- dimensional circle at which the two spheres in Figure 2-2 
intersect. This circle has insufficient phase space to lead to a singularity in the gap 
equation: the integrand is singular on this circle, but the integral does not diverge. 
Therefore, the LOFF state is not guaranteed to occur if one takes the weak coupling 
limit at fixed 6 fa In this respect, the LOFF state is like the BCS state at nonzero 
5 fa. for weak coupling, Ao — » and because the BCS state can only exist if it has 
Ao > y/25{i, it must vanish for couplings weaker than some threshold. We shall see, 
however, that at any fixed weak coupling, the LOFF state, like the BCS state, is 
guaranteed to occur at some 5 fa. the BCS state arises if Sfj, < 8 fa and the LOFF 
state arises if 5 fa < 5fi < 5 fa. 

One of the most striking features of the LOFF state is the spin structure of 
the condensate. The familiar "2SC" state pairs quarks of the same helicity and 
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opposite momentum, so the spins are antiparallel and the quarks are arranged in 
an antisymmetric combination to form spin singlet Cooper pairs. The LOFF state 
also pairs quarks of the same helicity, but now the quark momenta are no longer 
antiparallel, as can be seen from Figure 2-1. Therefore the LOFF Cooper pairs 
cannot be spin singlets: they are superpositions of both spin zero and spin one. This 
is revealed explicitly by evaluating the nonzero (ipip) expectation values in the LOFF 
state: 

(2.4) 

z(^ q |( C r 1 )^e a/3 3^ a (r)C , La 03 ^(r)|^ q ) = 2r^ r 

where i, j are flavor indices (1 = up, 2 = down), a, (3 are color indices, C = 27°7 2 , L = 
(1 — j5)/2 is the usual left-handed projection operator, and cr M „ = (i/2)^, 7„]. The 
constants and are left-handed J = and J = 1 condensates, respectively. T A 
and are defined analogously. The T's can be expressed in terms of the variational 
parameters of the LOFF wavefunction: 

T a = ^E sin ^(P) cos ^(P) e " L(P) sin( a " (p) + ad(p) )e-^) 

peV (2 5) 

pev 

Here V is the spatial volume of the system, a(p) are the polar angles of the quark 
momenta, as in Figure 2-1, and the dependence on the azimuthal angle <fi follows from 
our use of the spinor conventions described in Refs. [7, 9, 11]. The expressions for 

and r|| are the same as those in (2.5) except that </>(p) is replaced by 7r — <p(p). 
In Eq. (2.5) and throughout, {1/V) becomes J d 3 p/(2n) 3 in an infinite system. 

Once we have derived a gap equation by minimizing the free energy with respect 
to these variational parameters, we expect the condensates to be simply related to 
gap parameters occurring in the gap equation. We will see explicitly how Ta and Tb 
are determined in the next section. 

Notice that the condensates of Eq. (2.4) are plane waves in position space by virtue 
of the nonzero momentum 2q of a Cooper pair. Ta describes pairing which is antisym- 
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metric in color, spin, and flavor, while Tb describes pairing which is antisymmetric in 
color but symmetric in spin and flavor (in each case, Pauli statistics are obeyed). In 
the original LOFF condensate of electrons there can be no Tb, since electrons have 
no color or flavor, so that only the spin antisymmetric pairing is possible. 

The J = condensates (if) C Lip), (ipCRip) are Lorentz scalars (mixed under par- 
ity), while the J = 1 condensates (ipCLa 03 ip) , (ipCRa 03 ip) are 3- vectors (also mixed 
under parity) which point in the ^-direction, parallel to the spontaneously chosen 
direction q of the LOFF state. Because the ansatz contains a J = 1 component, 
it would be interesting to generalize it to include the possibility of LR pairing, in 
addition to LL and RR pairing. We discuss this further in Section 2.5. 



2.3 The gap equation and free energy 

Having presented a trial wavefunction for the LOFF state, we now proceed to min- 
imize the expectation value of the free energy (F) with respect to the variational 
parameters of the wavefunction (the #'s and £'s of equation (2.2)) to obtain a LOFF 
gap equation. The free energy is F = H — /j, u N u — /i^A^, where H is the Hamiltonian, 
and N u and A^ are the number operators for up and down quarks, respectively. We 
choose a model Hamiltonian which has a free quark term Hq and an interaction term 
Hj, and write the free energy as F = F + Hj, where F = H — fi u N u — /i^A^ is the 
free energy for noninteracting quarks. To describe the pairing interaction between 
quarks, we use an NJL model consisting of a four-fermion interaction with the color 
and flavor structure of one-gluon exchange: 

Hi = lj d 3 x [Ge^T^^T^) - Gm^T^^T^)} (2.6) 

where the T A are the color SU(3) generators, normalized so that tr(T A T B ) = 25 AB . 
Notice that we have relaxed some constraints on the spin structure of one-gluon 
exchange: we allow for the possibility of independent couplings Ge and Gm for 
electric and magnetic gluons, respectively. This spoils Lorentz boost invariance but 
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there is no reason to insist on boost invariance in a finite-density system. Indeed, 
in high density quark matter we expect screening of electric gluons but only Landau 
damping of magnetic gluons, and we might choose to model these effects by setting 
Ge <C Gm- We postpone a discussion of these issues and their implications for the 
LOFF state until Section 2.5. For now, we restrict ourselves to the case of Lorentz 
invariant single gluon exchange, by letting Ge = Gm = G > 0. 

We need to evaluate the expectation value of F in the LOFF state to obtain an 
expression for the free energy of the system in terms of the variational parameters of 
the ansatz. The noninteracting part of the free energy is simply 

(F > = 2 ^ +p\-^) + J2 2 (i<i - pi - ^) 

+ ^ 2(|q + p| + |q - p| - /i u - fi d ) sin 2 9 L (p) (2.7) 
+ (same, with L —> R). 

The first and second terms represent the contributions of the unpaired left-handed up 
and down quarks, respectively. The third term gives the (noninteracting) free energy 
of the left-handed quark pairs. The three terms are all repeated with L replaced by R 
to include the free energy for the right-handed quarks. The factors of two in equation 
(2.7) appear because there are two quark colors ("red" and "green") involved in the 
the condensate. The "blue" quarks do not participate in the pairing interaction and 
instead behave as free particles: the blue up and down quarks fill Fermi seas with 
Fermi momenta p F = fi u and p F = fid, respectively. Below, we will want to compare 
the free energy of the LOFF, BCS and normal states. Since at any given \l u and \x d 
the free energy of the spectator quarks is the same in all three states, we can neglect 
these blue quarks in the remainder of our analysis even though they do contribute to 
the total free energy. 

The expectation value of Hi gives the total binding energy of the pairing interac- 
tion: 

(H I ) = -\GV{\T A f + \T R A f) (2.8) 
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where the T^'s are the J = LOFF condensates defined in equations (2.5). These 
condensates are simply related to J = LOFF gap parameters defined as 




(2.9) 



The gap parameters A^ correspond to 1PI Green's functions and are the quantities 
which will appear in the quasiparticle dispersion relations and for which we will 
derive the self-consistency conditions conventionally called gap equations. We see 
from Eq. (2.8) that with G > the interaction is attractive in the J = channel and 
is neither attractive nor repulsive in the J = 1 channel. 

Our ansatz breaks rotational invariance, so once J = pairing occurs (r^ ^ 0) 
we expect that there will also be J = 1 pairing (r^ ^ 0). As we have seen, this 
arises even in the absence of any interaction in the J = 1 channel as a consequence 
of the fact that the momenta of two quarks in a Cooper pair are not anti-parallel 
if q 7^ 0. Because (H) is independent of Tb, the quasiparticle dispersion relations 
must also be independent of Tb- That is, the J = 1 gap parameter must vanish: 
Ab = 0. In Section 2.5, we shall see by direct calculation that A# is proportional 
to (Ge — Gm)Tb- In the present analysis with Ge = Gm, therefore, A# = while 



The £'s are chosen to cancel the azimuthal phases 0(p) in equations (2.5). By this 
choice we obtain maximum coherence in the sums over p, giving the largest possible 
magnitudes for the condensates and gap parameters. We have 



where y?£ and (fR are arbitrary p-independent angles. These constant phases do not 
affect the free energy — they correspond to the Goldstone bosons for the broken 
left-handed and right-handed baryon number symmetries — and are therefore not 
constrained by the variational procedure. For convenience, we set ipi = = and 
obtain condensates and gap parameters that are purely real. 

The relative phase tp^ — ifR determines how the LOFF condensate transforms 





(2.10) 
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under a parity transformation. Its value determines whether the J = condensate is 
scalar, pseudoscalar, or an arbitrary combination of the two and whether the J = 1 
condensate is vector, pseudovector, or an arbitrary combination. Because single gluon 
exchange cannot change the handedness of a massless quark, the left- and right- 
handed condensates in the LOFF phase are not coupled in the free energy of Eq. (2.8.) 
Our choice of interaction Hamiltonian therefore allows an arbitrary choice of (pi — ip^. 
A global U(1)a transformation changes Pl — Pr, and indeed this is a symmetry of our 
toy model. If we included [7(1 ^-breaking interactions in our Hamiltonian, to obtain 
a more complete description of QCD, we would find that the free energy depends on 
<Pl — <Pr, and thus selects a preferred value. For example, had we taken Hi to be 
the two-flavor instanton interaction as in Ref. [9, 10], the interaction energy would 
appear as Y L *T R + r^T^* instead of as in (2.8). This would enforce a fixed phase 
relation ipi — tp^ = 0, favoring condensates which are parity conserving [9, 10]. 

We now apply the variational method to determine the angles 9(p) in our trial 
wavefunction, by requiring that the free energy is minimized: d(F) /d9(p) = 0. This 
is complicated by the fact that the pairing region V and the blocking regions B u and 
Bd are themselves implicitly dependent on the 9 angles: these angles determine the 
extent of the LOFF pairing, and the phase space regions V, B u and Bd change when 
a condensate is present, as mentioned in Section 2.2. For now we simply ignore any 
6*-dependence of the phase space regions; our result will nevertheless turn out to be 
correct. Everything is the same for left and right condensates so we hereafter drop 
the L and R labels. Upon variation with respect to 9(p), we obtain 

„, , 2A A sm(f3 A (p)/2) 
tan2#(p) = , A , VMA ' 2.11 

|q + p| + |q-p| - a*u- v-d 

where (3 a (p) = &u(p) + «d(p) is the angle between the two quark momenta in a 
LOFF pair, as shown in Figure 2-1. Notice that the denominator on the right hand 
side of the above expression vanishes along the ellipsoidal surface of optimal LOFF 
pairing described in Section 2.2. When q = 0, the quark momenta are antiparallel so 
Pa(p) = 7T and Eq. (2.11) reduces to the simple BCS result: tan 29 = Aa/(\p\ — ft)- 
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With the 9 angles now expressed in terms of a gap parameter A^, we turn to 
the LOFF quasiparticle dispersion relations. They can be obtained by taking the 
absolute value of the expressions 



8fi + |(|q+p| - |q-p|) 



+ ^(h+p\ + h-p\-m 2 +^A 



isin 2 (|/? A (p)) 



(2.12) 



E 2 (p) = -5fx - |(|q + p| - |q - p|) 



isin 2 (i/3 A (p)) , 



whose meaning we now describe. For regions of p-space which are well outside both 
Fermi surfaces, E\ [E 2 ) is the free energy cost of removing a LOFF pair and adding 
an up quark with momentum q + p (a down quark with momentum q — p). For 
regions of p-space which are well inside both Fermi surfaces, E\ (E 2 ) is the free 
energy cost of removing a LOFF hole pair and adding a down hole with momentum 
q — p (an up hole with momentum q + p). Where the Fermi surfaces cross in p-space 
and pairing is maximal, both quasiparticles are equal superpositions of up and down. 
In the region of p-space which is well inside the up Fermi surface but well outside 
the down Fermi surface, E\ is negative, corresponding to a domain in which it is 
energetically favorable to have an unpaired up quark with momentum q + p rather 
than a (q + p, q — p) quark pair. Similarly, E 2 is negative where it is favorable to have 
an unpaired down quark with momentum q — p rather than a LOFF pair. Equations 
(2.12) allow us to finally complete our description of the LOFF phase by specifying 
the definitions of the phase space regions V, B u and Bd- The blocking region B u is 
the region where Ei(p) is negative, and unpaired up quarks are favored over LOFF 
pairs. Similarly Bd is the region where -E^p) is negative. The regions E\ < and 
E 2 < are shown as the shaded areas in Figure 2-2a for = 0, and in Figure 2-2b 
for A^ 7^ 0. LOFF pairing occurs in the region where E\ and E 2 are both positive: 



V = {p|£i(p) > and E 2 (p) > 0} 



(2.13) 
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corresponding to the entire unshaded regions of Figure 2-2. The actual quasiparticle 
dispersion functions are |i?i(p)| and l-E^p)): they are nonnegative everywhere, since 
they represent energies of perturbations of the LOFF state which is the presumed 
ground state of the system. 1 In the blocking regions, elementary excitations are 
created by replacing an unpaired quark with a quark pair, and vice versa in the 
pairing region. When q = 0, Eqs. (2.12) reduce to the more familiar BCS result: 
E {m (p) = ±5fi + V(|p|-£) 2 + Ai. 

With the boundaries of the blocking regions specified, one can verify by explicit 
calculation that the variation of these boundaries upon variation of the #'s does not 
change the free energy. This can be understood as follows. Notice that because we 
can create zero-energy quasiparticles on the boundaries of the blocking regions, there 
is no actual energy gap in the excitation spectrum of the LOFF state. The change 
in (F) due to variation of the boundaries of the blocking regions is zero because this 
variation simply creates zero-free-energy quasiparticles on these boundaries. This 
justifies our neglect of the ^-dependence of the phase space regions in the derivation 
of Eq. (2.11). 

Substituting the expression (2.11) for the 9 angles into the expression (2.5) for 
the Ta condensate, and using the relation A^ = GTa, we obtain a self-consistency 
equation for the gap parameter A^: 

1 = 2Gy- 2sin 2 (|/? A (p)) (2u) 

V Per J(\q + p| + |q - p| - 2/x)* + AA\ sin 2 (|^(p)) ' 



This can be compared to the BCS gap equation, obtained upon setting q = and 
eliminating the blocking regions: 

i = ggy-— 1 (2.i5) 

Note that in the LOFF gap equation (2.14), the gap parameter appears on the right 

1 Since the LOFF condensate contains pairs with momentum 2q, the momentum of its quasipar- 
ticle excitations is only defined modulo 2q. The momentum, modulo 2q, of a quasiparticle of energy 
|-Ei(p)| is p mod 2q. 



61 



hand side both explicitly in the denominator and also implicitly in the definition of 
the pairing region V, as given in (2.13). This means that if the q — > limit is taken 
at fixed the LOFF gap equation will only become the BCS gap equation if the 
blocking regions vanish in this limit. This happens if, as q — > 0, A^ tends to a 
limiting value which is greater than 5fi. For A^ < 5fi, the blocking region Bd does 
not disappear but instead becomes a spherical shell, as we can see by taking q — > 
in figure 2-2. Pairing is blocked in the region 



ji - yjstf - A \ < |p| < ft + yj 8tf - A\. (2.16) 

This is precisely the "breached pair" color superconductor discussed in chapter 1 [75, 
76]; the breach is just the blocking region 13 a- This second solution to the q = gap 
equation was first discovered by Sarma [74]. In the present calculation such states 
will always have higher free energy than both the unpaired state (A^ = 0) and the 
BCS state obtained by solving the gap equation (2.15) with no blocking regions.. But 
there may exist other scenarios in which the Sarma solution is favored: when one 
species is very heavy and the other very light, for example [76]. 

In the next section we will solve the LOFF gap equation (2.14) and determine the 
circumstances in which the LOFF state is the true ground state of the system. Once 
we have obtained a solution to the gap equation (2.14) for A^, the condensates are 
given by Ta = A^/G and 



^ B - 2 y 2A A sin(^(p))sin(^(p)) 

/ P ev J(\q + p| + |q - p| - 2/x)* + 4A^ sin 2 (|/5 A (p)) 



where Pb(p) = a«(p) — <^d(p)- (See Figure 2-1.) We now see explicitly that if 
the interaction is attractive in the J = channel, creating a nonzero Ta and A^, 
a nonzero J = 1 condensate r# is induced regardless of the fact that there is no 
interaction in the J = 1 channel. As a check, note that if q = 0, sin(|/?A(p)) = 1 
and sin(|/3e(p)) is given by the cosine of the polar angle of p. The right hand side of 
(2.17) therefore vanishes upon integration, and Tb vanishes when q = as it should. 
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It is now apparent that two features contribute to a nonzero Tg- The first is that the 
momenta in a quark pair are not antiparallel, which leads to the factors of sin(|/3^(p)) 
in Eq. (2.17). The second is that the pairing region is anisotropic, since if it were 
not the factor of sin(|/5s(p)) would ensure that the right-hand side of (2.17) vanishes 
upon integration. 

As written, the gap equations (2.14) and (2.15) are ultraviolet divergent. In QCD, 
of course, asymptotic freedom implies that the interaction between quarks decreases 
at large momentum transfer and we have not yet represented this fact in our toy 
model. In previous work [9, 11, 23], we chose to mimic the effects of asymptotic 
freedom (and to render the right hand side of the gap equation finite) by introducing 
a form factor associated with each fermion leg in the four-fermion interaction. This is 
not a good strategy when q^O. The two incident quarks carry momenta q + p and 
q — p while the outgoing quarks carry momenta q + p' and q — p'. Were we to cut 
off these four momenta with form factors on each leg, we would have a cutoff which 
depends explicitly on q. This is not a good representation of what happens in full 
QCD, in which the condition for when the interaction becomes weak is determined 
by the momentum p — p' transferred through the gluon and has nothing to do with 
q. For simplicity, we choose to introduce a hard cutoff in our NJL model, rather than 
a smooth form factor, and choose simply to cut off the momentum p. This is not 
equivalent to cutting off the momentum transfer, but has the desired feature of being 
a q-independent cutoff. That is, we limit the integration region to |p| < A in the BCS 
gap equation (2.15) and to {p G V and |p| < A} in the LOFF gap equation (2.14). In 
the BCS case, this criterion is equivalent to cutting off the momentum of each fermion 
leg. In the LOFF case, it is not equivalent and is more appropriate. The choice we 
have made is not the only q-independent cutoff one might try. For example, we have 
also obtained results upon cutting off momenta outside a large ellipsoid in p-space, 
confocal with the centers of the two Fermi spheres in Figure 2-2, but have found that 
this makes little difference relative to the simpler choice of the large sphere |p| < A. 
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2.4 Results 



We solve the gap equation (2.14) numerically (and analytically in the limit A^ <C 
5fi, q, A ) and calculate the LOFF state free energy as a function of <5/i and q, for 
given coupling G, average chemical potential fl, and cutoff A. We vary q to minimize 
the LOFF free energy, and compare it with that for the standard BCS pairing (2.15) 
to see which is favored. In this way we can map out the phase diagram for the three 
phases of pairing between the two species of quark: BCS, LOFF, and unpaired. 

Note that the solution to the gap equation, the LOFF gap parameter A^, is not a 
gap in the spectrum of excitations. The quasiparticle dispersion relations (2.12) vary 
with the direction of the momentum, yielding gaps that vary from zero (for momenta 
on the edge of the blocking regions in phase space) up to a maximum of A^. 

We will first discuss the range of 5fi in which there exists a LOFF state as a local 
energy minimum. Later we will go on to study the competition between LOFF and 
BCS, and see in what range of <5/i the LOFF state is the global minimum. We expect 
the BCS state to be preferred when the mismatch S/j, between the Fermi energies of 
the two species is small. When the mismatch is comparable to the BCS gap (5fi ~ Ao) 
we expect a transition to LOFF, and at larger 5/x we expect all pairing to cease. These 
expectations are largely borne out. 

In general we fix A = 1 GeV and p, — 0.4 GeV, and study different coupling 
strengths G which we parameterize by the physical quantity Ao, the BCS gap of 
Eq. (2.15) which increases monotonically with increasing G. When we wish to study 
the dependence on the cutoff, we vary A while at the same time varying the coupling G 
such that Ao is kept fixed. (This is in the same spirit as using a renormalization con- 
dition on a physical quantity — Ao — to fix the "bare" coupling — G.) We expect that 
the relation between other physical quantities and A will be reasonably insensitive 
to variation of the cutoff A. 

We wish to determine 5/i 2 , the boundary separating the LOFF phase and the 
normal phase. The LOFF to unpaired phase transition is second order, so it occurs 
where the solution A^ to the LOFF gap equation (2.14) is zero. Setting A^ = in 
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Figure 2-3: (a) The zero-gap curves for the LOFF state. To the right of a solid curve, 
there is no solution to the LOFF gap equation, to the left of the curve there is a 
solution, and on the curve the gap parameter is zero. The three curves are (from 
strongest to weakest coupling): A = 0.1,0.04,0.01 GeV. The region q < 5/i is 
complicated to describe [66], and solutions found in this region never give the lowest 
free energy state at a given 5fi. (b) Here, we choose Ao = 0.04 GeV and focus on 
the region near <5/i2, the maximum value of 5fi at which the LOFF state exists. The 
dashed curve shows the value of |q| which minimizes the free energy of the LOFF 
state at a given 5fi. 5/ii, discussed below, is also indicated. 

the gap equation (2.15) yields an analytical expression relating Sfi and q, for any given 
G and A. In Figure 2-3a we show the = curve for three couplings corresponding 
to A = 0.1 GeV (strong coupling), A = 0.04 GeV and A = 0.01 GeV (weak 
coupling). We have only drawn the zero-gap curve in the region where q > 5fi. We 
expect this to be the region of interest for LOFF pairing because when q > 5n the 
two spheres of Figure 2-2 do in fact intersect. We have verified that, as described in 
some detail in Ref. [66], there are regions of Figure 2-3a with q < 8/j, within which 
the LOFF gap equation (2.14) has (one or even two) nonzero solutions, but these 
solutions all correspond to phases whose free energy is either greater than that of 
the normal phase or greater than that of the BCS phase or both. Figure 2-3 shows 
that for a given coupling strength, parameterized by A , there is a maximum Sfi for 
which the LOFF state exists: we call it 5^2- For 5fi > <5/i2, the mismatch of chemical 
potentials is too great for the LOFF phase to exist. 

We see from Figure 2-3a that as the coupling gets weaker, <5/i 2 /Ao gets gradually 
larger. (Of course, <5/i 2 itself gets smaller: the quantities plotted are £/z/A and qj A .) 
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Note that in the A — >• limit, the zero gap curve is essentially that shown in the 
figure for A = 0.01 GeV, in agreement with the curve obtained at weak coupling by 
Fulde and Ferrell [66]. The fact that this curve ceases to move in the Ao — > limit 
means that 8/X2 — ► while 8/x 2 /Ao — > const in this limit. 

For 8fi — > Sfi2 from below, we see from Figure 2-3 that there is a solution to the 
LOFF gap equation only at a single value of q. For example, at A = 0.04 GeV we 
find q = 0.880A = 1.183 8fx 2 at 6fi 2 = 0.744A . (In agreement with Refs. [65, 66], 
in the weak coupling limit we find q = 0.906A = 1.20 8 fx 2 at 8/X2 = 0.754A .) For 
any value of 8 fx < 8fx 2 , solutions to the LOFF gap equation exist for a range of |q|. 
We must now find the value of |q| for which the free energy of the LOFF state is 
minimized. We obtain the free energy of the LOFF state at a point in Figure 2-3 by 
first solving the gap equation (2.14) numerically to obtain A^, and then using (2.9) 
and (2.11) to evaluate (Fq + Hi) given in (2.7) and (2.8). For each value of 8/x < 8/jL 2 
we can now determine which choice of q yields the lowest free energy. The resulting 
"best-g curve" curve is shown in Figure 2-3b for Ao = 0.04 GeV. 2 

Finally, for each point on the best-g curve we ask whether the LOFF free energy at 
that 8 [a and (best) q is more or less than the free energy of the BCS state at the same 
8 fx. In this way, we find 8/xi at which a first order phase transition between the LOFF 
and BCS states occurs. In Figure 2-4 we show the competition between the BCS and 
LOFF states as a function of the Fermi surface mismatch 8/x, for a fixed coupling 
corresponding to A = 40 MeV. The LOFF state exists for 8fx < 8fi 2 = 0.744A . At 
each 8 fx < 8 /X2, we plot the gap parameter and free energy characterizing the LOFF 
state with the best q for that 8/x. Although the BCS gap A is larger than the LOFF 

gap A a, as 8 fx increases we see from Eq. (1.16) that the BCS state pays a steadily 

2 As a check on our determination of the best q, we have confirmed that the total momentum of 
the LOFF state with the best q is zero, as must be the case for the ground state of the system at a 
given S/x (by a theorem attributed to Bloch [108]). This is a powerful check, because it requires the 
net momentum of the unpaired quarks in the blocking regions (which is in the negative z direction; 
see Figure 2-2) to be cancelled by the net momentum carried by the LOFF condensates. When, in 
future work, our ansatz is extended to describe a LOFF crystal rather than a single plane wave, this 
check will no longer be powerful. Once we go from T ~ exp(2iq • r) to T ~ cos(2q • r) or to a more 
involved crystalline pattern, the total momentum of the condensates and of the unpaired quarks will 
each be zero. 
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Figure 2-4: LOFF and BCS gaps and free energies as a function of 5fi, with coupling 
chosen so that Ao = 40 MeV and with fx = 0.4 GeV, A = 1 GeV. Free energies 
are measured relative to the normal state. At each Sfi we have varied q to find the 
best LOFF state. The vertical dashed line marks 8ll = 5fi 1 , the value of 8ll above 
which the LOFF state has lower free energy than BCS. The expanded inset (wherein 
s = 10~ 7 GeV 4 ) focuses on the region 5li 1 < 5li < 5[i 2 where the LOFF state has the 
lowest free energy. This figure should be compared with the sketch in figure 1-4. 



increasing free-energetic price for maintaining p u F = pp, whereas the LOFF state pays 
no such price. We now see that the LOFF state has lower free energy than the BCS 
state for 5fi > 5lli, in this case 8li\ = 0.7104A . At 5ll = 5lli, the gap parameter is 
A^ = 0.0078 GeV = 0.195A . (Had we calculated <5/i 1 by comparing the BCS free 
energy with that of the unpaired state instead of with that of the LOFF state, we 
would have obtained 5li\ = 0.71lAo. As the inset to Figure 2-4 confirms, the BCS 
free energy varies so rapidly that this makes an almost imperceptible difference. In 
later figures, we therefore obtain <5/ii via the simpler route of comparing BCS vs. 
normal.) At the coupling corresponding to A = 40 MeV, we have found that the 
LOFF state is favored over both the BCS state and the normal state in a "LOFF 
window" 0.710 < 5/i/A < 0.744. 

With solutions to the gap equation in hand, we can obtain the J = condensate 
Ta = GAa and the J = 1 condensate Tb given in Eq. (2.17). In Figure 2-5, we show 
both condensates within the LOFF window 5ll\ < 5li < 5ii2- We see first of all that 
fg ^ 0, as advertised. For the choice of parameters in Figs. 2-4 and 2-5 we find 
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Figure 2-5: The two LOFF condensates (J = 0) and T B (J = 1) for the same 
choice of parameters as in Figure 2-4. We focus on the region 5fi\ < 5fi < 5fi2- For 
reference, in the BCS phase T A = A /G = 0.00583 GeV 3 and T B = 0. 

Fb/^a essentially constant over the whole LOFF window, varying from 0.121 at 
to 0.133 at 5/i 2 . Increasing A tends to increase T b /Ta, as does decreasing A. Second 
of all, we see that the phase transition at 5fi = 5fi2, between the LOFF and normal 
phases, is second order in the mean-field approximation we employ throughout. 

It is interesting to explore how the width of the LOFF window depends on the 
strength of the coupling, and to confirm that it is insensitive to the cutoff. We do 
this in Figure 2-6, where we plot 6^2/ A (solid lines) and Sfi\/A (dashed lines). 
The LOFF state is favored for 5fii/A < 5fi/A < <5/i 2 /A , i.e. between the solid 
and dashed curves in Figure 2-6. In the weak coupling limit, the LOFF window 
tends to 0.707 < 5^/A < 0.754 and A A at 5^ tends to 0.23A , as in Refs. [65, 66]. 
Note that if one takes the weak-coupling limit Ao — > at fixed 5fi, neither BCS nor 
LOFF pairing survives because <5/i/A — > oo. However, for any arbitrarily small but 
nonzero coupling, the LOFF phase is favored within a range of 5fi. Figure 2-6 thus 
demonstrates that in an analysis of the LOFF state in the weak-coupling limit, it is 
convenient to keep S^/Aq fixed while taking Ao — > 0. We see from Figure 2-6 that 
strong coupling helps the BCS state more than it helps the LOFF state. When the 
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Figure 2-6: The interval of 5fi within which the LOFF state occurs, as a function of 
the coupling (parameterized as usual by the BCS gap A ). Below the solid line, there 
is a LOFF state. Below the dashed line, the BCS state is favored. The different lines of 
each type correspond to different cutoffs A = 0.8 GeV to 1.6 GeV. b\i\j A and 5fi 2 / Ao 
show little cutoff-dependence, and the cutoff-dependence disappears completely as 
A ,<5/i->0. 

coupling gets strong enough, there is no longer any window of Fermi surface mismatch 
8/jL in which the LOFF state occurs: the BCS state is always preferred. 

The different lines of each type in Figure 2-6 are for different cutoffs and show 
that there is in fact little sensitivity to the cutoff. The A dependence of 5/xi/Ao and 
S/j,2/A is mild for all values of A which are of interest, and is weakest for A — > 0. 
This is because in that limit pairing can only occur very close to the unblocked ribbon 
of the ellipsoid of Fig. 2b, along which the integrand in the gap equation is singular 
and pairing is allowed. Thus most of the pairing region V, and in particular the 
region near A, become irrelevant in this limit. 

The one physical quantity which we have explored which does turn out to depend 
qualitatively on A is the ratio T g /Fa- Those quarks with momenta as large as A which 
pair have momenta which are almost antiparallel, and so contribute much less to 
than to T^. For this reason, the ratio Tb/^a is sensitive to the number of Cooper 
pairs formed at very large p, and hence to the choice of A. As discussed above, 
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pairing far from the favored ribbon in phase space becomes irrelevant for A — > 0, 
and indeed in this limit we find that the A dependence of Tb/^a decreases. However, 
for Ao = 40 MeV we find that changing A from 1.2 GeV to 0.8 GeV increases V b/^a 
by more than 50%. 

We chose to show results for A = 40 MeV in Figure 2-4 because with this choice, 
the LOFF window occurs at values of 8 fx that are reasonable for quark matter in the 
interiors of compact stars. Consider fx = 400 MeV and m s = 300 MeV (recall that 
m s is a density-dependent effective mass which is significantly larger than the current 
quark mass). Substituting these numbers in equations (1.9) describing unpaired neu- 
tral quark matter, we find that the baryon number density is 4 times nuclear matter 
density and 8 fx = \{fXd — (J>u) = 28 MeV. (Had we chosen m s = 200 MeV, we would 
have obtained 8fx — 12 MeV and a baryon number density that is 5 times nuclear 
matter density.) Of course, neither 8 fx nor the value of Ao are accurately known for 
the quark matter which may exist within a compact star. Still, it seems possible that 
their ratio could be appropriate for the quark matter to be in the LOFF phase. If 
there is a range of radii within a compact star in which quark matter occurs with 
5 Hi < 8fx < 6/i<2, this quark matter will be a crystalline color superconductor. 

In Figure 2-4, the LOFF gap parameter A^ is 7.8 MeV at 8 fx = Sfii. It remains 
larger than typical neutron star temperatures T ns ~ 1 keV until very close to 8 fx = S/j.2- 
Similarly, the LOFF free energy, which is 4.8 x 10" 8 GeV 4 = 4.8 x (10 MeV) 4 at 
8 fx = 5fxi, is much larger than T 4 S throughout the LOFF window except very close 
to 8 /x = 8fX2- The LOFF gap and free energy are likely to be much larger for a 
condensate of multiple plane waves, as we will see in the next chapter. Furthermore, 
we will see in chapter 5 that the free energy of the LOFF state is of the right order 
to lead to interesting glitch phenomena. 

2.5 More general Hamiltonian and ansatz 

In Section 2.3, we introduced the four-fermion interaction Hamiltonian Hi of Eq. (2.6) 
with independent couplings Ge and Gm for the interactions which model the exchange 
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of electric and magnetic gluons. It proves convenient to use the linear combinations 



Ga — \{Ge + 3Gm) 
Gb = \{Ge — Gm) ■ 



(2.18) 



of the coupling constants in terms of which the expectation value of Hj in the LOFF 
state (2.2) becomes 



Thus, a positive coupling Ga describes an attractive interaction which induces a J = 
condensate Fa- As we have seen, in the LOFF state this is necessarily accompanied 
by a J = 1 condensate Tb- In our analysis to this point, we have set Ga = G > 
and Gb = 0. We now discuss the general case, in which Gb ^ 0. 

Before beginning, let us consider how to choose Gb/Ga in order for our model 
Hamiltonian to be a reasonable toy model for QCD at nonzero baryon density. At zero 
density, of course, Lorentz invariance requires Gb = 0. At high densities, on the other 
hand, electric gluons are screened while static magnetic gluons are not. (Magnetic 
gluons with nonzero frequency are damped.) We now know [12] that at asymptotically 
high densities it is in fact the exchange of magnetic gluons which dominates the pairing 
interaction. This suggests the choice Ge = 0, corresponding to Gb/Ga = — 1/3. At 
the accessible densities of interest to us, it is presumably not appropriate to neglect 
Ge completely. Note also that the four-fermion interaction induced by instantons 
in QCD only yields interactions in flavor-antisymmetric channels. It results in an 
attractive interaction in the J = channel and no interaction in the J = 1 channel. 
Thus, although the instanton interaction cannot be written in the form (2.6), for our 
purposes it can be thought of as adding a contribution to Ga, but none to Gb- Hence 
our model is likely to best represent high density QCD for a ratio of couplings lying 
somewhere in the range 



(Hi) 




rir + ir^D-iG^Or^ + irfi 2 ) 



(2.19) 




(2.20) 
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We plot our results over a wider range of couplings below. 

Once Gb ^ and there is an interaction in the J = 1 channel, we expect, in 
addition to the J = 1 condensate Tg, a J = 1 gap parameter Ag. The quasiparticle 
dispersion relations are then determined by A^ and As, which are defined as 

AA = GATA (2.21) 
Ab = Gb^b ■ 

Following through the variational calculation as in Section 2.3 leads to the coupled 
gap equations: 

2Ga sr-^ 2Sa(A a Sa + A B S B ) 



v tt v Vdq + pI + |q - pI - 2^) 2 + ^aS a + &bS b ) 2 

2Gb sr-^ 2Sb(AaSa + A b Sb) 

V t£ V(\V + P\ + |q - Pi - 2/i) 2 + 4(A A S A + Kb~Sb? (2.22) 

S A = sin(|^(p)) 

S B = sin(| / 9 B (p)) 

with /3a(p) = ck«(p) + ttd(p), /9s (p) = Q! u (p) — a^(p) defined in terms of the an- 
gles in Figure 2-1. The pairing region V is still defined by (2.13) but with new 
quasiparticle dispersion relations obtained from Eqs. (2.12) with A\S\ replaced by 
(A A S A + A B S B ) 2 . 

For Gb = 0, the coupled equations (2.22) reduce to Eqs. (2.14) and (2.17). Note 
that if, instead, G B > and Ga = 0, we find an attractive interaction in the J = 1 
channel in Eq. (2.18) and no interaction in the J = channel. Analysis of Eqs. (2.22) 
in this case yields a nonzero value of A#, while A^ = even though Ta ^ 0. The 
geometry of the LOFF pairs requires Ta ^ when Tb ^ 0- 

Rather than describing how every Figure in Section 2.4 changes when Gb ^ 0, we 
choose to focus on the question of how the interval of 5[i within which the LOFF state 
occurs (the LOFF window) chang function of Gb/Ga- To further simplify the 

presentation, we specialize to the weak-coupling limit in which Ao — > 0. This means 
that, as in Figure 2-6, the LOFF window is independent of the cutoff A. 
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Figure 2-7: The interval of 8fx in which the LOFF state is favored at weak coupling, 
as a function of the ratio of couplings Gb/Ga- Below the solid line, there is a LOFF 
state. Below the dashed line, the ordinary BCS state is favored. Gb = corresponds 
to the Lorentz-invariant interaction with Ge = Gm- QCD at high density is likely 
best described by a coupling in the range — | < Gb/Ga < 0. 

We show the dependence of the LOFF window on Gb/ Ga in Figure 2-7. The lower 
boundary 8 fx = 8fi l is, as in Section 2.4, the same (up to a very small correction) as 
the 8fi at which the BCS and normal states have equal free energies. We find the 
upper boundary 8fx = 8fi 2 by first dividing Eqs. (2.22) by A^ and then looking for a 
value of 8 fx at which Aa — > and Ab — > but A^/Ab remains nonzero. As before, 
this defines a zero-gap curve, and 8/X2 is the maximum value of Six reached by this 
curve. 

We find that the lower boundary 8fx 1 is completely unaffected by the value of 
Gb, since the BCS state is purely J = 0. So in the weak-coupling limit we obtain 
the result of Section 2.4, <5/i 1 /A = 0.707, independent of Gb/Ga- In contrast, 8/x 2 , 
the upper boundary of the LOFF window, increases with increasing Gb- This is 
understandable: the LOFF state always produces a J = 1 condensate, so we expect 
it to be fortified by Gb > and penalized by Gb < 0. There is no analogue of 
this behavior in an electron superconductor [65, 66], where there can be no J = 1 
condensate. Our J = 1 condensate affects the gap equation and free energy only if 
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Gb ^ 0; for this reason, our weak coupling results are in agreement with those of 
LOFF [65, 66] only if Gb = 0, as in Section 2.4. The effect of a coupling Gb in the 
physically interesting range (2.20) is to reduce the LOFF window, but only slightly. 

In the vicinity of 5fi 2 we should evaluate the competition between LOFF pairing 
and the formation of single-flavor (uu) and (dd) condensates. These are spatially- 
uniform spin-one condensates, as introduced in section 1.6 and described in more 
detail in the NJL survey of chapter 4. We defer a complete discussion of these con- 
densates until chapter 4, but here review some aspects that are relevant in the current 
context. Considering only the up quarks, the single-flavor spin-one condensates can 
involve LL/RR pairing, with a condensate structure (uCaosu), and LR pairing, with 
a condensate structure (uC^u). In the LL/RR channel, the Ga term in our model 
Hamiltonian gives no interaction, and the Gb term gives a repulsive interaction for 
Gb < 0. Thus, for reasonable Fermi couplings that are likely to best represent high 
density QCD (as in equation (2.20)), our model Hamiltonian does not predict any 
pairing in this channel 3 . In the LR channel, the model Hamiltonian gives an attrac- 
tive interaction if Ge + Gm > and is independent of the linear combination of 
couplings Ge — 3Gm- If we solve the gap equation for the (uC^u) (LR) condensate 
with G E = G M = G, A = 40 MeV, /i u = 0.4 GeV, and A = 1 GeV, we find a 
gap of 8 keV and a free energy which is about five orders of magnitude smaller than 
that of the LOFF phase. (If we choose Ge = and Gm > 0, the interaction is still 
attractive but the gap is even smaller.) Therefore, even though for 8/j, > 5fi2 we 
expect single-flavor pairing and consequent (uCj^u) and (dCj^d) condensates, the 
resulting condensation energy is so small that it is a good approximation to neglect 

these condensates in the evaluation of 8{J,2, as we have done. 

3 This conclusion is modified at asymptotically high density, where it has been shown by Schafcr 
[63] that long-range single-gluon exchange does in fact induce pairing in the LL/RR channel. (The 
long-range interaction emphasizes near-collinear scattering which is attractive for both electric and 
magnetic gluons.) Simultaneously, the long-range gluon interaction considerably widens the window 
in which the LOFF plane wave phase can occur [46, 47]. Both effects must be taken into account in 
the comparison between LOFF and single-flavor phases in the high density limit 
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2.6 Conclusions 



We have studied the formation of a plane-wave LOFF state involving pairing between 
two flavors of quark whose chemical potentials differ by 28/j,. This state is character- 
ized by a gap parameter and a diquark condensate, but not by an energy gap in the 
dispersion relation. In the LOFF state, each Cooper pair carries momentum 2q with 
|q| « 1.25/i. The condensate and gap parameter vary in space like a plane wave with 
wavelength 7r/|q|. 

We focused primarily on an NJL-type four-fermion interaction with the quantum 
numbers of single gluon exchange. In the limit of weak coupling (BCS gap A <C fi) 
the LOFF state is favored for values of 8/j, which satisfy < 5fi < 5fi 2 , where 
5fi 1 /A = 0.707 and <5/i 2 /A = 0.754. The LOFF gap parameter decreases from 
0.23Ao at 5fi = Sn 1 to zero at 5fi = 5fi 2 - These are the same results found by LOFF 
in their original analysis. Except for very close to 5fi 2 , the critical temperature 
above which the LOFF state melts will be much higher than typical neutron star 
temperatures. At stronger coupling the LOFF gap parameter decreases relative to 
Ao and the window of Sh/Aq within which the LOFF state is favored shrinks. The 
window grows if the interaction is changed to weight electric gluon exchange more 
heavily than magnetic gluon exchange. 

Because it violates rotational invariance by involving Cooper pairs whose mo- 
menta are not antiparallel, the quark matter LOFF state necessarily features nonzero 
condensates in both the J = and J = 1 channels. Both condensates are present 
even if there is no interaction in the J = 1 channel. In this case, however, the J = 1 
condensate does not affect the quasiparticle dispersion relations; that is, the J = 1 
gap parameter vanishes. If there is an attraction in the J = 1 channel (as, for exam- 
ple, if the strength of the electric gluon interaction is increased) the size of the LOFF 
window increases. 

The single plane wave state is a rather feeble state, as we can see from figure 2-4: 
it only prevails as the favored ground state in a tiny window of 5fi, and its gap and 
free energy are small compared to those of the BCS phase. But our results for the 
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single-plane-wave state are just the starting point for an exploration of much more 
complicated and much more robust LOFF crystals involving multiple plane waves. 
We embark on this exploration in the next chapter. 
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Chapter 3 

Crystalline Superconductivity: 
Multiple Plane Waves 

3.1 Overview 

In this chapter we explore LOFF crystalline states that are superpositions of multiple 
plane waves in position space [43] . In momentum space, each plane wave corresponds 
to another "pairing ring" on each Fermi surface, as in figure 1-3. These multiple- 
plane-wave states are much more robust than the single-plane-wave state that we 
studied in chapter 2, for the simple reason that adding more plane waves utilizes 
more of the Fermi surface for pairing, with a corresponding gain in condensation 
energy. 

In studying the multiple-plane-wave states, we do not attempt to write down 
variational wave functions as we did for the single-plane-wave state in chapter 2. 
Rather, in section 3.2, we use a Nambu-Gorkov formalism to directly obtain an infinite 
set of coupled gap equations for the gap parameters {A q }. In principle, finding the 
best LOFF crystal requires exhaustively searching the infinite-dimensional parameter 
space of these crystal order parameters, to find the unique solution of the coupled gap 
equations that is a global minimum of the free energy. This is an infinite task. To 
organize this search, we construct a Ginzburg-Landau free energy function f2({A q }) 
which is valid in the vicinity of the LOFF second-order critical point 5fi 2 - In section 
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3.2 we show how the Nambu-Gorkov method can be used to systematically calculate 
the coefficients in the Ginzburg-Landau potential. Then in section 3.3 we use this 
Ginzburg-Landau potential to study a large catalog of candidate crystal structures. 
It turns out that the potential is unstable and therefore predicts a strong first-order 
phase transition, rather than a second-order transition as obtained for the single- 
plane-wave case, and is therefore not under quantitative control. Nevertheless, the 
Ginzburg-Landau approach is useful and indeed quite powerful because it organizes 
the calculation in such a way that simple qualitative lessons emerge which tell us what 
features make a particular crystal structure energetically favorable or unfavorable. 
This narrows our search, remarkably, to a uniquely favored crystal structure. In 
section 3.4 we discuss this structure, an FCC crystal consisting of eight plane wave 
condensates with wave vectors pointing towards the eight corners of a cube. We 
predict that the phase is quite robust, with gaps comparable in magnitude to the 
BCS gap that would form if the Fermi momenta were degenerate. 

3.2 Methods 

3.2.1 The gap equation 

We study the crystalline superconducting phase in a toy model for QCD that has two 
massless flavors of quarks and a pointlike interaction. The Lagrange function is 

£ = ^ + ^ - |A(VT^)(</TW0 (3.1) 
o 

where /& = 7°(/2 — T^dfi). The r's are Pauli matrices in flavor space, so the up and 
down quarks have chemical potentials as in (1.9). The vertex is r A = 7^T a so that 
our pointlike interaction mimics the spin, color, and flavor structure of one-gluon 
exchange. (The T a are color SU(3) generators normalized so that tr(T a T b ) = 25 ab .) 
We denote the coupling constant in the model by A. 

It is convenient to use a Nambu-Gorkov diagrammatic method to obtain the gap 
equation for the crystalline phase. Since we are investigating a phase with spatial 
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inhomogeneity, we begin in position space. We introduce the two- component spinor 
ty(x) = (ip(x),il! T (x)) and the quark propagator iS(x,x') = (x)^f (x')) , which has 
"normal" and "anomalous" components G and F, respectively: 

iG(x,x') iF(x,x')\ ( mx)j(x')) (^(x)f(x')) 
tS{x,x) = \ _ =1 _ • (3.2) 

iF(x,x') iG(x,x') J \(ij T (x)ij(x')} (^ T {x)^ T (x')) 

The conjugate propagators F and G satisfy 

iG(x,x') = ^°{iG{x\x))^° (3.3) 
iF(x,x') = 7°(2F(x / ,x)) t 7°- (3-4) 

The gap parameter A(x) that describes the diquark condensate is related to the 
anomalous propagator F by a Schwinger-Dyson equation 

A(x) = i-XT A F(x, x)T T A (3.5) 

illustrated diagrammatically in Fig. 3-1. In our toy model, we are neglecting quark 
masses and thus the normal part of the one-particle-irreducible self-energy is zero; 
the anomalous part of the 1PI self energy is just A(x). The crystal order parameter 
A(x) defined by (3.5) is a matrix in spin, flavor and color space. In the mean- 
field approximation, we can use the equations of motion for ty{x) to obtain a set 
of coupled equations that determine the propagator functions in the presence of the 
diquark condensate characterized by A(x): 



i0+ £ A(x) \ I G{x,x') F(x,x') 
A(x) (ifi-fif) \F(x,x') G(x,x') 



H 




: 









5 {4 \x-x') (3.6) 



where A(x) = 7°A(x)"'"7 . Any function A(x) that solves equations (3.5) and (3.6) 
is a stationary point of the free energy functional Q[A(x)]; of these stationary points, 
the one with the lowest Q describes the ground state of the system. Our task, then, 
is to invert (3.6), obtaining F in terms of A(x), substitute in (3.5), find solutions for 
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Figure 3-1: The Schwinger-Dyson graph for the LOFF gap parameter A. The black 
dot is the pointlike interaction vertex and the double line represents the full anomalous 
propagator F, which is given in terms of A in Fig. 3-3. 



A(x), and then evaluate Q for all solutions we find. 

There are some instances where analytic solutions to equations (3.5) and (3.6) 
can be found. The simplest case is that of a spatially uniform condensate. Transla- 
tional invariance then implies that the propagators are diagonal in momentum space: 
S(p,p') = S(p)(2ti) 4: 8^ (p — p'). In this case, Eqs. (3.6) immediately yield 



which is easily inverted to obtain S, which can then be substituted on the right-hand 
side of Eq. (3.5) to obtain a self-consistency equation (i.e. a gap equation) for A. 
The solution of this gap equation describes the familiar "2SC" phase [6, 7, 9, 10, 18], 
a two-flavor, two-color BCS condensate, 



where T 2 , T2, and C75 indicate that the condensate is a color antitriplet, flavor singlet, 
and Lorentz scalar, respectively. 1 The remaining factor A , which without loss of 
generality can be taken to be real, gives the magnitude of the condensate. In order 
to solve the resulting gap equation for Ao, we must complete the specification of 
our toy model by introducing a cutoff. In previous work [9, 10, 11, 18], it has been 
shown that if, for a given cutoff, the coupling A is chosen so that the model describes a 

1 The QCD interaction, and thus the interaction in our toy model, is attractive in the color and 
flavor antisymmetric channel and this dictates the color-flavor pattern of (3.8). Our toy model 
interaction docs not distinguish between the Lorentz scalar (3.8) and the pseudoscalar possibility. 
But, the instanton interaction in QCD favors the scalar condensate. 




(3.7) 



A = T 2 r 2 C 75 A 
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reasonable vacuum chiral condensate, then at /x ~ 400 — 500 MeV the model describes 
a diquark condensate which has A of order tens to 100 MeV. Ratios between physical 
observables depend only weakly on the cutoff, meaning that when A is taken to vary 
with the cutoff such that one observable is held fixed, others depend only weakly 
on the cutoff. For this reason, we are free to make a convenient choice of cutoff so 
long as we then choose the value of A that yields the "correct" A . Since we do not 
really know the correct value of A and since this is after all only a toy model, we 
simply think of Ao as the single free parameter in the model, specifying the strength 
of the interaction and thus the size of the BCS condensate. Because the quarks near 
the Fermi surface contribute most to pairing, it is convenient to introduce a cutoff u 
defined so as to restrict the gap integral to momentum modes near the Fermi surface 
(||p| — p,\ < uj). In the weak coupling (small A) limit, the explicit solution to the gap 
equation is then 

A = 2Lue- n2/2Xp2 . (3.9) 

This is just the familiar BCS result for the gap. (Observe that the density of states 
at the Fermi surface is Nq = 2/x 2 /vr 2 .) We denote the gap for this BCS solution by 
A , reserving the symbol A for the gap parameter in the crystalline phase. We shall 
see explicitly below that when we express our results for A relative to A , they are 
completely independent the cutoff uj as long as A //i is small. 

The BCS phase, with Ao given by (3.9), has a lower free energy than unpaired 
quark matter as long as S/j < 5 fa = Aq/a/2 [90]. The first-order unpairing transition 
at S/j, = 8 fa is the analogue in our two-flavor toy model of the unlocking transition 
in QCD. For 5fi > 5 fa, the free energy of any crystalline solution we find below must 
be compared to that of unpaired quark matter; for 5fi < 5 fa, crystalline solutions 
should be compared to the BCS phase. We shall work at 5fi > 5 fa. 

The simplest example of a LOFF condensate is one that varies like a plane wave: 
A(x) = Aexp(— i2q ■ x). The condensate is static, meaning that q = (0, q). We shall 
denote |q| by go- In this condensate, the momenta of two quarks in a Cooper pair is 
(p + — P + q) f° r some p, meaning that the total momentum of each and every pair 
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is 2q. See Refs. [42, 45] for an analysis of this condensate using the Nambu-Gorkov 
formalism. Here, we sketch the results. If we shift the definition of \P in momentum 
space to ^q(p) = (i>(j> + q), i^{—p + ?)), then in this shifted basis the propagator is 
diagonal: 

iS q (p : p') = (* ? (p)* s (pO> = iS q (p)5\p - p 1 ) (3.10) 
and the inverse propagator is simply 




See Refs. [42, 45] for details and to see how this equation can be inverted and substi- 
tuted into Eq. (3.5) to obtain a gap equation for A. This gap equation has nonzero 
solutions for 5fi < 8/j, 2 — 0.7544A , and has a second-order phase transition at 
8[i — Sfi 2 with A ~ (Sfi 2 — 5jj) 1 / 2 . We rederive these results below. 

If the system is unstable to the formation of a single plane-wave condensate, we 
might expect that a condensate of multiple plane waves is still more favorable. Again 
our goal is to find gap parameters A(x) that are self-consistent solutions of Eqs. (3.5) 
and (3.6). We use an ansatz that retains the Lorentz, flavor, and color structure of 
the 2SC phase: 

A(x) = T 2 r 2 C 75 A(x) (3.12) 

but now A(x) is a scalar function that characterizes the spatial structure of the 
crystal. We write this function as a superposition of plane waves: 

A(x) = J2 A n e ~ t2q ' X ( 3 - 13 ) 
q 

where, as before, q = (0,q). The {A q } constitute a set of order parameters for the 
crystalline phase. Our task is to determine for which set of q's the A q 's are nonzero. 
Physically, for each A q ^ the condensate includes some Cooper pairs for which the 
total momentum of a pair is 2q. This is indicated by the structure of the anomalous 
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Figure 3-2: A process whereby a quark with momentum p scatters by interactions 
with two plane- wave condensates and acquires a momentum p — 2qi + 2q 2 . 

propagator F in momentum space: Eqs. (3.5) and (3.13) together imply that 

F(p,p') = -i(^(p)^(- p >)) = Y,F< i (p)(^) 4 5 {4) ( P -p'-2q) (3.14) 

q 

and 

A q = ^A J ^f A Up)^l (3-15) 

where A q = T 2 T 2 C'y 5 A (l . Eq. (3.15) yields an infinite set of coupled gap equations, 
one for each q. (Note that each F q depends on all the A q 's.) It is not consistent 
to choose only a finite set of A q to be nonzero because when multiple plane-wave 
condensates are present, these condensates induce an infinite "tower" (or lattice) of 
higher momentum condensates. This is easily understood by noting that a quark 
with momentum p can acquire an additional momentum 2q 2 — 2qi by interacting 
with two different plane-wave condensates as it propagates through the medium, as 
shown in Fig. 3-2. Note that this process cannot occur when there is only a single 
plane- wave condensate. The analysis of the single plane- wave condensate closes with 
only a single nonzero A q , and is therefore much easier than the analysis of a generic 
crystal structure. Another way that this difficulty manifests itself is that once we move 
beyond the single plane-wave solution to a more generically nonuniform condensate, 
it is no longer possible to diagonalize the propagator in momentum space by a shift, 
as was possible in Eqs. (3.10, 3.11). 

3.2.2 The Ginzburg-Landau approximation 

The infinite system of equations (3.15) has been solved analytically only in one 
dimension, where it turns out that the gap parameter can be expressed as a Ja- 
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cobi elliptic function that, as promised, is composed of an infinite number of plane 
waves [87]. In three dimensions, the crystal structure of the LOFF state remains 
unresolved [92, 109, 107]. In the vicinity of the second-order transition at <5/i 2 , how- 
ever, we can simplify the calculation considerably by utilizing the smallness of A to 
make a controlled Ginzburg-Landau approximation. This has the advantage of pro- 
viding a controlled truncation of the infinite series of plane waves, because near 5fi 2 
the system is unstable to the formation of plane-wave condensates only for q's that 
fall on a sphere of a certain radius qo, as we shall see below. This was in fact the 
technique employed by Larkin and Ovchinnikov in their original paper [65], and it 
has been further developed in Refs. [92, 109, 107]. As far as we know, though, no 
previous authors have done as complete a study of possible crystal structures in three 
dimensions as we attempt. Most have limited their attention to, at most, structures 
1, 2, 5 and 9 from the 23 structures we describe in Fig. 3-9 and Table 3.1 below. As 
far as we know, no previous authors have investigated the crystal structure that we 
find to be most favorable. 

The authors of Refs. [92, 109, 107] have focused on using the Ginzburg-Landau 
approximation at nonzero temperature, near the critical temperature at which the 
LOFF condensate vanishes. Motivated by our interest in compact stars, we follow 
Larkin and Ovchinnikov in staying at T = while using the fact that, for a single 
plane-wave condensate, A — > for 5fi — > 5fi 2 to motivate the Ginzburg-Landau 
approximation. The down side of this is that, in agreement with previous authors, 
we find that the T = phase transition becomes first order when we generalize 
beyond a single plane wave. In the end, therefore, the lessons of our Ginzburg- 
Landau approximation must be taken qualitatively. We nevertheless learn much that 
is of value. 

To proceed with the Ginzburg-Landau expansion, we first integrate equations (3.6) 
to obtain 

G(x,x') = G {0) (x,x') - / d*z G (0) ( x, z)A(z)F(z, x 1 ) (3.16) 
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Figure 3-3: The diagrammatic expression for the full anomalous propagator F, and 
the first three terms in the series expansion in powers of A. 

F(x,x') = - J d 4 z G (0) (x,z)A(z)G(z,x') (3.17) 

where = (ift + Z^)" 1 , G^ = — ^l) t )~ 1 . Then we expand these equations 
in powers of the gap function A(x). For F(x,x') we find (suppressing the various 
spatial coordinates and integrals for notational simplicity) 

F = -QV*AG<® -PAG^A^AGC) 

(3.18) 

-G^AG^AG^AG^AG^AG^ + 0(A 7 ) 

as expressed diagrammatically in Fig. 3-3. We then substitute this expression for F 
into the right-hand side of the Schwinger- Dyson equation (actually the conjugate of 
equation (3.15)). After some spin, color, and flavor matrix manipulation, the result 
in momentum space is 



A* = --fn( q )A;--A £ j (qiq2 q 3 q)A; i A q2 A* 3 5 q 



J qi-q2+q3-q 

qi,q2,q3 

2A/i 2 



E ^(qiq2q 3 q4q 5 q)A; i A q2 A; 3 A q4 A; 5 5 qi 

-q2+q3-q4+q5-q 
qi,q2,q3,q4,q5 

+C(A 7 ) (3.19) 



as shown in Fig. 3-4. The prefactors have been chosen for later convenience. The 
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Figure 3-4: The series expansion corresponding to Eqs. (3.19) and (3.20). This dia- 
grammatic equation is obtained by substituting the series expansion of Fig. 3-3 into 
the Schwinger-Dyson equation of Fig. 3-1. 



functions II, J, and K corresponding to the three graphs in Fig. 3-4 are given by: 

n(q) = -ij 2 f 7047^-^r 1 (^+2^+^rv 

</(qiq 2 q 3 q4) = -i^ J J^klM ~ fidY\t + Hi + D~ l 

x(^+ Hi - 2^ 2 - ^-\i> + Hi - Hi + Hz + Au)~ lf f 

n 2 f d 4 p 

^(qiq2q3q4q 5 qe) = J j^ni^- fidY 1 (tf + Hi + fin)" 1 

x(^+ 2i x - 24 2 - i/i d y\i> + Hi - H 2 + 2(/ 3 + KY 1 

x^+2^-2^ + 2^-2^-^)- 1 

x(^+ 2^ - 2^ 2 + 2# 3 - 2# 4 + 2^ 5 + ^ U )"V- 

(3.20) 

We shall see that 8/jL and |q| are both of order A which in turn is of order A . This 
means that all these quantities are much less than p in the weak coupling limit. Thus, 
in the weak coupling limit we can choose the cutoff u such that 5fi, |q| <C u <C p,. In 
this limit, J and K are independent of the cutoff u, as we shall see in appendix A 
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where we present their explicit evaluation. In this limit, 



n(q) 
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log 



log 



7T 



2|q| "°Vlq| — 2 \q z 
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a 



2|q| 
|q| 5fi 



|q| — 5ji 



A 2 



2 bg V 4(q 2 - 8f 



(3.21) 



where we have used the explicit solution to the BCS gap equation (3.9) to eliminate 
the cutoff u in favor of the BCS gap Ao, and where the last equation serves to define 
a. Note that a depends on the cutoff uo only through Ao, and depends only on the 
ratios |q|/A and 5/i/A . 

It will prove convenient to use the definition of a to rewrite the Ginzburg-Landau 
equation Eq. (3.19) as 







-q2+q:i-q 



a(|q|)A q + Yl Ji^^KA^KAi- 

qi,q2,q3 

+ ^(qiq 2 q3q4q 5 q)A; i A q2 A; 3 A q4 A; 5 5 qi 

-q2+q3-q4+qs- 



qi,q2,q3,q4,q5 
+0(A 7 ). 



(3.22) 



To learn how to interpret a, consider the single plane-wave condensate in which 
A q 7^ only for a single q. If we divide equation (3.19) by A*, we see that the 
equation II = — 7r 2 /2A/i 2 , which is to say a = 0, defines a curve in the space of 
(|q|, 5fi) where we can find a solution to the gap equation with A q — > 0, with |q| on 
the curve and for any q. This curve is shown in Fig. 3-5. We shall see below that 
when only one A q is nonzero, the "J sum" and U K sum" in (3.19) are both positive. 
This means that wherever a < 0, i.e. below the solid curve in Fig. 3-5, there are 
solutions with A q ^ for these values of |q|, and wherever a > 0, i.e. above the solid 
curve, there are no single plane-wave solutions to the gap equation. The solid curve 
in Fig. 3-5 therefore marks the boundary of the instability towards the formation of a 
single plane-wave condensate. The highest point on this curve is special, as it denotes 
the maximum value of 5fi for which a single plane-wave LOFF condensate can arise. 



^7 




0.905 

Figure 3-5: Along the solid curve, a(|q|,<5/i) = 0. The maximum 8{i reached by this 
curve is 5fi = 0.754A = Sfi 2 , which occurs at |q| = 0.9051A = 1.19975fi 2 . Along 
the upper and lower dashed curves, a = +0.1 and a = —0.1, respectively. 

This second-order critical point occurs at (|q|,<fyt) = (qo,5fi 2 ) with 5[i 2 — 0.7544A 
and qo/8fi 2 ~ 1.1997, where A is the BCS gap of Eq. (3.9). 

As S/i — » b[i 2 from above, only those plane waves lying on a sphere in momentum 
space with |q| = q are becoming unstable to condensation. If we analyze them one by 
one, all these plane waves are equally unstable. That is, in the vicinity of the critical 
point 5fi 2 , the LOFF gap equation admits plane- wave condensates with A q ^ for 
a single q lying somewhere on the sphere |q| = go- For each such plane wave, the 
paired quarks occupy a ring with opening angle ipo — 2 cos _1 (<5/i/|q|) ~ 67.1° on each 
Fermi surface, as shown in Fig. 1-3. 

3.2.3 The free energy 

In order to compare different crystal structures, with (3.22) in hand, we can now 
derive a Ginzburg-Landau free energy functional fi[A(x)] which characterizes the 
system in the vicinity of 5fi 2 , where A — > 0. This is most readily obtained by noting 
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that the gap equations (3.22) must be equivalent to 



<9A q 



(3.23) 



because solutions to the gap equations are stationary points of the free energy. This 
determines the free energy up to an overall multiplicative constant, which can be 
found by comparison with the single plane-wave solution previously known. The 
result is 



= a(q ) A q A q + ^ Yl J (qiq2q 3 q4)A; i A q2 A; 3 A q4 5 qi „ q2+q3 

q, |q|=go qi-q4, |q*|=go 



3 ^(qiq2q3q4q5q6)A qi A q2 A q3 A q4 A q5 A q6 5 qi _ q2+q3 _ q4+q5 _ q6 
qi-qe, |q»|=<?o 

C(A 8 ) (3.24) 



where N = 2/i 2 / 7r 2 and where we have restricted our attention to modes with |q| = g , 
as we now explain. Note that in the vicinity of |q| = go an d 5 fx = 5 fa, 

a ~ ( 5 »; 5 »A . (3.25) 



5 fa 

We see that for a > (that is, for 5fi > 5 fa) A q = is stable whereas for a < (that 
is, 8/jL < 5 fa), the LOFF instability sets in. In the limit a — >• 0~, only those plane 
waves on the sphere |q| = go are unstable. For this reason we only include these plane 
waves in the expression (3.24) for the free energy. Notice that equation (3.24), which 
we have derived starting from the gap equations (3.15), is the same as equation (1.22) 
which was obtained from generic arguments of translational and rotational symmetry. 
The added power of our gap equation derivation is that it enables us to calculate the 
Ginzburg-Landau coefficients a, J, and K from the microscopic theory. 

We shall do most of our analysis in the vicinity of 5fi = 5 fa, where we choose 
|q| = <Zo = 1-19975 fa as just described. However, we shall also want to apply our 
results at 5fi > 5 fa. At these values of 8 fa we shall choose |q| in such a way as to 
minimize a(|q|), because this minimizes the quadratic term in the free energy Q and 



s9 



thus minimizes the free energy in the vicinity of A — > 0, which is where the Ginzburg- 
Landau analysis is reliable. As Fig. 3-5 indicates, for any given 5/x the minimum value 
of a is to be found at |q| = 1.19975/i. Therefore, when we apply Eq. (3.24) away from 
<5/i 2 , we shall set go = 1.19975//, just as at 5ji = 5ji 2 - As a consequence, the opening 
angle of the pairing rings, ipo = 2 cos _1 (5yu/|q|), is unchanged when we move away 
from 5/i = 5[i 2 ■ 

3.3 Results 
3.3.1 Generalities 

All of the modes on the sphere |q| = go become unstable at <5/i = <5/z 2 . The quadratic 
term in the free energy includes no interaction between modes with different q's, 
and so predicts that A q ^ for all modes on the sphere. Each plane-wave mode 
corresponds to a ring of paired quarks on each Fermi surface, so we would obtain 
a cacophony of multiple overlapping rings, favored by the quadratic term because 
this allows more and more of the quarks near their respective Fermi surfaces to pair. 
Moving beyond lowest order, our task is to evaluate the quartic and sextic terms in the 
free energy. These higher order terms characterize the effects of interactions between 
A q 's with differing q's (between the different pairing rings) and thus determine how 
condensation in one mode enhances or deters condensation in other modes. The 
results we shall present rely on our ability to evaluate J and K, defined in Eqs. (3.20). 
We describe the methods we use to evaluate these expressions in appendix A and 
focus here on describing and understanding the results. We shall see, for example, 
that although the quadratic term favors adding more rings, the higher order terms 
strongly disfavor configurations in which A q 's corresponding to rings that intersect 
are nonzero. Evaluating the quartic and sextic terms in the free energy will enable 
us to evaluate the free energy of condensates with various configurations of several 
plane waves and thereby discriminate between candidate crystal structures. 

A given crystal structure can be described by a set of vectors Q = {q a , q^, ■ ■ ■}, 
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Figure 3-6: Rhombic and hexagonal combinations of q's. On the left is a rhombus 
with qi— q2+q3— q4 = 0. On the right is a hexagon with qi~ q2+q3~ q4+q5~ q6 = 0. 
The edges have equal lengths (|qj| = go)- The shapes are in general nonplanar. 

specifying which plane wave modes are present in the condensate, and a set of gap 
parameters {A qo , A q6 , • • •}, indicating the amplitude of condensation in each of the 
modes. Let us define Q as the group of proper and improper rotations that preserve 
the set Q. We make the assumption that Q is also the point group of the crystal 
itself; this implies that A q = A q / if q' is in the orbit of q under the group action. For 
most (but not all) of the structures we investigate, Q has only one orbit and therefore 
all of the A q 's are equal. 

For a given set Q, the quartic term in the free energy (3.24) is a sum over all 
combinations of four q's that form closed "rhombuses", as shown in Fig. 3-6. The four 
q's are chosen from the set Q and they need not be distinct. By a rhombus we mean a 
closed figure composed of four equal length vectors which will in general be nonplanar. 
A rhombus is therefore characterized by two internal angles (ip, x) with the constraint 
< ip + x < n - Each shape corresponds to a value of the J function (as defined in 
equations (3.20)); the rotational invariance of the J function implies that congruent 
shapes give the same value and therefore </(qiq2q3q4) = J^jX)- So, each unique 
rhombic combination of q's in the set Q that characterizes a given crystal structure 
yields a unique contribution to the quartic coefficient in the Ginzburg-Landau free 
energy of that crystal structure. The continuation to next order is straightforward: 
the sextic term in the free energy (3.24) is a sum over all combinations of six q's 
that form closed "hexagons", as shown in Fig. 3-6. Again these shapes are generally 
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nonplanar and each unique hexagonal combination of q's yields a unique value of the 
K function and a unique contribution to the sextic coefficient in the Ginzburg-Landau 
free energy of the crystal. 

As we previewed in chapter 1, when all of the A q 's are equal, we can evaluate 
aggregate quartic and sextic coefficients (3 and 7, respectively, as sums over all rhombic 
and hexagonal combinations of the q's in the set Q: 

/5 = ^J(D), 7 = ^A'(o). (3.26) 
□ O 

Then, for a crystal with P plane waves, the free energy has the simple form 

^ = PaA 2 + -/?A 4 + -7A 6 + 0(A 8 ) (3.27) 
Nq 2 3 

and we can analyze a candidate crystal structure by calculating the coefficients f3 and 
7 and studying the resultant form of the free energy function. 

If (3 and 7 are both positive, a second-order phase transition occurs at a = 0; near 
the critical point the value of 7 is irrelevant and the minimum energy solution is 

A Hnr) • (3 - 28) 

for Sfx < Sfx 2 (i.e. a < 0). 

If f3 is negative and 7 is positive, the phase transition is in fact first order and 
occurs at a new critical point defined by 

3d 2 

a = a, = ^ . (3.29) 

In order to find the S/x* corresponding to a*, we need to solve 

a (|q|, 'fyO = a (1.19975/4,, 5/4.) = a* = , rrt ■ (3.30) 

lb^7 

Since a* is positive, the critical point 5/4* at which the first-order phase transition 
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occurs is larger than 8fjL 2 . If a* is small, then Sfi^ ~ (1 + a*)5/i 2 . Thus, a crystalline 
color superconducting state whose crystal structure yields a negative (3 and positive 
7 persists as a possible ground state even above 5fi 2 , the maximum 8/j, at which the 
plane- wave state is possible. At the first-order critical point (3.29), the free energy 
has degenerate minima at 



If we reduce Sfi below the minimum with A ^ deepens. Once 5/j, is reduced 
to the point at which the single plane wave would just be starting to form with a 
free energy infinitesimally below zero, the free energy of the crystal structure with 
negative j3 and positive 7 has 



at S/j, = 8(i 2 - 

Finally, if 7 is negative, the order A 6 Ginzburg-Landau free energy is unbounded 
from below. In this circumstance, we know that we have found a first order phase 
transition but we do not know at what 8fi^ it occurs, because the stabilization of the 
Ginzburg-Landau free energy at large A must come about at order A 8 or higher. 

3.3.2 One wave 

With these general considerations in mind we now proceed to look at specific examples 
of crystal structures. We begin with the single plane- wave condensate (P = 1). The 
quartic coefficient of the free energy is 



A = 0, 




(3.31) 




(3.32) 



(3 = J = J(0,0) 



1 



1 



0.569 



(3.33) 



4 q2 - 8fi 



- ~ 4 

2 — ~ 



8/j, 



2 ' 
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Figure 3-7: Two different "hexagonal" shapes (as in Fig. 3-6) that can be constructed 
from two vectors q a and q&. These shapes correspond to the functions Ki(ip) and 
K 2 (ip) in Eq. (3.36). 



and the sextic coefficient is 



„ „, 1 q 2 + 35/i 2 1.637 

7 = K = K (qqqqqq) = ^7-^, - (3-34) 



yielding a second-order phase transition at a = 0. These coefficients agree with those 
obtained by expanding the all-orders-in-A solution for the single plane wave which 
can be obtained by variational methods [66, 89, 41] or by starting from (3.11), as 
in Refs. [42, 45]. The coefficient (3 in (3.33) was first found by Larkin and Ovchin- 
nikov [65]. 

3.3.3 Two waves 

Our next example is a condensate of two plane waves (P = 2) with wave vectors q<1 
and q ^ and equal gaps A qa = A q6 = A. The most symmetrical arrangement is an 
antipodal pair (q^ = — q a ), which yields a cosine spatial variation A(x) ~ cos(2q a -x). 
We will find it useful, however, to study the generic case where q a and q^ have the 
same magnitude but define an arbitrary angle tp. We find that the quartic coefficient 
is 

PW) = 2J + 4J(^,0) (3.35) 

and the sextic coefficient is 

j(ip) = 2K + UK^tfj) + QK 2 (tP) (3.36) 



94 



WW 

10 




^ = ~ 67.1° 



(3(tfj ) ~ -1.138/5/i 2 



60°^@e e ^JL20° 150° 180° 



50 r 




20 
10 



30° 



60 c 



7(^0) ^ +0.249/V 



90° 120° 150° 180° 

Figure 3-8: f3(ip) and 7(^)5 the quartic and sextic coefficients in the Ginzburg-Landau 
free energy for a condensate consisting of two plane waves whose wave vectors define 
an angle ip. 

where K^ijj) = -^(qaqaqaqaqbqfc) and K 2 {ip) = if (q a q a q 6 q a q a q b ). (K\ and K 2 arise 
from the "hexagonal" shapes shown in Fig. 3-7.) The functions (3{ip) and are 
plotted in Fig. 3-8. These functions manifest a number of interesting features. Notice 
that the functions are singular and discontinuous at a critical angle ip = ipo ~ 67.1°, 
where is the opening angle of a LOFF pairing ring on the Fermi surface. For 
the two-wave condensate we have two such rings, and the two rings are mutually 
tangent when ip = ipQ. For ip < ipo, both j3 and 7 are large and positive, implying 
that an intersecting ring configuration is energetically unfavorable. For ip > ipo 
the functions are relatively flat and small, indicating some indifference towards any 
particular arrangement of the nonintersecting rings. There is a range of angles for 
which j3 is negative and a first-order transition occurs (note that 7 is always positive). 
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The favored arrangement is a pair of adjacent rings that nearly intersect (ip = ipo + e). 

It is unusual to find coefficients in a Ginzburg-Landau free energy that behave 
discontinuously as a function of parameters describing the state, as seen in Fig. 3-8. 
These discontinuities arise because, as we described in the caption of Fig. 1-3, we are 
taking two limits. We first take a double scaling weak coupling limit (as discussed 
on page 29) in which A , A, 5fi, |q| <C uj <C ft while 5[i/A , |q|/A and A/A 
(and thus the angular width of the pairing bands) are held fixed. Then, we take the 
Ginzburg-Landau limit in which 5/x/Ao — > 5fi 2 /^-o and A/Aq — > and the pairing 
bands shrink to rings of zero angular width. In the Ginzburg-Landau limit, there 
is a sharp distinction between ip < 67.1° where the rings intersect and ip > 67.1° 
where they do not. Without taking the weak coupling limit, the plots of f3(ip) and 
7("0) would nevertheless look like smoothed versions of those in Fig. 3-8, smoothed 
on angular scales of order Sfx 2 /u 2 (for (3) and 5^/uj a (for 7). However, in the weak 
coupling limit these small angular scales are taken to zero. Thus, the double scaling 
limit sharpens what would otherwise be distinctive but continuous features of the 
coefficients in the Ginzburg-Landau free energy into discontinuities. 

3.3.4 Crystals 

From our analysis of the two- wave condensate, we can infer that for a general multiple- 
wave condensate it is unfavorable to allow the pairing rings to intersect on the Fermi 
surface. For nonintersecting rings, the free energy should be relatively insensitive to 
how the rings are arranged on the Fermi surface. However, Eqs. (3.26) suggest that 
a combinatorial advantage is obtained for exceptional structures that permit a large 
number of rhombic and hexagonal combinations of wave vectors. That is, if there are 
many ways of picking four (not necessarily different) wave vectors from the set of wave 
vectors that specify the crystal structure for which qx — q 2 + q3 — q4 = 0, or if there 
are many ways of picking six wave vectors for which qi — q 2 + q3 — q4 + qs — q6 = 0, 
such a crystal structure enjoys a combinatorial advantage that will tend to make the 
magnitudes of (5 or 7 large. For a rhombic combination qi — q 2 + q3 — q4 = 0, the four 
q's must be the four vertices of a rectangle that is inscribed in a circle on the sphere 
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|q| = q . (The circle need not be a great circle, and the rectangle can degenerate 
to a line or a point if the four q's are not distinct). For a hexagonal combination 
qi — q2 + q3 — q4 + q5 — q6 = 0, the triplets (qiqsqs) and (q2q4qe) ar e vertices of two 
inscribed triangles that have a common centroid. In the degenerate case where only 
four of the six q's are distinct, the four distinct q's must be the vertices of an inscribed 
rectangle or an inscribed isosceles trapezoid for which one parallel edge is twice the 
length of the other. When five of the six q's are distinct, they can be arranged as a 
rectangle plus any fifth point, or as five vertices of an inscribed cuboid arranged as 
one antipodal pair plus the three corners adjacent to one of the antipodes. 

We have investigated a large number of different multiple-wave configurations 
depicted in Fig. 3-9 and the results are compiled in Table 3.1. The name of each 
configuration is the name of a polygon or polyhedron that is inscribed in a sphere 
of radius go; the P vertices of the given polygon or polyhedron then correspond to 
the P wave vectors in the set Q. With this choice of nomenclature, keep in mind 
that what we call the "cube" has a different meaning than in much of the previous 
literature. We refer to an eight plane-wave configuration with the eight wave vectors 
directed at the eight corners of a cube. Because this is equivalent to eight vectors 
directed at the eight faces of an octahedron — the cube and the octahedron are dual 
polyhedra — in the nomenclature of previous literature this eight-wave crystal would 
have been called an octahedron, rather than a cube. Similarly, the crystal that we 
call the "octahedron" (six plane waves whose wave vectors point at the six corners of 
an octahedron) is the structure that has been called a cube in the previous literature, 
because its wave vectors point at the faces of a cube. 

Because the LOFF pairing rings have an opening angle ipo ~ 67.1°, no more than 
nine rings can be arranged on the Fermi surface without any intersection [94, 95]. 
For this reason we have focussed on crystal structures with nine or fewer waves, but 
we have included several structures with more waves in order to verify that such 
structures are not favored. We have tried to analyze a fairly exhaustive list of candi- 
date structures. All five Platonic solids are included in Table 3.1, as is the simplest 
Archimedean solid, the cuboctahedron. (All other Archimedean solids have even 
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Figure 3-9: Stereographic projections of the candidate crystal structures. The points 
(•) and circles (O) are projections of q's that are respectively above and below the 
equatorial plane of the sphere |q| = q . 
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Table 3.1: Candidate crystal structures with P plane waves, specified by their sym- 
metry group Q and Foppl configuration. Bars denote dimensionless equivalents: 
j3 = PS/J 2 , 7 = 7<5/i 4 , = Vl/(5nlN ) with N = 2/i 2 /7r 2 . f2 min is the (dimensionless) 
minimum free energy at 8fi = Sn 2 , obtained from (3.32). The phase transition (first 
order for (3 < and 7 > 0, second order for (3 > and 7 > 0) occurs at <5/v 
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c Minimum 7 at 62 = 7r — 63 ~ 70.5° (a cube with one vertex removed). 
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more vertices.) We have analyzed many dihedral polyhedra and polygons: regular 
polygons, bipyramids, prisms, 2 antiprisms, 3 and various capped or augmented poly- 
hedra 4 . For each crystal structure we list the crystal point group Q and the Foppl 
configuration of the polyhedron or polygon. The Foppl configuration is a list of the 
number of vertices on circles formed by intersections of the sphere with consecutive 
planes perpendicular to the principal symmetry axis of the polyhedron or polygon. 
We use a modified notation where a or a indicates that the points on a given circle are 
respectively eclipsed or staggered relative to the circle above. Note that polyhedra 
with several different principal symmetry axes, namely those with T, O, or I sym- 
metry, have several different Foppl descriptions: for example, a cube is (44) along a 
fourfold symmetry axis or (1331) along a threefold symmetry axis. (That is, the cube 
can equally be described as a square prism or a bicapped trigonal antiprism. This 
should make clear that the singly capped trigonal antiprism of Fig. 3-9 and Table 3.1 
is a cube with one vertex removed.) 

We do not claim to have analyzed all possible crystal structures, since that is 
an infinite task. However, there are several classic mathematical problems regarding 
extremal arrangements of points on a sphere and, although we do not know that our 
problem is related to one of these, we have made sure to include solutions to these 
problems. For example, many of the structures that we have evaluated correspond 
to solutions of Thomson's problem [110, 94] (lowest energy arrangement of P point 
charges on the surface of a sphere) or Tammes's problem [95, 94] (best packing of P 
equal circles on the surface of a sphere without any overlap). In fact, we include all 
solutions to the Thomson and Tammes problems for P < 9. Our list also includes 
all "balanced" configurations [111] that are possible for nonintersecting rings: a bal- 
anced configuration is a set Q with a rotational symmetry about every q e Q; this 

2 A trigonal prism is two triangles, one above the other. A cube is an example of a square prism. 

3 An antiprism is a prism with a twist. For example, a square antiprism is two squares, one 
above the other, rotated relative to each other by 45°. The octahedron is an example of a trigonal 
antiprism. 

4 Capping a polyhedron adds a single vertex on the principal symmetry axis of the polyhedron (or 
polygon). Thus, a capped square is a square pyramid and an octahedron could be called a bicapped 
square. Augmenting a polyhedron means adding vertices on the equatorial plane, centered outside 
each vertical facet. Thus, augmenting a trigonal prism adds three new vertices. 
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corresponds to an arrangement of particles on a sphere for which the particles are in 
equilibrium for any two-particle force law. 

For each crystal structure, we have calculated the quartic and sextic coefficients (3 
and 7 according to Eqs. (3.26), using methods described in appendix A to calculate 
all the J and K integrals. To further discriminate among the various candidate 
structures, we also list the minimum free energy f2 min evaluated at the plane-wave 
instability point 5fi = 5fi 2 where a = 0. To set the scale, note that the BCS state 
at 5/1 — has f^cs = — /-^Ao/V 2 , corresponding to f^Bcs = —0.879 in the units of 
Table 3.1. For those configurations with (3 > and 7 > 0, f2 min = at 5fi = 5fi 2 

and ^min 

< for 5fi < 5/i 2 , where a < 0. Thus, we find a second-order phase 
transition at 5fi = 5fi 2 . F° r those configurations with (3 < and 7 > 0, at 5fi = 5fi 2 
the minimum free energy occurs at a nonzero A with f2 m i n < 0. (The value of A 
at which this minimum occurs can be obtained from (3.32).) Because Q m i n < at 
Sfi = Sfi 2 , if we go to 5fi > 5fi 2 , where a > 0, we lift this minimum until at some 
Sfj,^ it has Q = and becomes degenerate with the A = minimum. At <5/i = S/i^, 
a first-order phase transition occurs. For a very weak first-order phase transition, 
Six* ~ 5ji 2 ~ 0.754A . For a strong first-order phase transition, 5^ ^> 5fi 2 and the 
crystalline color superconducting phase prevails as the favored ground state over a 
wider range of 5 fx. 

3.3.5 Crystal structures with intersecting rings 

There are seven configurations in Table 3.1 with very large positive values for 7. 
These are precisely the seven configurations that have intersecting pairing rings: the 
hexagon, hexagonal prism, capped square prism, bicapped square antiprism, icosa- 
hedron, cub octahedron, and dodecahedron. The first two of these include hexagons, 
and since ipo > 60° the rings intersect. The last four of these crystal structures have 
more than nine rings, meaning that intersections between rings are also inevitable. 
The capped square prism is an example of a nine-wave structure with intersecting 
rings. It has a 7 which is almost two orders of magnitude larger than that of the 
augmented trigonal prism and the capped square antiprism which, in contrast, are 
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nine wave structures with no intersecting rings. Because of their very large 7's all 
the structures with intersecting rings have either second-order phase transitions or 
very weak first-order phase transitions occurring at a 5fi^ ~ <5/i 2 . At 5fi = 5fi 2 , all 
these crystal structures have f2 m i n very close to zero. Thus, as our analysis of two 
plane waves led us to expect, we conclude that these crowded configurations with 
intersecting rings are disfavored. 

3.3.6 "Regular" crystal structures 

At the opposite extreme, we see that there are several structures that have negative 
values of 7: the square, square pyramid, square antiprism, trigonal prism, capped 
trigonal antiprism, and cube. Our analysis demonstrates that the transition to all 
these crystal structures (as to those with (3 < and 7 > 0) is first order. But, we 
cannot evaluate f2 min or 5[i^ because, to the order we are working, Q is unbounded 
from below. For each of these crystal structures, we could formulate a well-posed (but 
difficult) variational problem in which we make a variational ansatz corresponding to 
the structure, vary, and find fi min without making a Ginzburg-Landau approximation. 
It is likely, therefore, that within the Ginzburg-Landau approximation Q will be 
stabilized at a higher order than the sextic order to which we have worked. 

Of the sixteen crystal structures with no intersecting rings, there are seven that 
are particularly favored by the combinatorics of Eqs. (3.26). It turns out that these 
seven crystal structures are precisely the six that we have found with 7 < 0, plus 
the octahedron, which is the most favored crystal structure among those with 7 > 0. 
As discussed earlier, more terms contribute to the rhombic and hexagonal sums in 
Eqs. (3.26) when the q's are arranged in such away that their vertices form rectangles, 
trapezoids, and cuboids inscribed in the sphere |q| = go- Thus the square itself 
fares well, as do the square pyramid, square antiprism, and trigonal prism which 
contain one, two, and three rectangular faces, respectively. The octahedron has three 
square cross sections. However, the cube is the outstanding winner because it has 
six rectangular faces, six rectangular cross sections, and also allows the five-corner 
arrangements described previously. The capped trigonal antiprism in Table 3.1 is a 
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cube with one vertex removed. This seven-wave crystal has almost as many waves 
as the cube, and almost as many combinatorial advantages as the cube, and it turns 
out to have the second most negative 7. 

With eight rings, the cube is close to the maximum packing for nonintersecting 
rings with opening angle ipo ~ 67.1°. Although nine rings of this size can be packed on 
the sphere, adding a ninth ring to the cube and deforming the eight rings into a cuboid, 
as we have done with the capped square prism, necessarily results in intersecting rings 
and the ensuing cost overwhelms the benefits of the cuboidal structure. To form a 
nine-ring structure with no intersections requires rearranging the eight rings, spoiling 
the favorable regularities of the cuboid. Therefore a nine-ring arrangement is actually 
less favorable than the cuboid, even though it allows one more plane wave. We see 
from Table 3.1 that 7 for the cube is much more negative than that for any of the 
other combinatorially favored structures. The cube is our winner, and we understand 
why. 

To explore the extent to which the cube is favored, we can compare it to the 
octahedron, which is the crystal structure with 7 > for which we found the strongest 
first-order phase transition, with the largest 5^ and the deepest Q m i n . The order A 6 
free energy we have calculated for the cube is far below that for the octahedron at 
all values of A. To take an extreme example, at <5/i = 5/i 2 the octahedron has fl min = 
—13.365 at A = 1.263<5/i 2 = 0.953Ao whereas for the cube we find that Q = —2151.5 
at this A. As another example, suppose that we arbitrarily add +^800A 8 /5/i 8 to the 
Q of the cube. In this case, at 5fi = 5fi 2 we find that the cube has n min = —32.5 at 
A = 0.656A and is thus still favored over the octahedron, even though we have not 
added any A 8 term to the free energy of the octahedron. These numerical exercises 
demonstrate the extraordinary robustness of the cube, but should not be taken as 
more than qualitative. We do not know at what A and at what value f2 m i n the true 
free energy for the cube finds its minimum. However, because the qualitative features 
of the cube are so favorable we expect that it will have a deeper f2 m i n and a larger 5fi^ 
than the octahedron. Within the Ginzburg-Landau approximation, the octahedron 
already has A = 0.953A and a deep f2 min = —13.365, about fifteen times deeper 
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than Q BCS = -0.879 for the BCS state at 5fi = 5 . 

Even if we were to push the Ginzburg-Landau analysis of the cube to higher order 
and find a stable f2 m i n , we would not be able to trust such a result quantitatively. 
Because it predicts a strong first-order phase transition, the Ginzburg-Landau ap- 
proximation predicts its own quantitative demise. What we have learned from it, 
however, is that there are qualitative reasons that make the cube the most favored 
crystal structure of them all. And, to the extent that we can trust the quantitative 
calculations qualitatively, they indicate that the first-order phase transition results 
in a state with A and fi m i n comparable to those of the BCS phase, and occurs at a 
5fi^ > 5fi 2 . 

3.3.7 Varying continuous degrees of freedom 

None of the regularities of the cube which make it so favorable are lost if it is deformed 
continuously into a cuboid, slightly shorter or taller than it is wide, as long as it is 
not deformed so much as to cause rings to cross. Next, we investigate this and some 
of the other possible continuous degrees of freedom present in a number of the crystal 
structures we have described above. 

So far we have neglected the fact, mentioned at the start of this section, that some 
of the candidate structures have multiple orbits under the action of the point group 
Q . These structures include the square pyramid, the four bipyramids, and the five 
capped or augmented structures listed in Table 3.1; all have two orbits except for the 
three singly capped crystal structures, which have three orbits. For these multiple- 
orbit structures each orbit should have a different gap parameter but in Table 3.1 
we have assumed that all the gaps are equal. We have, however, analyzed each of 
these structures upon assuming different gaps, searching for a minimum of the free 
energy in the two- or three-dimensional parameter space of gaps. In most cases, 

the deepest minimum is actually obtained by simply eliminating one of the orbits 

5 Of course, the BCS state isotropically utilizes the entire Fermi surface for pairing, so we will not 
find a crystalline state which has a condensation energy greater than that of the cube. The numbers 
arc unreliable here, but can be taken to indicate that the condensation energy could be comparable 
to that of the BCS phase. 
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from the configuration (i.e. let A = for that orbit); the resultant structure with 
one less orbit appears as another structure in Table 3.1. For example, the bicapped 
square antiprism has two orbits: the first is the set of eight q's forming a square 
antiprism, the second is the antipodal pair of q's forming the two "caps" of the 
structure. Denote the gaps corresponding to these two orbits as Ai and A 2 . This 
structure is overcrowded with intersecting rings, so it is not surprising to find that 
a lower-energy configuration is obtained by simply letting A 2 = 0, which gives the 
"uncapped" square antiprism. Configurations with fewer orbits are generally more 
favorable, with only three exceptions known to us: the trigonal bipyramid is favored 
over the triangle or the antipodal pair; the square pyramid is favored over the square 
or the point; and the capped trigonal antiprism is favored over any of the structures 
that can be obtained from it by removing one or two orbits. For these configurations, 
Table 3.1 lists the results for Ai = A 2 (= A3); the numbers can be slightly improved 
with Ai 7^ A 2 (7^ A3) but the difference is unimportant. 

For some configurations in Table 3.1 the positions of the points are completely fixed 
by symmetry while for others the positions of the points can be varied continuously, 
while still maintaining the point-group symmetry of the structure. For example, 
with the square pyramid we can vary the latitude of the plane that contains the 
inscribed pyramid base. Similarly, with the various polygonal prism and antiprism 
structures (and associated cappings and augmentations), we can vary the latitudes 
of the inscribed polygons (equivalently, we can vary the heights of these structures 
along the principal symmetry axis). For each structure that has such degrees of 
freedom, we have scanned the allowed continuous parameter space to find the favored 
configuration. Table 3.1 then shows the results for this favored configuration, and the 
latitude angles describing the favored configuration are given as footnotes. However, 
if the structure always has overlapping rings regardless of its deformation, then either 
no favorite configuration exists or the favorite configuration is a degenerate one that 
removes the overlaps by changing the structure. There are two instances where this 
occurs: the capped square prism can be deformed into a square pyramid by shrinking 
the height of the square prism to zero, and the bicapped square antiprism can be 



105 




Figure 3-10: The sextic free energy coefficient 7 for the square antiprism as a function 
of the polar (latitude) angle 9 of the top square facet. (The polar angle of the bottom 
square facet is it — 9.) The inset plot shows the detail in the range of 9 where no 
rings intersect. Solid and dashed vertical lines indicate the positions of primary and 
secondary singularities as discussed in the text (other secondary singularities occur, 
but are not discernible on the plot). 

deformed into an antipodal pair by moving the top and bottom square faces of the 
antiprism to the north and south poles, respectively. For these structures, Table 3.1 
just lists results for an arbitrarily chosen nondegenerate configuration. 

A typical parameter scan is shown in Fig, 3-10, where we have plotted 7 for 
the square antiprism as a function of the polar angle 9 of the top square facet (the 
polar angle of the bottom square facet is ir — 9). As we expect, 7 is very large in 
regions where any rings intersect, and we search for a minimum of 7 in the region 
where no rings intersect. The plot has a rather complicated structure of singularities 
and discontinuities; these features are analogous to those of Fig. 3-8, and as there 
they ciriSG clS cL result of the double limit we are taking. Primary singularities occur 
at critical angles where pairing rings are mutually tangent on the Fermi surface. 
Secondary singularities occur where rings corresponding to harmonic q's, obtained by 
taking sums and differences of the fundamental q's that define the crystal structure, 
are mutually tangent. Such q's arise in the calculation of J and K because these 
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Figure 3-11: The sextic free energy coefficient 7 for a scan that linearly interpolates 
from the cube (twist angle <p = 0) to the square antiprism of Table 3.1 (twist angle 
(f) = 45°). Dashed vertical lines indicate secondary singularities. 

calculations involve momenta corresponding to various diagonals of the rhombus and 
hexagon in Fig. 3-6. 

In addition to varying the latitudes of the Foppl planes in various structures, we 
varied "twist angles". For example, we explored the continuous degree of freedom 
that turns a cube into a square antiprism, by twisting the top square relative to 
the bottom square by an angle <fi ranging from 0° to 45°. In Fig. 3-11, we show a 
parameter scan in which we simultaneously vary the twist angle <fi and the latitudes 
of the square planes in such a way that the scan interpolates linearly from the cube 
to the most favorable square antiprism of Table 3.1. In this parameter scan, we find 
a collection of secondary singularities and one striking fact: 7 is much more negative 
when the twist angle is zero {i.e. for the cube itself) than for any nonzero value. 
For the cube, 7 = — 459.2/5/i 4 , whereas the best one can do with a nonzero twist 
is 7 = — 64.2/5/i 4 , which is the result for an infinitesimal twist angle. Thus, any 
nonzero twist spoils the regularities of the cube that contribute to its combinatorial 
advantage, and this has a dramatic and unfavorable effect on the free energy. 

Finally, we have scanned the parameter space of a generic cuboid to see how this 
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compares to the special case of a cube. That is, we vary the height of the cuboid 
relative to its width, without introducing any twist. This continuous variation does 
not reduce the combinatorial advantage of the crystal structure. As expected, there- 
fore, we find that as long as the cuboid has no intersecting rings, it has a free energy 
that is very similar to that of the cube itself. Any cuboid with intersecting rings 
is very unfavorable. In the restricted parameter space of nonintersecting cuboidal 
arrangements, the cuboid with the most negative free energy is a square prism with 
a polar angle of 51.4° for the top square face. (For this polar angle the pairing rings 
corresponding to the four corners of the square are almost mutually tangent.) This 
prism is slightly taller than a perfect cube, which has a polar angle of 54.7°. The free 
energy coefficients of the best cuboid are (3 = — 111.563/(5/i) 2 , 7 = — 463.878/(<5/i) 4 . 
These coefficients differ by less than 1% from those for the cube, given in Table 3.1. 
There is no significant difference between the cube and this very slightly more favor- 
able cuboid: all the qualitative arguments that favor the cube favor any cuboid with 
no intersecting rings equally well. We therefore expect that if we could determine 
the exact (rather than Ginzburg-Landau) free energy, we would find that the favored 
crystal structure is a cuboid with a polar angle somewhere between 51.4° and 56.5°, 
as this is the range for which no rings intersect. We expect no important distinc- 
tion between the free energy of whichever cuboid in this narrow range happens to be 
favored and that of the cube itself. 

3.4 Conclusions 

We have argued that the cube crystal structure is the favored ground state at zero 
temperature near the plane-wave instability point 5fi = 5n 2 - By the cube we mean 
a crystal structure constructed as the sum of eight plane waves with wave vectors 
pointing towards the corners of a cube. The qualitative points (which we have demon- 
strated in explicit detail via the analysis of many different crystal structures) that 
lead us to conclude that the cube is the winner are: 

• The quadratic term in the Ginzburg-Landau free energy wants a |q| such that 
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the pairing associated with any single choice of q occurs on a ring with opening 
angle ipo ~ 67.1° on each Fermi surface. 

• The quadratic term in the Ginzburg-Landau free energy favors condensation 
with many different wave vectors, and thus many different pairing rings on 
the Fermi surfaces. However, the quartic and sextic terms in the free energy 
strenuously prohibit the intersection of pairing rings. No more than nine rings 
with opening angle 67.1° can be placed on the sphere without overlap. 

• The quartic and sextic terms favor regular crystal structures, for example those 
that include many different sets of wave vectors whose tips form rectangles. 
None of the nine-wave structures with no intersections between pairing rings 
are regular in the required sense. The cube is a very regular eight-wave crystal 
structure. 

Quantitatively, we find that a cube (actually, a cuboid that is only slightly taller 
than it is wide) has by far the most negative Ginzburg-Landau free energy, to sextic 
order, of all the many crystal structures we have investigated. For the cube structure, 
the eight q vectors are the eight shortest vectors in the reciprocal lattice of a face- 
centered-cubic crystal. Therefore, we find that A(x) exhibits face-centered-cubic 
(FCC) symmetry. The explicit form of the gap function A(x) is shown in equation 
(1.25), and a unit cell of the FCC crystal is shown in figure 1-6. 

Our Ginzburg-Landau analysis predicts a first-order phase transition to the cu- 
bic crystalline color superconductor at some 5fi = 5fi^. The fact that we predict 
a first-order phase transition means that the Ginzburg-Landau analysis cannot be 
trusted quantitatively. Furthermore, at order A 6 , which is as far as we have gone, 
the Ginzburg-Landau free energy for the cube is unbounded from below. We there- 
fore have no quantitative prediction of 8{i* or the magnitude of A. The best we can 
do is to note that the cube is significantly favored over the octahedron, for which 
the order A 6 Ginzburg-Landau analysis predicts ~ 3.6Ao and predicts that at 
8/jL = 8fi 2 = 0.754Ao, the gap is A ~ 0.95Ao and the condensation energy is larger 
than that in the BCS state by a factor of about fifteen. As we have warned repeatedly, 
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these numbers should not be trusted quantitatively: because the Ginzburg-Landau 
approximation predicts a strong first-order phase transition, it predicts its own break- 
down. We have learned several qualitative lessons from it, however: 

• We have understood the qualitative reasons that make the cube the most favored 
crystal structure of them all. It has the largest number of pairing rings that 
can be "regularly" arranged on the Fermi surface without ring overlaps. 

• The cube structure has an unstable Ginzburg-Landau free energy. As a result 
the Ginzburg-Landau analysis cannot give quantitative results, but the large 
instability suggests a robust crystalline phase. The gap parameter A for the 
crystalline phase could be comparable to the gap A for the BCS phase. Sim- 
ilarly, the condensation energy by which the crystalline phase is favored over 
unpaired quark matter at 5fi ^ could be comparable to that for the BCS 
phase at 5/i = 0. 

• Although we began our analysis in the vicinity of the second-order plane-wave 
instability point <5/i 2 , the transition to the crystalline phase is strongly first- 
order and occurs at a 5fi^ ^> 5fi 2 . The emergence of a first-order transition 
from a study of a second-order point has a precedent in the Ginzburg-Landau 
analysis for the liquid-solid transition [93]. 

• We learn that the crystalline color superconductivity window S/j, 1 < 5fi < 5fi^ 
is large. Because 6fi 2 is n °t much larger than Sn 1 , the original LOFF window 
S/j, 1 < 5fi < 5fi 2 wherein the single plane-wave condensate is possible is quite 
narrow. We have learned, however, that <5/x_ ^> 5fi 2 . Furthermore, because the 
condensation energy of the crystalline phase is so robust, much greater than 
that for the single plane wave and likely comparable to that for the BCS phase, 
the value of 5fi 1 , the location of the transition between the crystalline phase 
and the BCS phase, will be significantly depressed. So the LOFF window is 
considerably widened in both directions. 

• Given the robustness of the FCC crystalline phase in our two flavor model, in 
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real QCD with three flavors of quarks we can expect that the crystalline phase 
will occupy a large regions of the (T,/j,) phase diagram of figure 1-1. 
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Chapter 4 

Single Color and Single Flavor 
Color Superconductivity 

4.1 Overview 

There are nine types of quark (3 colors, 3 flavors) in dense quark matter. In general 
we expect that all nine quarks will find attractive channels in which to form Cooper 
pairs, because any pairing lowers the free energy by BCS condensation. The color- 
flavor-locked (CFL) phase directly accomplishes this feat, but for non-CFL phases of 
color superconductivity that occur at intermediate densities as in figure 1-1, it is not 
obvious how the system contrives to pair all nine quarks. Here we speculate on some 
of the contrivances that might occur. 

In the 2SC phase (which seems unlikely to occur in neutral quark matter), only 
two colors and two flavors pair, leaving five "orphaned" quarks: the blue up and down 
quarks, and all three colors of strange quark. In this context it has been proposed that 
the blue up and down quarks could pair together in a single-color condensate [9, 98], 
and all the strange quarks could pair together in a single-flavor condensate which 
involves all three colors [63, 8]. 

Single-flavor pairing can also occur in the limit of very large separation between 
the Fermi surfaces for the different flavors. In this context the system will abandon 
inter-species pairing and form (uu), (dd), and (ss) condensates. Each could involve 
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three colors (e.g. in a color-spin- locked (CSL) phase [63]) and therefore all nine quarks 
are paired. 

In the LOFF crystalline phase, so far we have only discussed the two-flavor sit- 
uation of pairing between up and down quarks. This is a "2SC-LOFF" phase: just 
like the spatially-uniform 2SC phase, it pairs only two colors and two flavors and 
leaves five orphaned quarks. In a real three-flavor context, the crystalline phase could 
involve ud, us, and ds pairs, thus involving all nine quarks in a "CFL-LOFF" phase. 
Alternatively, the three-flavor LOFF state might involve only ud and us pairs, i.e. it 
might only pair quarks with adjacent Fermi surfaces (recall that in neutral quark 
matter the Fermi momenta are split with p s F < p F < p d F ). In this case, the ud pairing 
could involve r and g quarks, and the us pairing could involve r and b quarks, and 
two quarks are orphaned: the blue down quark and the green strange quark. These 
lonely quarks might resort to self-pairing if the QCD interaction permits an attraction 
in this channel (we will see that it does in some NJL models, if the quark is heavy 
enough) . 

In this chapter we survey a large catalog of different color-flavor-spin channels for 
diquark condensation [62]. This catalog includes the single-color and/or single-flavor 
pairing channels that are mentioned in the various scenarios described above [8, 9, 
63, 98]. We consider only translationally invariant phases, but many of the conden- 
sates in our survey have non-zero spin, and therefore spontaneously break rotational 
invariance. In section 4.2 we identify the attractive channels using NJL models with 
four-fermion interactions based on instantons, magnetic gluons, and combined electric 
and magnetic gluons. In section 4.3 we solve the NJL mean-field gap equations for the 
channels that are attractive. In section 4.4 we investigate the quasiquark dispersion 
relations for the spin-one condensates. The quasiquark energy gaps are anisotropic 
and gapless modes can occur at the poles or at the equator of the Fermi surface. 
Finally in section 4.5, we summarize the results, suggest directions for further study, 
and speculate about implications for the physics of compact stars. 
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4.2 Mean-field survey of quark pairing channels 



4.2.1 Calculation 

To see which channels are attractive we perform a mean-field calculation of the pairing 
energy for a wide range of condensation patterns. We write the NJL Hamiltonian in 
the form 

H Hf Iee -\- H interaction 

#frcc = ^-filo + m)^ (4.1) 

n interaction — W a Vj W 7 Wl n ia (3b kc Sd > 

where color indices are a, /3, 7, 5, flavor indices are i, j, k, I, spinor indices are a, b, c, d. 
The four-fermion interaction is supposed to be a plausible model of QCD, so in the 
interaction kernel 7i we include three terms, with the color-flavor-spinor structure of 
a two-flavor instanton, electric gluon exchange, and magnetic gluon exchange, 

inst 

Klec = \G E SjSi 5 ab 5 cd 1^5} - Sffi) 

Hmst = -|Gj e ik e jl J ([70(1 + 75)U[7o(l + 7is)]«i + [7o(l - 75)] a6 [7o(l - 7s)]«i) 
1(35^-^) 

(4.2) 

We consider condensates that factorize into separate color, flavor, and Dirac ten- 
sors (i.e. that do not show "locking") and calculate their binding energy by contracting 
them with (4.2). 

There is no Fierzing ambiguity in this procedure. For a given pairing pattern X, 
the condensate is 

(^V)ipi = Ww . (4.3) 
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We can then calculate the interaction ( "binding" ) energy of the various condensates, 

H = -Y, A(X) 2 (d£>G B + Sg>G M + S£2G 7 ) (4.4) 
x 

(X) 

The binding strengths interaction &i ye the strength of the self-interaction of the con- 
densate X due to the specified part of the interaction Hamiltonian. 

4.2.2 Properties of the pairing channels 

In Table 4.1 we list the the simple (translationally invariant, factorizable) channels 
available for quark pairing. The meanings of the columns are as follows. 

1. Color : two quarks either make an antisymmetric color triplet (which requires 
quarks of two different colors) or a symmetric sextet (which can occur with 
quarks of two different colors, and also if both quarks have the same color). 
For the 3^ we use <£f s = e 135 in Eq. (4.3). For the 65- we use a single-color 
representative <£P S = S^^S 6,1 in Eq. (4.3). 

2. Flavor : two quarks either make an antisymmetric flavor singlet (which requires 
quarks of two different flavors) or a symmetric triplet (which can occur with 
quarks of two different flavors, and also if both quarks have the same flavor). 
For the 1a we use ^7 = and for the 3s we use $ji = crj ; in (4.3). 

3. Spin, parity: since the chemical potential explicitly breaks the Lorentz group 
down to three-dimensional rotations and translations, it makes sense to classify 
condensates by their total angular momentum quantum number j and parity. 

4. Dirac : This column gives the Dirac matrix structure T bd used in (4.3), so the 
condensate is ip T Tip. We also designate each condensate as "LL" (even number 
of gamma matrices, so pairs same-chirality quarks) or "LR" (odd number of 
gamma matrices, so pairs opposite-chirality quarks). 

5. BCS-enhancement : Condensates that correspond to pairs of particles or holes 
near the Fermi surface have a BCS singularity in their gap equation that guar- 
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antees a solution, no matter how weak the coupling. To see which condensates 
have such a BCS enhancement, we expanded the field operators in terms of 
creation and annihilation operators (see Appendix B.3). The order of the co- 
efficient of the a(p)a(— p) and 6^(p)6^(— p) terms is given in the table. 0(1) 
means BCS-enhanced, means not BCS-enhanced. In the C7075 condensate 
the coefficient goes to zero as the quark mass goes to zero (hence it is labelled 
u O(m)" in the table) meaning that the channel loses its BCS enhancement in 
the chiral limit. This is discussed further in Appendix B.3. 

6. Binding strength: For each channel we show the binding strength for the instan- 
ton interaction, the full (electric plus magnetic) gluon, which could reasonably 
be used at medium density, and for the magnetic gluon alone, which is known 
to dominate at ultra-high density [6, 12]. Channels with a positive binding 
strength and BCS enhancement will always support pairing (the gap equation 
always has a solution, however weak the couplings Gi,Ge,Gm)- Other things 
being equal, the pairing with the largest binding strength will have the lowest 
free energy, and is the one that will actually occur. 

It may seem strange that there are entries in the table with angular momentum 
j = 1 and an antisymmetric Dirac structure (C7375), and with j = but a symmetric 
Dirac structure (C70). If all the angular momentum came from spin this would be 
impossible. But even though there are no explicit spatial derivatives in the diquark 
operators, there can still be orbital angular momentum. In Appendix B.3 the angular 
momentum content of the particle-particle component of the condensates is analyzed 
into its spin and orbital content. We see, for example, that C7375 has an antisymmet- 
ric space wavefunction (/ = 1) and a symmetric spin wavefunction (s = 1), combined 
to give an antisymmetric j '• = 1. 

4.2.3 Results 

The results of the binding strength calculation are shown in Table 4.1. The first block 
is antisymmetric in flavor and color, and so describes pairing of two flavors and two 
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Table 4.1: Binding strengths of di quark channels in NJL models in the mean-field 
approximation. The first 6 columns specify the channels, and the last 3 columns give 
their attractiveness in NJL models with various types of four-fermion vertex: 2-flavor 
instanton, single gluon exchange, single magnetic gluon exchange (expected to dom- 
inate at higher density). See equations (4.3) and (4.4) and subsequent explanation. 

colors. The second block is for two flavors and one color, the third for one flavor and 
two colors, and the final block is for one color and one flavor. 

Certain features can be easily understood: the flavor-symmetric condensates all 
have zero instanton binding energy, because the instanton vertex is flavor-antisymmetric 
in the incoming quarks. The gluonic vertices give the same results for C75 as for C, 
and for C 00375 as for Co"o3> because the gluonic interaction is invariant under U(1)a 
transformations, under which the LL condensates transform into each other (C75 C 
and Ccr 03 75 ^ C 003) while the LR condensates are invariant. We see that there are 
many attractive channels: 

1) Two colors and two flavors (3a, 1a,- • •)■ 
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The strongly attractive channel (3^,1^,0,+) (C75) is the 2SC and CFL quark 
Cooper pairing pattern, and has been extensively studied. The gap is large 
enough that even species with different masses, whose Fermi momenta are quite 
far apart, can pair (hence the CFL phase which pairs red and green u and d, 
red and blue u and s, and green and blue d and s in this channel). Its parity 
partner (3^,1^,0,— )(C) is disfavored by instantons, and is therefore unlikely 
to occur at phenomenologically interesting densities. The additional channel 
(3a,^a,^,— XC7375) is more weakly attractive and also breaks rotational invari- 
ance, and is therefore expected to be even less favored. This is confirmed by gap 
equation calculations (Fig. 4-2) which show that its gap is smaller by a factor 
of 10 to 100. 

2) One color, two flavors (65,1,4,- • •)• 

It is generally emphasized that the quark-quark interaction is attractive in the 
color-antisymmetric 3a channel. But, as we see in table 4.1, the color-symmetric 
(65,1^,1,+) (CV03) is attractive for instantons and the magnetic gluon four- 
fermion interaction. The instanton gives it a gap of order 1 MeV (Fig. 4-1), 
while the gluon interaction gives a small gap of order 1 eV (Fig. 4-2). This 
channel was originally suggested for pairing of the blue up and down quarks 
that are left out of 2SC [9], and is discussed in more detail in Ref. [98]. Its gap 
is small, so it could only pair quarks of similar mass, i.e. the light quarks, but 
in a real-world uniform phase such pairing will not occur either, because charge 
neutrality causes the up and down chemical potentials (and hence Fermi mo- 
menta) to differ by tens of MeV, which is larger than the gap. In a non-uniform 
mixture of two locally charged phases [112], however, it is conceivable that the 
up and down Fermi momenta could be similar enough to allow pairing in this 
channel. The parity partner (65,1^,1,— ){Ca^ r )^) is disfavored by instantons. 
The channel (65, 1^,0,— )(C*7o) is attractive, but has no particle-particle compo- 
nent, and presumably only occurs for sufficiently strong coupling. Solving the 
gap equations for reasonable coupling strength we find no gap in this channel. 
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3) Two colors and one flavor (3a,3s,. ■ •)■ 

The only attractive channel is (3a,3s,1,+) (C73). This is a pairing option for 
red and green strange quarks when the up and down quarks are paired in the 
2SC state. . We have solved the relevant gap equation (Figs. 4-2,4-3) and find 
gaps of about 1 MeV or less. If three colors are available then a competing 
possibility is to lock the colors to the spin (CSL), so the condensate is a linear 
combination of C'-fi and Ca 0i with a color structure that is correlated with the 
spatial direction, e.g. red and green quarks pair in the z direction, red and blue 
in the y direction, green and blue in the x direction. This leaves an unbroken 
global SO (3) of spatial rotations combined with color rotations, so the gap 
is isotropic, which helps to lower the free energy [63, 64]. Note also that the 
channels (3^,35, l,+)(C<7o3) and (3^,3^,1,— )(C<r 03 'y 5 ) which are repulsive in the 
NJL model become attractive at asymptotic density when the gluon propagator 
provides a form factor that strongly emphasizes small-angle scattering [63]. 

4) One color and one flavor (65,35,. . .). 

There is an attractive channel here, the (65,35,0,+) (C7075). It loses its particle- 
particle component as the quark mass goes to zero, making it very weak for up 
and down quarks, but stronger for strange quarks (Fig. 4-3). It is suitable for 
the blue strange quarks when red and green strange quarks have paired in the 
(3 A ,35,1,+)(C7 3 ) channel. 
Many of the attractive channels have repulsive partners with the same symmetries, 
so a condensate in the attractive channel will automatically generate a small addi- 
tional one in the repulsive channel. For example, the (3/i,l/i,0,+)(C75) can generate 
(3^,1^,0,+) (C7075). This was discussed in Ref. [35], where the induced (C7075) con- 
densate (there called "k") in 2+1 flavor CFL was calculated and found to be small. 
In Ref. [113] it was observed that if all three quarks are massive then this condensate 
may be important. In the context of CFL the (3a,1aA+)(C75) can also generate 
(65,35,0,+)(C75) [11, 114], since they both break the full symmetry group down to 
the same subgroup. 
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Figure 4-1: Gap parameters in the attractive channels for an NJL interaction based on 
the two-flavor instanton. Since the instanton interaction requires two quark flavors, 
we take the quarks to be massless, which is a good approximation for the u and d. 
The cutoff is A = 800 MeV. 

4.3 Gap calculations for the attractive diquark chan- 
nels 

For the attractive channels we performed uncoupled gap equation calculations, and 
obtained the dependence of the quark pairing on fi. The amount of pairing is given 
by the gap parameter A(yu), which occurs in the self energy (See Appendix B.l) as 

A«f ab (p) = A(/i) C^^T ab , (4.5) 

with color matrix C, flavor matrix J 7 , and Dirac structure T ab . Note that A(yu) is a 
gap parameter, not the gap. It sets the scale of the gap in the quasiparticle excitation 
spectrum, but as we will see in Section 4.4 the gap itself often depends on the direction 
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Figure 4-2: Gap parameters in the attractive channels for an NJL interaction based 
on magnetic-gluon exchange. We show the one-flavor and two-flavor channels, for 
massless quarks. The cutoff is A = 800 MeV. 

of the momentum. 

The 4-fermion interactions that we use are nonrenormalizable, so our gap equa- 
tions (shown explicitly in Appendix B.2) involve a 3-momentum cutoff A, which 
represents the decoupling of our interactions at higher momentum, due to instanton 
form factors, effective gluon masses, etc. The usual procedure for NJL model calcu- 
lations is to calibrate the coupling strength for each cutoff A by known low-density 
physics such as the size of the chiral condensate. However, it is well known that this 
leads to an approximately cutoff-independent maximum gap (as a function of /i) in 
the ijjC'jsif) channel, so we used that criterion directly as our calibration condition, 
setting the maximum gap to 100 MeV. 

The results of our calculations, for cutoff A = 800 MeV, are plotted in Figs. 4- 
1, 4-2, and 4-3. For other cutoffs the overall shape of the curves is very similar. 
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Figure 4-3: Gap parameters in the attractive channels for an NJL interaction based 
on magnetic-gluon exchange, for quarks of mass 250 and 350 MeV, a reasonable range 
of values for the strange quark at medium density. We show only the single-flavor 
channels. The cutoff is A = 800 MeV. 

Because we use a sharp cutoff A, the gap falls to zero when /i reaches A (see, e.g., 
Eq. (B.16)). We show gap plots for the instanton interaction (Fig. 4-1) and magnetic 
gluon interaction. The full electric + magnetic gluon gives results that are similar to 
those for the magnetic gluon, but with no gap in the (6s,1a,1,+)(Co"o3) channel. 

For the magnetic gluon, we show a gap plot for massless quarks (Fig. 4-2) which in- 
cludes the two-flavor channels (3^,1^,0,+) (C75), (3a,1a,1,— XC7375) and (6s,1a,1,+) (CV 03 ) 
which could sustain u-d pairing, as well as the single-flavor channel (3a,3s,1,+) (C73) 
which could sustain u-u or d-d pairing. 

We also show a gap plot (Fig. 4-3) for a single quark flavor with mass of 250 or 
350 MeV. This is appropriate for the strange quark, whose effective mass is expected 
to lie between 150 MeV and about 400 MeV (see Ref. [26], Fig. 1). 
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The relative sizes of the gaps in the different channels reflect the pairing strengths 
given in Table 4.1. We see that the Lorentz scalar (3^,1^,0,+) (C75) (solid line) is 
dominant. The j = 1 channels have much smaller gap parameters. The (3^,3^,1,+) (C73) 
gap parameter (dash-dot line in Figs. 4-2 and 4-3) rises to a few MeV with the mag- 
netic gluon interaction. The (65,1^,1,+) (^003) gap parameter (dashed line) rises to 
about 1 MeV with an instanton interaction, but only 1 eV with the magnetic gluon 
interaction. It should be remembered, however, that the temperature of a compact 
star can be anything from tens of MeV at the time of the supernova to a few eV after 
millions of years, so gaps anywhere in this range are of potential phenomenological 
interest. 

The (65,35,0,+) (C7075) channel (dotted line), which is the only attractive channel 
for a single color and flavor of quark, is highly suppressed for massless quarks at high 
density but reaches about 10 keV for strange quarks (m = 350 MeV, Fig. 4-3). This 
is because its particle-particle component goes to zero as m — > (Eq. (B.33) and 
Table 4.1). 

Up to this point we have not mentioned the j — 1, rrij — ±1 channels (e.g. ipC , ~y±ip = 
■0C(7i ± ry 2 )?/>). We have only discussed the j — l,rrij — channels (e.g. tfjCj^ip). 
That is because rotational invariance of the interaction Hamiltonian that we are us- 
ing guarantees that changing rrij from to ±1 will not affect the binding energy and 
gap equation. This can be seen by considering the form of the binding energy. From 
Eq. (4.1) it is 

E B ~ (^y ac (^) bd H abcd . (4.6) 

Note that it is quadratic in the diquark condensate, with one of the factors be- 
ing complex conjugated. So if we have some 3-vector condensate, for example <fi = 
Si 0i(V'7i'0)) then its binding energy is 

^cxl^p + l^ + l^l 2 (4.7) 

It is clear that the rrij = condensate 0, = (0, 0, 1) has the same binding energy 
as the rrij = ±1 condensate (pi = (1/a/2)(1, ±2, 0). We have explicitly solved the 
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gap equations for the rrij = ±1 condensates, and find their solutions identical to the 
corresponding rrij = condensates. However, the quasiquark excitations in the two 
cases are quite different, and we proceed to study these in the next section. 

4.4 Quasiquark dispersion relations 

The physical behavior of quark matter will be dominated by its lowest energy exci- 
tations. As well as Goldstone bosons that arise from spontaneous breaking of global 
symmetries, there will be fermionic excitations of the quarks around the Fermi sur- 
face. In the presence of a diquark condensate, the spectrum of quark excitations is 
radically altered. Instead of arbitrarily low energy degrees of freedom, associated with 
the promotion of a quark from a state just below the Fermi surface to just above it, 
there is a minimum excitation energy (gap), above which the excitation spectrum is 
that of free quasiquarks, which are linear combinations of a particle and a hole. 

The dispersion relations of the quasiparticles can be calculated straightforwardly 
by including a condensate of the desired structure in the inverse propagator S^ 1 , 
shown in Eq. (B.3). Poles in the propagator correspond to zeros in S 1-1 , so the 
dispersion relations are obtained by solving det S~ 1 (po, p, fi, A, m) = for the energy 
Po as a function of the 3-momentum p of the quasiparticle, quark chemical potential 
/i, gap parameter A, and quark mass m. 

The gap is by definition the energy required to excite the lowest energy quasiquark 
mode. Isotropic condensates have a uniform gap, but one of the most interesting 
features of j > condensates is that they are not in general fully gapped: the gap goes 
to zero for particular values of momentum p, which correspond to particular places 
on the Fermi surface. This means that transport properties such as viscosities and 
emissivities, which are suppressed by factors of exp(— A/T) in phases with isotropic 
quark pairing, may not be so strongly suppressed by a j > condensate. In Figs. 4-4 
and 4-5 we show the variation of the gap over the Fermi surface by plotting the energy 
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of the lowest excitation function of angle, 



E glip (9) = mm\E l (p,9)\ (4.8) 

p,i 

where Ei(p, 9) is the energy of the i th quasiquark excitation with momentum (psin9 cos0,psin$sin0,pc 
For the plots we take fi = 500 MeV and A = 50 MeV, with quark mass m = (Fig. 4- 
4) or m = 250 MeV (Fig. 4-5). 




Figure 4-4: Energy gap in units of the gap parameter as a function of polar angle on 
the Fermi surface for rotational symmetry breaking phases with massless quarks, at 
fi = 500 MeV, gap parameter A = 50 MeV. 7+ = 71 + 272, 0+ = a i + ^02- 



1) C73 condensate: j — l,rrij — 0. 
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Polar angle 



Figure 4-5: Energy gap in units of the gap parameter as a function of polar angle 
on the Fermi surface for rotational symmetry breaking phases at fi = 500 MeV, with 
gap parameter A = 50 MeV. The quarks have mass m = 250 MeV. 7+ = 71 + 272, 

0"+ = O"01 + «C02- 

There is one quasiquark excitation with energy less than the gap parameter A. 

E{pf = (Vp 2 + mV/^iM 2 + A^ 
yUeffW 2 = /i 2 + A 2 cos 2 (£) (4.9) 
A cS (6) 2 = A 2 (sin 2 (£) W//i cff cos 2 (#)) 

From the expression for A c g(#) we see that for massless quarks the gap goes to zero 
for momenta parallel to the z-axis, i.e. at the poles on the Fermi surface (solid curve 
in Fig. 4-4). Massive quarks retain a small gap of order mA/fi at the poles (solid 
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curve in Fig. 4-5). 

2) C(ti ± 272) condensate: j = 1, rrij = ±1. 

There are two quasiquark excitations with energy less than 2 A, 

E{p) 2 = 2A 2 + m 2 + fi 2 + p 2 ± (4A 4 + 4fi 2 {p 2 + m 2 ) + 2A 2 p 2 (1 - cos(2fl)) 

V 1 (4.10) 

±4A 2 /i V /4m 2 + 2p 2 (l + cos (20))) 2 



For this condensate the effective gap again goes to zero at the poles, but in this case 
it remains zero even in the presence of a quark mass (dotted curve in Figs. 4-4,4-5). 

3) CV03 condensate: j = l,rrij = 0. 

There is one quasiquark excitation with energy less than A, its dispersion relation is 

[98] 

E{P? = WP 2 + ™V//4r ± + A 2 ff 

McffW 2 = /i 2 + A 2 sin 2 (£) (4.11) 
A cS (6) 2 = A 2 (cos 2 (0)+m 2 //i 2 ff sin 2 (0)) 

This is related to the dispersion relation for the C73 condensate by 9 — > tt/2 — 6: for 
massless quarks the quasiquarks are gapless around the equator of the Fermi sphere 
(dash-double-dot curve in Fig. 4-4) and in the presence of a quark mass they gain 
a small gap of order Am/ /j, (dash-double-dot curve in Fig. 4-5). The equator is a 
larger proportion of the Fermi surface than the poles, so in this case we might expect 
a greater effect on transport properties. 

4) C(eroi ± 2002) condensate: j = l,rrij = ±1. 

There are two quasiquark excitations with energy less than A. They have rather 
complicated dispersion relations. Going to the massless case, and assuming E, (p — 
/j) /i, which will be true for the low-energy quasiquark degrees of freedom that we 
are interested in, we find 



E(p) = (A 2 11 + A 2 p cos(0) ± 77/ a/2)/ (2/i 2 ) 

rj 2 = (8/i 4 (/i - p) 2 + 8A 2 /i 2 (/i 2 - cos(0) V) + 2A 4 (^ + cos(^)p) 2 ) 
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In this case there is a region near the poles, 9 < A//i, where the gap is zero (dashed 
curve in Fig. 4-4). This is because at those angles E(p) has zeros at two values of p 
close to fi. When 9 « A/^u those two zeros merge and disappear from the real p axis. 
The presence of a quark mass m > A wipes out this effect, but there is still no gap 
at the poles on the Fermi surface (dashed curve in Fig. 4-5). 

In comparing Figs. 4-4 and 4-5 it is interesting to note that introducing a mass 
for the quark opens up a gap whenever the gap lines intersect each other at a non- 
zero angle (after one includes the mirror-image negative-energy gap curves for the 
quasiholes). This occurs at zero energy at the poles for C73 and at the equator for 
Ca 03 . It occurs at non-zero energy for C(ji ± 172). The case of C(o" i ± io- Q2 ) is 
similar, but it is not obvious from the gap plot for reasons described above. 

We see that the J 7^ phases show a rich variety of quasiquark dispersion relations. 
For massless quarks they are all gapless in special regions of the Fermi surface, and 
for massive quarks the rrij = ±1 condensates remain gapless for momenta parallel 
to the spin. It follows that for these phases the quasiquark excitations will play an 
important role in transport properties, even when the temperature is less than the 
gap parameter. 

Moreover, different condensates (rrij = ±1 vs. rrij = 0) that because of rotational 
invariance of the Hamiltonian have exactly the same binding energy and gap equation, 
nevertheless have completely different energy gaps over the Fermi surface. They will 
therefore behave quite differently when exposed to nonisotropic external influences, 
such as magnetic fields or neutrino fluxes, and also in their coupling to external 
sources of torque, e.g. via electron-quasiquark scattering. All these influences are 
present in compact stars, and it will be interesting and complicated to sort out which 
is favored under naturally occurring conditions. And it should not be forgotten that 
these conditions vary with the age of the star. 
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4.5 Conclusions 



In this chapter we have completed a mean-field survey of 24 different channels for 
diquark condensation, using an NJL model with Fermi couplings for instantons and 
pointlike gluons. Our catalog of channels is by no means comprehensive, but it does 
exhaust the possibilities for structures that can be factorized into separate tensors 
for color, flavor, and spin (notably, this excludes interesting "locking" phases such as 
color-flavor-locking or color-spin- locking). 

As we promised in section 4.1, we can identify single-color and/or single-flavor 
phases that are useful for pairing "orphaned" quarks in the various pairing scenarios 
described at the start of this chapter. In fact our survey reveals five channels that 
are attractive and therefore susceptible to BCS condensation in the NJL model. In 
the same order that they are listed in table 1.1, these channels are: 



The first is well-known, the second is an obscure (and never-favored) spin-one cousin of 
the 2SC phase, and the last three are the channels that are good for pairing orphaned 
quarks. For these attractive channels we have calculated gaps and the results are 
shown in Figs. 4-1, 4-2, and 4-3. (Note that Fig. 4-2 also shows two more channels 
which are parity-partners of the 2SC and single-color channels. These channels are 
attractive for gluons, but they are parity odd and disfavored by instantons.) 

The gap of the 2SC phase has been calculated, by various methods, to be in the 
range of 10 to 100 MeV. In the present work we have assumed a maximum 2SC gap 
of 100 MeV and used this assumption to calibrate the NJL models. Then we find, 
as shown in Figs. 4-1, 4-2, and 4-3, that the "exotic" phases have gaps that are no 
larger than about 1 MeV and can be orders of magnitude smaller than this. Notice 
that while the 2SC gap curve is fairly flat on these plots, the gap curves for the 



single-color and single-flavor 



single-color 



single-flavor 



2SC 



spin-one 2SC 



(3 A ,1 A ,0,+)(C75) 
(3 A ,1 A ,1-)(C 7 37 5 ) 

(6 S ,lA,l,+)(C<703) 

(3 A ,3 5 ,1,+)(C7 3 ) 
(6 5 ,3 5 ,0,+)(C 7 o75) 
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exotic phases are quite steep. Because these are semilog plots this implies that the 
exotic gaps are drastically /i-dependent. They can change by more than two orders 
of magnitude when the chemical potential is increased from 400 MeV to 500 MeV. 

Moreover, the gap values are very sensitive to the choice of NJL model. For 
example, the gap for the single-color phase is seen to differ by six orders of magnitude 
when it is calculated with an instanton interaction (dashed line in Fig. 4-1) versus a 
magnetic gluon interaction (dashed line in Fig. 4-2). In both cases the NJL model 
is calibrated to give the same value (100 MeV) for the 2SC gap. Both the instanton 
and magnetic gluon interactions are reasonable models for some physics at moderate 
densities, but they yield very different results for the gap in this single-color phase. 

Buballa et al [98] observed similar difficulties in their NJL investigation of the 
single-color state. The model uncertainties cannot be resolved by appealing to a per- 
turbative calculation at asymptotically high densities, because the long-range mag- 
netic gluon which dominates at these densities [12] is always repulsive in the color- 
symmetric channel [63] and therefore does not predict a single-color pairing state. 

The calculation for the single-flavor phase is also very model-sensitive. However, in 
this case the state can be studied with a perturbative, model-independent calculation 
in the high density limit [63, 64]. In this regime we find the simple result that the 
spin-one (single-flavor) gap is A = cA , where A is the spin-zero (2SC) gap and c a 
constant (independent of n) that is somewhere between 0.002 and 0.01 (the constant 
is known exactly and depends on the particular spin ansatz, e.g. polar phase or color- 
spin-locking phase). The spin-one gap therefore varies less drastically with [i than in 
NJL models because it is just a constant fraction of the spin-zero gap. If we assume 
(without justification) that this relation holds at moderate densities, we predict a 
spin-one gap of 200 keV to 1 MeV, for a spin-one gap of 100 MeV. For fi of 400 to 
500 MeV, the NJL calculation predicts a much smaller gap that ranges from 0.1 to 
10 keV in this [i interval. 

In the high-density limit, the single-flavor condensate also includes an admixture 
of the C(To3 spin channel. This channel is repulsive in the NJL model but attractive 
for long-range magnetic gluons. In the high-density limit the color structure alone 
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dictates the attractiveness of the pairing channel: color-symmetric pairing is repul- 
sive and color-antisymmetric pairing is attractive, independent of color or spin [63]. 
But this simple conclusion is unlikely to persist at moderate densities when various 
in-medium effective interactions (like those of our NJL model) can supersede this 
perturbative result. 

The preceding discussion is intended to convey that our NJL calculations should 
be interpreted conservatively. They are useful because they indicate which pairing 
channels are most likely to be attractive at moderate densities that are of interest 
to us. But they provide only very rough estimates for the magnitudes of the gaps 
in the various channels. In our summary of pairing channels shown in table 1.1, we 
have estimated rough upper limits for the gaps by taking the peak values of these 
gaps from figures Figs. 4-1, 4-2, and 4-3. These numbers should be interpreted as 
reasonable upper bounds for gaps which could be orders of magnitude smaller. 

Although many of the channels we have studied have very small gaps and therefore 
very small critical temperatures, they will be phenomenologically relevant if they are 
the best pairing option available for some of the quarks. Since the temperature of a 
compact star falls to tens of eV when its age reaches millions of years, pairing can occur 
in such channels late in the star's life. The corresponding quasiquark excitations will 
become massive and this can suppress their participation in transport processes. With 
this in mind it would be useful to extend our analysis of the exotic condensates to finite 
temperature. Recent work indicates that the usual BCS relationship between the 
critical temperature and the gap parameter may be modified for spin-one condensates 
[98, 64]. 

The spin-one condensates are gapless at special 3-momenta (or nearly so with 
nonzero quark masses). Recall that the (j = l,rrij = ±1) condensates (^C^iir^)^ 
or i/jC(ctoi ±i<Jo2)?P) always have gapless quasiquarks at the poles of the Fermi sphere. 
The j = 1, rrij = condensates have gapless regions for massless quarks, at the poles 
(ipC^ip) or around the equator (ipCao^ip), but if the quarks are massive then the 
quasiquarks have a minimum gap of order mA/fi (Figs. 4-4, 4-5). The gapless or 
nearly gapless modes may continue to play a role in transport processes, even when 
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the temperature is much less than the gap parameter. It would be very useful to 
develop a transport theory for the spin-one condensates, to determine how the light 
modes contribute to neutrino emission/absorption via URCA processes or otherwise, 
and how they affect specific heat, viscosity, conductivities, etc. A natural first step 
would be to write down an effective theory, which would contain the lowest quasiquark 
modes and Goldstone bosons arising from the breaking of rotational symmetry (which 
could be called "spin waves" by analogy with helium-3 [97]). 

If spin-one color condensates occur in compact stars, they are likely to have in- 
teresting and complicated spatial and rotational textures like those of the various 
spin-one helium-3 phases. The j = 1 condensates can carry angular momentum sim- 
ply by aligning themselves in large domains, without involving any superfluid vortices, 
but it seems they will typically occur in conjunction with other phases that are su- 
perfluid. It would be interesting to investigate how the angular momentum is carried 
in this situation. 

Finally, we note that since the single-flavor color superconductor has a small gap 
and also exhibits an electromagnetic Meissner effect [96] (unlike the CFL phase [59]), 
it could have an upper critical magnetic field smaller than that of some pulsars. We 
can estimate the critical field by assuming that the critical magnetic energy density 
is comparable to the condensation energy of the color superconductor, i.e. B^/8n ~ 
/i 2 A 2 /7r 2 . By this estimate, pulsars with typical fields of 10 12 — 10 13 Gauss could 
destroy single-flavor pairing if the gap is smaller than 0.1 — 1 keV, and magnetars 
with fields of order 10 15 — 10 16 Gauss could destroy condensates with gaps of as large 
as 0.1 - 1 MeV. 
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Chapter 5 



Applications and Outlook 

5.1 Overview 

In this final chapter we take a closer look at two physical settings in which the crys- 
talline phase might occur with observable consequences. The first setting is astro- 
physical: the quark matter core of a rotating compact star, where a layer of crystalline 
quark matter can pin vortices and cause glitches. The second setting is terrestrial: an 
ultracold trapped gas of fermionic atoms, where a spin imbalance can induce the for- 
mation of a crystalline superfluid. In the last section, we close by describing unsolved 
problems and directions for future research. 

5.2 Pulsar glitches 

Many pulsars have been observed to glitch. Glitches are sudden jumps in rotation 
frequency Q which may be as large as AQ/Q ~ 1CT 6 , but may also be several orders of 
magnitude smaller. The frequency of observed glitches is statistically consistent with 
the hypothesis that all radio pulsars experience glitches [115]. Glitches are thought to 
originate from interactions between the rigid crust, somewhat more than a kilometer 
thick in a typical neutron star, and rotational vortices in the neutron superfluid. 
The inner kilometer of the crust consists of a rigid lattice of nuclei immersed in 
a neutron superfluid [116]. Because the pulsar is spinning, the neutron superfluid 
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(both within the inner crust and deeper inside the star) is threaded with a regular 
array of rotational vortices. As the pulsar's spin gradually slows due to emission 
of electromagnetic radiation, these vortices must gradually move outwards since the 
rotation frequency of a superffuid is proportional to the density of vortices. Deep 
within the star, the vortices are free to move outwards. In the crust, however, the 
vortices are pinned by their interaction with the nuclear lattice. What happens next 
varies from model to model. Perhaps the vortices exert sufficient force on the crust 
to tear it apart, resulting in a sudden breaking and rearrangement of the crust and 
a change in the moment of inertia [117]. Perhaps a large cluster of vortices within 
the inner crust builds up enough outward pressure to overcome the pinning force, 
suddenly becomes unpinned, and moves macroscopically outward [118, 119, 120, 121, 
122, 123, 124]. This sudden decrease in the angular momentum of the superffuid 
within the crust results in a sudden increase in angular momentum of the rigid crust 
itself, and hence a glitch. Perhaps, due to interactions between neutron vortices 
and proton flux tubes, the neutron vortices pile up just inside the inner crust before 
suddenly coming unpinned [125]. Although the models differ in important respects, 
all agree that the fundamental requirements are the presence of rotational vortices in 
a superffuid, and the presence of a rigid structure which impedes the motion of these 
vortices (by a pinning force of suitable magnitude) and which encompasses enough of 
the volume of the pulsar to contribute significantly to the total moment of inertia. 1 

Although it is premature to draw quantitative conclusions, it is interesting to 
speculate that some glitches may originate not at the crust, but deep within a pulsar 
which has a color superconducting quark matter core. If some region of the core is 
in the crystalline phase, because this phase is a superffuid it will be threaded by an 
array of rotational vortices, and these vortices may be pinned in the crystal by the 
periodic spatial modulation of the diquark condensate. The basic reasoning is that 

because the diquark condensate must vanish in the core of a rotational vortex, the 

1 The first model of glitches which was proposed [126] relies on the cracking and settling of the 
neutron star crust ( "starquakes" ) as the neutron star spins down. This model does not require the 
presence of rotational vortices. However, this model fails to explain the magnitude and frequency of 
glitches in the Vela pulsar [123, 124]. 
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vortex might prefer to reside at at a node of the LOFF crystal. It is interesting to note 
that enhanced pinning of magnetic flux tubes has been proposed as an experimental 
signature of the LOFF phase in an electron superconductor [127]. (However, the 
analogy may be misleading, because the underlying pinning mechanisms may differ 
for quark matter versus electrons in a solid, as we will discuss.) 

Before we delve into an assessment of the feasibility of quark matter glitching, it is 
worth noting that perhaps the most interesting consequence is for compact stars made 
entirely of quark matter (also commonly called "strange quark matter" because it will 
contain strange quarks in addition to up and down quarks). The work of Witten [128] 
and Farhi and Jaffe [129] suggests that strange quark matter may be energetically sta- 
ble relative to nuclear matter even at zero pressure. Then it might occur that observed 
compact stars are strange quark stars [130, 131] rather than neutron stars. This has 
recently been suggested for certain accreting compact stars [132], although the evi- 
dence is far from unambiguous [133]. 2 A conventional neutron star may feature a core 
made of strange quark matter, but strange quark stars are made (almost) entirely of 
quark matter with either no hadronic matter content at all or with a thin crust, of 
order one hundred meters thick, which contains no neutron superfluid [131, 137] (the 
nuclei in this thin crust are supported above the quark matter by electrostatic forces; 
these forces cannot support a neutron fluid). Because of the absence of superfluid 
neutrons, and because of the thinness of the crust, no successful models of glitches 
in the crust of a strange quark star have been proposed. Since pulsars are observed 
to glitch, the apparent lack of a glitch mechanism for strange quark stars has been 
the strongest argument that pulsars cannot be strange quark stars [138, 139, 140]. 
This conclusion must now be revisited: the quark matter in a strange star would be 
a color superconductor, and glitches may originate from pinning of vortices within a 
layer of the strange star which is in a crystalline color superconducting state 3 . 

2 Recently, the isolated compact star RXJ185635-3754 has also been proposed as a candidate 
strange star from a simple black-body fit to its observed thermal radiation [134], but this result is 
contradicted: the spectrum deviates from black-body, and when atmospheric effects are taken into 
account the inferred radius is consistent with conventional stellar models [135, 136]. 

3 Madsen's conclusion [141] that a strange quark star is prone to r-mode instability due to the 
absence of damping must also be revisited, since the relevant fluid oscillations may be damped within 
or at the boundary of a region of crystalline color superconductor. 
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Crystalline quark matter glitching is more interesting for strange stars than for 
other compact stars, because it is the only way for strange stars to glitch. Certainly, 
crystalline quark matter glitching could also occur in a more conventional neutron 
star with a quark matter core, but the core would be one of two potential locations 
where glitches may originate. Any attempt to observe a quark matter glitch would 
therefore require glitching models that are sophisticated enough to enable differentia- 
tion between the signatures of a quark matter core glitch and a nuclear matter crustal 
glitch. 

In this context the reader may be concerned that a glitch deep within the quark 
matter core of a neutron star may not be observable. However, due to electron 
scattering off vortices, the rotation of the superfluid interior of the star is coupled 
to the rotation of the electron plasma on a very short time scale (on the order of 
seconds), and the rotation of the electron plasma is subsequently coupled to the 
rotation of the outer crust on a similarly short time scale. Therefore the crust will 
rapidly respond to a change in the rotation of the interior superfluid [142]. This 
rapid coupling of the superfluid to the crust, due to the fact that the electron plasma 
penetrates throughout the star, is usually invoked to explain that the core nucleon 
superfluid speeds up quickly after a crustal glitch: the only long relaxation time is 
that of the vortices within the inner crust. Here, we invoke it to explain that the 
crust speeds up rapidly after a core glitch has accelerated the superfluid just outside 
the LOFF layer. After a glitch in the LOFF region, the only long relaxation times 
are those of the vortices in the LOFF region and in the inner crust. 

Our proposed quark matter glitch mechanism requires that the crystalline color 
superconducting phase occupies enough of the interior of a compact star to appre- 
ciably contribute to the star's moment of inertia. To address this question we need 
to complete a realistic three-flavor analysis of the crystalline state, and map out 
the n interval that the crystalline state will occupy in the real QCD phase diagram 
(figure 1-1). The extent of this interval will determine the extent of any volume 
of crystalline color superconductivity inside a compact star. So far, of course, we 
have only studied the crystalline state in a two-flavor model. In this model, we saw 
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that the plane wave crystalline state is favored only when 5n/A is in a small in- 
terval 5/i 2 /A ] ~ [0.707,0.754] (where A is the gap in the homogeneous 
BCS phase, i.e. 2SC), but when we extended our analysis to multiple plane wave 
crystals we found that the crystalline phase could occur over a much larger interval 
[5/i^/Ao, Sn^/A } (as sketched in figure 1-4). Extending the analysis to three flavors, 
we expect that the 5/j,/A interval of the two-flavor problem will correspond to an 
interval of m^/4//A in the three-flavor problem. If the plane wave crystal were fa- 
vored, it would only occupy a small sliver in the QCD phase diagram. The much 
more robust FCC crystal can occupy a much larger region, perhaps even the entire 
interval between the hadronization and unlocking transitions. As a function of in- 
creasing depth in a compact star, /i increases, m s decreases, and Ao also changes, and 
the crystalline phase will occur wherever m 2 s /Ajj, A is in the (large) interval where the 
crystalline phase is favored. Therefore while a quantitative analysis is yet to be done, 
we have reason to believe that the crystalline phase could encompass a significant 
volume of the quark matter core of a compact star. 

The next issue to resolve is whether the crystalline phase will actually pin vortices. 
A real demonstration of pinning, and a real calculation of the pinning force, will 
require explicitly constructing a vortex in the crystalline phase, and then observing 
how the energy of the vortex varies as it is moved across the crystal. As we shall 
discuss below, however, it is difficult to construct a vortex solution on top of the 
existing modulation of the crystalline condensate. The difficulty is related to the fact 
that the vortex is itself a deformation of the crystal, because both vortex and crystal 
are made of the same diquark condensate. It is easier to understand pinning in the 
nuclear crust, because there are two components involved: the first component (the 
lattice of nuclei) is a substrate that causes pinning in the second component (the 
neutron superfluid). It is the spatial variation of the substrate that pins superfluid 
vortices, not a crystalline variation of the superfluid order parameter itself. In other 
words, the superfluid does not "self-pin". In an electron superconductor, magnetic 
flux tubes are pinned to a substrate of random magnetic impurities. Enhancement of 
pinning in an electronic LOFF phase could occur by self-pinning of the LOFF crystal; 
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however, an alternative proposal is that the the LOFF phase does not self-pin but just 
makes it easier for the magnetic flux tubes to be pinned by the magnetic impurities. 4 
There is no "substrate" underlying quark matter, so if pinning does occur in the 
crystalline color superconductor, it must occur by self-pinning of superfluid vortices 
by the crystalline modulation of the same superfluid. 

Supposing that there is pinning, we can make an order-of-magnitude estimate of 
the pinning force which sidesteps the complications involved in explicitly constructing 
a vortex solution (at the end of this section we will return to the vortex problem and 
discuss some preliminary investigations). We perform a calculation similar to that 
done by Anderson and Itoh [118] for the pinning of neutron vortices in the inner 
crust of a neutron star. In that context, the pinning calculation has since been made 
quantitative [145, 122, 124]. We will attempt an estimate for both the plane wave 
LOFF phase and the more robust multiple-plane wave state (i.e. the FCC crystal). 
For the former, we can use numbers for the gap and free energy as in Figure 2-4. 
That is, we assume that at fi — 400 MeV we have a BCS gap A = 40 MeV and 
a BCS free energy of f2 = 2.6 x 10 7 MeV 4 at 5jj, = 0; then at <5/x 1 the gap for the 
plane-wave LOFF state is A pw ~ 8 MeV ~ 0.2Ao and the phase is favored over the 
normal state by a free energy f2 pw ~ 5 x 10 4 MeV 4 . The periodicity of the crystal 
is b pw = 7r/(2|q|) ~ 9 fm, and the thickness of a rotational vortex is given by the 
correlation length £ pw ~ 1/ A p w ~ 25 fm. For the FCC state, we cannot yet calculate 
numbers for the gap and free energy but we will use Af cc ~ Ao/2 = 20 MeV and 
fifcc ~ ^bcs/4 = 6.5 x 10 6 MeV 4 , reflecting our prediction that the FCC phase is 
robust. The node spacing is bi cc = \/37r/(2|q|) ~ 15 fm (i.e. half the lattice constant, 
equation (1.26)), and the thickness of a rotational vortex is £f cc ~ 1/Af cc ~ 10 fm. 

All these numbers are quite uncertain, but we will use them for the present. In the 

4 If the electronic LOFF condensate were to have a standing wave variation cos(2q • r) rather 
than an FCC structure, the condensate would vanish on nodal planes that are perpendicular to the 
magnetic flux tubes [127, 143] (the arrangement is similar to that of the vortex state in a layered 
high-T c cupratc superconductor with a magnetic field perpendicular to the Cu02 planes [144]). The 
appearance of nodal planes leads to a segmentation of vortices; the individual vortex segments are 
more flexible and are more readily pinned to random impurities in the substrate. In particular, the 
vortices would not be pinned to nodal sites of the LOFF crystal, in fact there would be no spatial 
variation of the LOFF crystal transverse to the flux tubes. 
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context of crustal neutron superfluid vortices, there are three distinct length scales: 
the vortex thickness £, the lattice spacing between nuclei b, and R, the radius of the 
individual nuclei. (The condensate vanishes within regions of size R separated by 
spacing b.) In the LOFF phase, the latter two length scales are comparable: since 
the condensate varies like a sum of plane waves it is as if R ~ b. The fact that these 
length scales are similar in the LOFF phase will complicate a quantitative calculation 
of the pinning energy; it makes our order of magnitude estimation easier, however. 
The pinning energy is the difference between the energy of a section of vortex of length 
b which is centered on a node of the LOFF crystal versus one which is centered on a 
maximum of the LOFF crystal. It is of order 

, 4 MeV plane wave 
E p ~ Qb 3 ~ { (5.1) 

3 GeV FCC crystal 
The resulting pinning force per unit length of vortex is of order 

E p I (5 MeV)/(100 fm 2 ) plane wave 
b [ (1 GeV) /( 100 fm 2 ) FCC crystal 

The fact that b and £ are comparable length scales will make a complete pinning force 
calculation more difficult and is likely to yield an f p which is significantly less than 
that we have obtained by dimensional analysis [122, 124]. Therefore these figures 
should be interpreted as upper-bound estimates. Note that our estimate of f p is 
quite uncertain both because it is only based on dimensional analysis and because 
the values of A, b and Q are known to only within an order of magnitude at best. For 
comparison, we present the corresponding numbers for the pinning of crustal neutron 
vortices: the pinning energy of neutron vortices in the inner crust is [122] 

E p « 1 - 3 MeV (5.3) 
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and the pinning force per unit length is [122, 123] 




1 - 3 MeV 



(5.4) 



(25- 50 fm)(4-20 fm) ' 



where the form of this expression is appropriate because £ < b. Of course it is 
premature to compare our crude results (5.1, 5.2) for quark matter pinning to the 
results (5.3, 5.4) of serious calculations for the pinning of crustal neutron vortices as 
in Refs. [145, 122, 124]. Nevertheless, we observe that the results are comparable for 
the case of a plane wave LOFF state, while pinning in the more robust FCC state 
could be orders of magnitude stronger than nuclear pinning (again, however, we are 
likely overestimating the pinning force so the numbers should be interpreted as upper 
bounds). 

The nuclear pinning force (5.4), when applied in a stellar model, does yield crustal 
glitches in accord with those observed in pulsars. It remains to be seen whether a 
much stronger pinning force, occurring in a crystalline region of the quark matter 
core, would also yield glitches consistent with observation. If the pinning force is 
large, then the restraint of vortices could be limited by the critical shear stress of the 
crystal: if the critical shear stress is exceeded, then the vortex can be released by 
a crystal dislocation. A dimensional analysis, like that which we have done for the 
pinning force, can only predict that the pinning force and the critical shear stress are 
of comparable magnitude, and no conclusion can be reached without a more careful 
investigation. 

We now return to the problem of constructing a vortex in the crystalline phase. 
Vortices are usually constructed beginning with a Ginzburg-Landau free energy func- 
tional written in terms of A(x) and VA(x). Instead, we have constructed a Ginzburg- 
Landau potential (equation 3.24) that is written in terms of the momentum modes 
A q of the order parameter . In principle, this contains the same information, but it 
is not well-suited to the analysis of a localized object like a vortex. In position space, 
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the Ginzburg-Landau free energy should look like 



fi[A(x)] = A(x)* [a + C(V 2 + 4g 2 ) 2 ] A(x) + C(A 4 ) 



(5.5) 



where C is a positive constant and a is the same parameter that appears in equa- 
tion (3.24), i.e. it is negative below the second-order point 5fi 2 - The quadratic term 
successfully reproduces the plane-wave instability of the LOFF phase: just below the 
second-order point, all the modes on the sphere |q| = go become unstable. Adding 
quartic and higher terms to equation (5.5) is difficult, however, because these terms 
are Fourier transforms of the complicated J and K functions defined in chapter 3. 
Terms with high powers of spatial derivatives will be required to reproduce the ef- 
fects of J and K and somehow yield a functional that is minimized by an FCC crystal 
(equation 1.25). So it is not at all clear that the usual methodology for construct- 
ing a vortex starting from a position-space Ginzburg-Landau potential is the right 
approach. 

As always, it should be easier to understand the physics of a vortex far from 
its core. The natural expectation is that far from its core, a vortex will be described 
simply by multiplying the A(x) of (1.25) by exp[z'0(x)], where 0(x) is a slowly varying 
function of x that winds once from to 2n as you follow a loop encircling the vortex 
at a large distance. In a uniform superfluid, this slowly varying phase describes a 
particle-number current flowing around the vortex. But in the crystal, the resulting 
particle-number current is 



and the current does not flow in a large loop because it vanishes at the nodal planes 
of the crystal (see figure 1-6). The dilemma is that the FCC crystal structure divides 
all of space into small cells bounded by intersecting nodal planes, and a supercurrent 
cannot flow across a nodal plane. In the presence of a vortex, the crystal structure 
must change to accommodate the rotational supercurrent flow. 

One obvious possibility is that the condensate just reduces to a one-dimensional 
standing wave cos(2q • x) with no variation in the transverse direction. This con- 




(5.6) 
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densate is just a stack of layers separated by nodal planes perpendicular to q, and 
vorticity is easily achieved (the supercurrent just circulates in each transverse layer). 
If this is what happens, then the vortices will not be pinned because they can move 
freely in the transverse directions. 

However, it may not be necessary to change the crystal structure so drastically. 
Suppose we instead consider changing the different A q 's by small phases, i.e. we 
consider condensates of the form 

A(x) = A q e i2q - X = Ae^e l2q ' x (5.7) 
q q 

(we assume that all the A q 's keep the same magnitude A). These are small distortions 
of the FCC ground state but perhaps these small distortions can change the nodal 
structure to accommodate vortices. For small phases, we can expand the free energy 
quadratically: 

fi({<M) = Vo + Yl ^q M qq'0q' + W 4 ) (5-8) 
q,q' 

Odd powers of 4> are not allowed because the free energy should be symmetric under 
A(x) — > A(x)*. The matrix M can be diagonalized to find eight normal modes for 
the phase angles. Four of these normal modes will have zero eigenvalue. One mode 
is a common phase for all eight q's (i.e. the superfluid mode). Recalling that the 
eight q's correspond to eight corners of a cube, the other three zero modes assign 
a common phase e tip to the four corners of one face of the cube, and an opposite 
phase e~ tip to the four corners of the opposite face. These are the three phonon zero 
modes; they correspond to translations in the x, y, or z directions by a displacement 
V5p/(2|q|). There is one nonzero "tetrahedral" mode which assigns a common phase 
to four tetrahedral corners of the cube (i.e. four nonadjacent corners) and the opposite 
phase to the other four corners. Finally, there are three "skew" modes which have 
the same nonzero eigenvalue of M; these modes assign a common phase to the four 
corners that are the endpoints of two opposite edges of the cube, and an opposite 
phase to the other four corners. 
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Figure 5-1: The effect of distorting the crystal by a phase shift, as in equation (5.7). 
On the left is the FCC crystal (phase shift 0); on the right is the distorted crystal 
(with a "skew" phase shift of magnitude 7r/8). The gap function A(x) vanishes on 
the gray nodal planes and on the heavy black nodal lines. The curved surfaces are 
contours where |A(x)| = 2.6A (note that for the undistorted crystal, |A(x)| is simple 
cubic even though A(x) is FCC). 

It is these last three "skew" modes that might be of interest. In figure 5-1 we see 
what happens to the crystal if a small phase is applied for one of these modes. The 
result is that only one set of parallel nodal planes remains, while the other nodal planes 
vanish and are replaced by nodal lines. With this deformation it should be possible 
to have circulating supercurrents in the layers between the remaining nodal planes, 
and the cost of the nodal deformation that accommodates this vorticity is small when 
the skew phase angle is small. Moreover, there is still a transverse variation of the 
condensate, and this could cause vortex pinning. In particular, the vortices might 
pin at the new nodal lines. In the deformed crystal, a current can be turned on by 
letting the different A q 's have different magnitudes. The next step is to slowly vary 
these magnitudes so that the current circulates in a large loop. This investigation is 
proceeding at the time of this writing. 
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5.3 Ultracold Fermi gases 



One of the most exciting prospects of the LOFF crystalline phase is the possibility 
that it could be observed as a "crystalline superfTuid" in a trap of ultracold fermionic 
atoms [106]. This could provide a laboratory test of our prediction that the LOFF 
phase has an FCC crystal structure. Even more significantly, the LOFF phase should 
be of interest to experimenters because it might in fact be easier to observe than the 
uniform phase: while the onset of uniform BCS superfluidity in the atomic system 
could be subtle to detect [146], the onset of crystallization could have a dramatic 
observational signature. Indeed, it might be possible to literally see the crystal struc- 
ture, as we will discuss. 

In section 1.7.2, we described how the two-flavor NJL model that we have used to 
investigate the crystalline color superconducting phase might be a surprisingly good 
model for the analysis of LOFF pairing in an atomic system. This is why we expect 
our prediction of an FCC crystal structure to apply in the atomic context. In this 
context, the two species that pair are the two hyperfine states of the fermionic isotope 
that is being trapped (experimenters have used 6 Li or 40 K [101]). The interaction 
between atoms, accomplished by a Feshbach resonance [104], is in the s-wave channel 
and is well-described as pointlike: the range R of the interatomic potential is quite 
small (R ~ 10 — 100 A, a typical van der Waals scale) compared to the atomic spacing 
achieved in the atom trap (on the order of 0.1 /mi), so the system is dilute. On the 
other hand, the scattering length a is quite large, of approximately the same order 
as the atomic spacing (the Duke group reports a ~ 0.5 fim [102]), which means that 
the system is strongly coupled [147]. The s-wave scattering length a is proportional 
to the four-fermi coupling constant g by the relation g = 4irh 2 a/m. 

The atoms are confined in a harmonic trap that is "cigar-shaped" with a very 
large aspect ratio (~30 for the CO2 laser trap used by the Duke group [102], and ~70 
for the magnetic trap used by the MIT group [148]). The transverse size of the cloud 
is a few to ten microns, while the longitudinal size is several hundred microns. The 
number density in the cloud varies somewhat with depth in the harmonic potential, 
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but a typical density is n = 10 14 cm -3 (for total number density of atoms in both 
hyperfine states). For now we assume that both states A and B are equally populated, 
i.e. ua = ub = n/2. Then both species have a common Fermi wave vector kp{= 
Pf/K) = (3n 2 nY^ 3 ~ 1.4 x 10 5 cm -1 , a Fermi energy Ep ~ 7 x 1CT 10 eV, and a Fermi 



Evaporative cooling has been used to obtain a highly-degenerate Fermi gas with a 
temperature perhaps as small as ~0.2 Tp [102] (there is some uncertainty in this figure 
because it is not easy to measure the temperature in the high-degeneracy regime). It 
is difficult to achieve temperatures very much smaller than Tp because the evaporative 
cooling efficiency is diminished by Pauli blocking [103]. The critical temperature for 
BCS superfluidity is [149] 



but this result is valid only in the weak-coupling regime {Jtzpa <C 1) where it predicts 
an exponential suppression of T c relative to Tp. The atomic systems are deliberately 
strongly coupled [kpa > 1) to avoid this exponential suppression. The strong coupling 
is achieved by tuning a magnetic field to the vicinity of a Feshbach resonance to obtain 
a very long scattering length (the Duke group reports kpa = 7.3 for their experimental 
conditions). In the strong coupling regime, the value of T c is not well known, but 
it is expected to be a substantial fraction of the Fermi temperature [105, 147]; a 
recent calculation [150] predicts that T c could be as large as 0.4Tp near the Feshbach 
resonance. 

It was recently proposed [151] that the superfluid phase transition can be identified 
from the free expansion of the cigar-shaped gas cloud after the trap is switched off. 
In particular, if the gas is in the normal (non-superfluid) state before the trap is 
switched off, then it will just expand ballistically and acquire a spherical shape. On 
the other hand, if the gas is in the superfluid state, then it will expand according to a 
hydrodynamic equation for the superfluid mode, leading to a pressure gradient and an 
enhanced expansion in the radial direction and an inversion of the cloud shape, i.e. the 



temperature Tp = Epfks — 10 /iK. The interatomic spacing is / ~ n 



1/3 



0.2 /im. 




(5.9) 
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cigar flattens into a pancake. A recent investigation by the Duke experimenters seems 
to produce this behavior, leading them to speculate that they may have observed the 
superfluid phase [102]. However, they also concede that collisional hydrodynamics of a 
non-superfluid gas could create the same inversion of the cloud shape. In the strongly- 
coupled system it is hard to rule out the influence of collisional hydrodynamics, so this 
may obscure any signature for the onset of superfluidity that involves the behavior of 
the expanding gas [146]. 

The crystalline superfluid might reveal itself in a more obvious way. To assess 
how we might observe the crystalline phase, we first make a rough estimate of the 
lattice spacing for the LOFF crystal in the atomic system. We write the BCS gap as 

A = CE F = Cf^ (5.10) 
2m 

where C is a number of order 1/10 to 1/2 that is poorly known. The LOFF phase is 
likely to occur when the separation of the Fermi energies of the two hyperfine states 
(A and B) is comparable to the value of the BCS gap, i.e. 

5E F = ^-^^v F 5p F ^A (5.11) 
2m 2m 

where v F is the Fermi velocity and Sp F = p F A — Pfb- To achieve a Fermi surface 
separation of this magnitude, the populations of the two hyperfine states A and B 
should be adjusted so that there is a fractional population difference 

N A -N B ^35p F ^3A _3 

N A + N B ~ 2 p F ~ 4 E F 4 y ' } 

where we have substituted equations (5.10) and (5.11). The momentum of a LOFF 
Cooper pair is 

g = 0-GSpF ~ 0.3Cp F , (5.13) 
again substituting equations (5.10) and (5.11). Finally, the spacing between nodal 
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planes of the FCC crystal is (see equation (1.26)) 



f*)*-„(3u (5.14) 



2g \C J k F \C 

where I = n _1//3 is the interatomic spacing. The interatomic spacing is about 0.2 /im, 
so we expect a nodal spacing 6=1 — 6 /im, for C = 1/2 — 1/10. 

All of these numbers are quite uncertain, first of all because they are derived from 
a weak-coupling analysis, whereas the atomic system is strongly coupled; and second 
of all because we have obtained our results at zero temperature, whereas the atomic 
system is likely to be observed near its critical temperature. Taking these numbers 
as rough estimates, however, it is interesting to consider how one might observe 
the crystalline phase. Our estimate of the nodal spacing suggests that the cigar- 
shaped cloud could have on the order of 1 to 10 nodes along the transverse width, 
and a few hundred nodes in the longitudinal direction. There will be a periodic 
modulation of the atom density in the crystalline superfluid: the density is depleted 
at the nodal regions of the crystal, in the same way that density is depleted in the 
core of a superfluid BCS vortex [152] (where the condensate also vanishes). The ratio 
of density modulation to total density is 

|Mx)|/n~|A(x)| 2 /^<C 2 (5.15) 

where C is the number from equation (5.10). Notice that since the modulation varies 
as |A| 2 it will actually be simple cubic rather than FCC. The density depletion is 
relatively modest (of order 10% or less). However, it was recently proposed [153] 
that density depletion in a vortex core might be enhanced at strong coupling, beyond 
the expectation of equation (5.15), to give a density reduction even as large as 70%. 
Although it has not been worked out, a similarly large depletion might occur at the 
nodal sites of the crystalline superfluid, leading to a dramatic density variation that 
could be easy to detect. 

The nodal spacing is perhaps too small to resolve by direct absorption imaging of 
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the in situ trapped gas (although it is worthwhile to investigate whether the crystal 
could be detected by Bragg scattering of infrared light). However, it might be possible 
to magnify the crystal structure by turning off the trapping potential and letting the 
gas expand until the structure can be seen. This is precisely the method that allowed 
the striking visualization of the vortex lattice in a rotating Bose-condensed gas [154]. 
In the BEC context, the expanding BEC cloud with its lattice of vortices is described 
by a time- dependent Gross-Pitaevskii equation [155]. For a cigar-shaped cloud rotated 
along its axis, the axial expansion is slower and the radial expansion is just a dilation of 
the vortex lattice. For a pancake-shaped cloud rotated along a normal axis, the axial 
expansion is faster and the radial expansion is not a dilation: the vortex core radius 
increases faster than the cloud radius, so as the vortex lattice expands the vortex holes 
become proportionally larger. These results for the expanding BEC vortex lattice are 
inspiration for an analogous calculation for the expanding LOFF crystal, to see how 
the density profile is deformed or dilated by the expansion. This is likely to require 
solving a time-dependent Ginzburg-Landau equation for the space-time evolution of 
the LOFF order parameter. 

If the gas does not expand like a fluid, but instead just expands ballistically when 
the trap is turned off, then the expanded gas is unlikely to bear a resemblance of the 
original crystal structure. In this scenario, the crystalline superfluid could be revealed 
by a measurement of the directional distribution of momenta in the expanded gas. 
The distribution should have an anisotropy reflecting the preferred arrangement of 
pairing rings on the Fermi surface, i.e. reflecting the fact that the FCC LOFF phase 
chooses to pair atoms in Cooper pairs with center-of-mass momenta pointing toward 
the eight corners of a cube. While a polar variation in the momentum distribution 
could always occur as a result of the initial axial confinement of the gas, any azimuthal 
variation would be an unmistakable signature of a crystalline phase. 
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5.4 Outlook 



Our investigation of the crystalline color superconducting phase and other non-CFL 
phases of quark matter has, not surprisingly, raised as many new questions as it 
has answered old ones. In this final section we enumerate some of the interesting 
directions for future work 

1. Surely the single most important problem to address is the fact that we have 
been unable to calculate a gap and free energy for the FCC crystal state. Our 
Ginzburg-Landau analysis is a powerful method that has taught us what fea- 
tures make a particular crystal structure energetically favorable or unfavorable, 
and therefore has led us to the FCC structure as the likely ground state. How- 
ever, the Ginzburg-Landau potential is unstable and therefore does not predict 
values for the gap and free energy. Given that, and given the prediction of a 
strong first-order phase transition, the Ginzburg-Landau approximation should 
now be discarded. What should be taken from our work is the prediction that 
the structure (1.25) is favored. Although A could be estimated by going to 
higher order in the Ginzburg-Landau approximation, a much better strategy is 
to do the calculation of A upon assuming the crystal structure (1.25) but with- 
out requiring A to be small. To do this will require constructing the anomalous 
propagator F(p,p') (3.14) by a resummation of the infinite geometric series 
shown diagrammatically in figure 3-3, rather than just keeping the first few 
terms in the geometric series as we have done in our Ginzburg-Landau analy- 
sis. The only difficulty is that F is not diagonal in momentum space; rather, 
it should be described as a vector of propagators _F q , one for each plane wave 
mode in the crystal, as in equation (3.14). Exactly resumming the geometric 
series of figure 3-3 will generate an infinite number of nonzero A q 's and cor- 
responding -Fq's, one for each q in the reciprocal lattice of the chosen crystal 
structure. A strategy that will make this calculation tractable is a truncation 
of this reciprocal lattice (analogous to the method of Ref. [72] to study chiral 
crystal structures), to lowest order just including the fundamental q vectors 
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(i.e. the eight q vectors that contribute to (1.25)), to next order including the 
next set of q vectors at larger radius in the reciprocal lattice space, etc. This 
calculation is proceeding at the time of this writing. 

2. If the crystalline color superconducting phase appears in nature, it will involve 
all three flavors of quark. We need to do a three-flavor analysis of the crystalline 
phase, including a realistic strange quark mass and imposing charge neutrality 
and weak equilibrium. Possible pairing strategies for the three flavor system 
include a CFL-LOFF phase, involving ud, us, and ds crystalline condensates, 
or a double-2SC-LOFF phase, involving only ud and ds pairings (i.e. pairings 
between species with adjacent Fermi surfaces). The three-flavor analysis is 
necessary to determine whether the crystalline phase is a superfluid, and it 
will also determine electromagnetic properties. Given the robustness of the 
crystalline color superconducting phase in the two-flavor system that we have 
studied, and its ability to occur over a large interval in Sli, we expect that 
crystalline color superconductivity will be a generic feature of nature's QCD 
phase diagram, occurring wherever quark matter that is not color-flavor-locked 
is to be found. 

3. We need to understand how to rotate a chunk of crystalline color supercon- 
ducting quark matter, i.e. we need to determine the structure of a rotational 
vortex in the crystalline phase. A determination of the vortex structure will 
resolve whether pinning can occur, and it is the first step toward a subsequent 
real calculation of the pinning force. As we have seen in section 5.2, even the 
long-distance behavior of the vortex is interesting to work out. 

4. Also related to pinning, we need to calculate the shear modulus for the crystal. 
Pinning of vortices in a crystalline color superconductor is likely to require 
both a nonzero pinning force and a nonzero shear modulus, because vortices 
can become unpinned by crystal dislocation if a critical shear stress is exceeded. 
The shear modulus could be determined from the phonon dispersion relation, 
which contains information about the elastic moduli of the crystal [48]. 
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5. The crystalline, single-flavor, and breached-pair color superconducting phases 
all have gapless quasiquark modes. If any of these phases were to appear in a 
compact star, these gapless modes might accommodate direct URCA reactions 
and play a critical role in the entire cooling history of the star. Therefore it is 
important to work out just how these special modes might contribute to cooling 
processes. 

6. In the atomic LOFF context, it will be crucial to study the free expansion of the 
atomic cloud after the trapping potential is shut off. For the scenario of fluid-like 
expansion, a calculation should be done in position space to determine how the 
crystal structure dilates or deforms as the cloud expands. For the case of ballistic 
expansion, a calculation should be done in momentum space to determine the 
anisotropy of the momentum distribution. It will also be useful to generalize 
our analysis to nonzero temperature, because the system temperature is likely 
to be near the critical temperature. 

7. Our analysis does not apply to QCD at asymptotically high density, where 
the QCD coupling becomes weak. In this regime, quark-quark scattering is 
dominated by gluon exchange and because the gluon propagator is almost un- 
screened, the scattering is dominated by forward scattering. This works in favor 
of crystalline color superconductivity [46], but it also has the consequence of 
reducing qo/Sfi and hence reducing ip . The authors of Ref. [46] find qo/Sfj, re- 
duced almost to 1, meaning i/jq reduced almost to zero. However, the authors 
of Ref. [47] find qo/Sfj, ~ 1.16 at asymptotically high density, meaning that 
ip ~ 61°. If the opening angle of the pairing rings on the Fermi surface does 
become very small at asymptotic densities, then the crystal structure there is 
certain to be qualitatively different from that which we have found. At present, 
the crystal structure at asymptotic densities is unresolved. This is worth pursu- 
ing, since it should ultimately be possible to begin with asymptotically free QCD 
(rather than a model thereof) and calculate the crystal structure at asymptotic 
density from first principles. (At these densities, the strange quark mass is ir- 
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relevant and a suitable Sfi would have to be introduced by hand.) Although 
such a first-principles analysis of the crystalline color superconducting state has 
a certain appeal, it should be noted that the asymptotic analysis of the CFL 
state seems to be quantitatively reliable only at densities that are more than 
fifteen orders of magnitude larger than those reached in compact stars [56]. At 
accessible densities, models like the one we have employed are at least as likely 
to be a good starting point. 



154 



Appendix A 

Evaluating J and K Integrals 



In this Appendix, we outline the explicit evaluation of the loop integrals in Eqs. (3.20) 
that occur in II, J and K. For all loop integrals, the momentum integration is 
restricted to modes near the Fermi surface by a cutoff uj <C p, meaning that the 
density of states can be taken as constant within the integration region: 

d 4 p = dp \p\ 2 d\p\ / dp^p 2 dp / ds dp (A.l) 

J — oo J Ji—uj J Air J — oo J —u> J Air 

where s = \p\ — p. Each integral is further simplified by removing the antiparticle 
poles from the bare propagators and that appear in the integrand. (We 
can disregard the antiparticles because their effect on the Fermi surface physics of 
interest is suppressed by of order A/ p.) To see how to remove the antiparticle poles, 
consider the propagator (tf + 2q( + that appears in the II integral. Recall that 

\i>u — ft — $H and we work in the limit where |q|,5/x <C ou <C |p| , /i. We are only 
interested in the behavior of the propagator in the vicinity of the particle poles where 
p° ~ ±(|p| — /2) <C /2. Therefore we can factor the denominator and drop subleading 
terms proportional to p°, Sfi, or |q| when they occur outside of the particle pole: 

1 = (p° + /i M )7°-(P + 2q)-7 

i> + 24+ fi u {p° + fi u -\p + 2q|)(p0 + p u + |p + 2q|) 

„ m° - p • 7 

~ (p° + p-6p- |p| -2q-p)(2/i) 
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We simplify all of the propagators in this way. In the numerator of each integrand we 
are then left with terms of the form 7 At 7 Q 7 /3 • ■ ■ 7 M . After evaluating these products of 
gamma matrices, the II integral can be written as 

n(q) = J ^ ^ J^da [(p° + s - - * - */* - 2q • p)] _1 • (A.3) 

This integral is straightforward to evaluate: Wick rotate p° — > 2^4, do a contour 
integration of the p^ integral, and then do the remaining simple integrals to obtain 
Eq. (3.21). 

By power counting, we see that while the II integral has a logarithmic dependence 
on the cutoff u, the J and K integrals have 1/uj 2 and l/u A cutoff dependences, 
respectively. We can therefore remove the cutoff dependence in the J and K integrals 
by taking the limit u/5fj,, ou/\q\ — > 00. Then the J and K integrals depend only on 
5fi and the q's and take the form 

^(qiq2qaq4) = 



1=1 



^(qiq2q3q4q5qe 

dp dp r+ao 



(A.4) 



where we have introduced new vectors 



ki = 0, k 2 = qi - q 2 , k 3 = qi - q 2 + q 3 - q 4 

£\ = qi, -^2 = qi-q2 + q3, £3 = qi - q2 + q3 - q4 + qs 

Notice that these vectors are the coordinates of vertices in the rhombus and hexagon 
shapes of Fig. 3-6. In particular, (k]k 2 ) and (€i£ 2 ) are the pairs of endpoints for the 
solid and dashed diagonals of the rhombus figure, while (kik 2 k 3 ) and [£\£i£z] are the 
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triplets of vertices of the solid and dashed triangles in the hexagon figure. 

We now introduce Feynman parameters to combine the denominator factors in 
Eqs. (A. 4). Two sets of Feynman parameters are used, one set for the factors involving 
kj's and one set for the factors involving £^s. For the J integral the result is 

J = / dxi dx 2 S(xi +X2 — 1) / dyi dy 2 5(yi + y 2 - 1) 
Jo Jo 

x / — — ds (s - Sfi + ip 4 + 2k • p)~ 2 + S/I- ip 4 + 2£ ■ p)^ 2 (A. 5) 

J 2lT 4:71 

where k = ^Vrcjkj, £ = YliVi^i- Next, we do the s integral by contour integration, 
followed by the p and p 4 integrals. For the P4 integral, noting that the s integration 
introduces a sign factor sgn(j>4) and that the integrand in (A. 5) depends only on ip^, 
we use 

/+oo POO 
dp 4 sgn(p 4 ) (•■■)= 2 Re / dp 4 ( • • ■ ) (A.6) 
-oo J e 

where e is an infinitesimal positive number. The final result is 



J = - Re j dxi dx 2 S(Y,x - 1) / dy 1 dy 2 S(^2y 
4 Jo Jo 



\k-£\ 2 -5fi 2 + 



(A.7) 



where 6fi + = 5fi + ie. We include the infinitesimal e is so that the integral is well- 
defined even when |k — £\ = 5fi is encountered in the integration region. This is a 
"principal value" specification for a multidimensional integral that is not Riemann- 
convergent. A similar analysis for the K integral gives the result 



K = - Re / dx 1 dx 2 dx 3 <5(£> - 1) / dy x dy 2 dy 3 5(J2u ~ 1) 
° Jo Jo 



|k-£| 2 + 3^/i 2 

{\k-£\ 2 -5n 2 + f 
(A.8) 

For the case of a single plane wave, where all the q,'s are equal (qi = q 2 = • • • = q), 
notice that \t\ = k 2 = k 3 = and £\ = £ 2 = £ 3 = q. Then the integrands in (A.7) 
and (A.8) are constants and we immediately obtain the results of (3.33) and (3.34), 
respectively. 

Finally, we must integrate the Feynmann parameters. For the J integral, two of 
the integrals can be done using the delta functions, a third can be done analytically, 
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and the final integration is done numerically, using an integration contour that avoids 
the singularity at |k — £\ = 5fi. For the K integral, we do the X3 and y 3 integrals 
using the delta functions, and then make a linear transformation of the remaining 
integration variables x±, x 2 , yi, y 2 , introducing new variables 

n = ^2 ai i x J + bijVj + c i> i = 1, . . . , 4 (A. 9) 

j j 

with ctij, bij and q chosen such that 

\k-£\ 2 = r 2 1 +rl + rl (A.10) 

While such a transformation puts the integrand in a convenient simple form, it com- 
plicates the description of the integration region considerably. Therefore we use a 
Fourier-Motzkin elimination procedure [156] to express the four-dimensional integra- 
tion region as a sum of subregions, for each of which we have an iterated integral with 
"affine" limits of integration: 

/ dxi I dx 2 / dy 1 / dy 2 ( • • • ) 
Jo Jo Jo Jo 

^[nf/tr^^l^l^-'l (A - n) 

A li=\ \ Ju lo +Ej<i% r i ) 

For each subregion A, we can immediately do the r 4 integration since the integrand is 
independent of r 4 . We are left with a three-dimensional integral over the volume of a 
polyhedron with six quadrilateral faces. We then apply the divergence theorem to turn 
the three-dimensional integral into a sum of surface integrals over the faces. For each 
surface integral, we convert to plane polar coordinates (p, 0) so that |k — £\ 2 = p 2 + d 2 , 
where d is the distance from the origin (n, r 2 , r^) = (0, 0, 0) to the plane of integration. 
Now the integration can be done because the integrand is independent of (p. Finally, 
the p integration is done numerically, using a deformed integration contour that avoids 
the singularity at |k — £\ = Sfi. 
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Appendix B 



Spin-One Calculations 



B.l Calculational details 

To allow the possibility of quark pairing, we use 8-component Nambu-Gorkov spinors, 

• =| m I (B.l) 

with 

* = (H P ),V(-P)) ■ (B.2) 
In Minkowski space the inverse quark propagator for massive fermions takes the form, 

S~\p)={ Mo (B.3) 

V A (^+m-/i7o) T / 

where 

A = 7o At 7o . (B.4) 

The gap matrix A is a matrix in color, flavor and Dirac space, multiplied by a gap 
parameter also denoted as A, 

Aff ab = A(/j)C af3 T ij T ab . (B.5) 
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The relation between the proper self energy and the full propagator is, 



/^-m + u7o \ / A \ 

S- 1 = S, 1 + £ = ( P + (B.6) 

V (tf+m-^o) T ' V A / 

where Sq 1 is the inverse propagator in the absence of interactions. The gap is deter- 
mined by solving a self-consistent Schwinger-Dyson equation for S. For a 4-fermion 
interaction modelling single gluon exchange, this takes the form 

S = J J^-VfSWV** (B.7) 

where V A is the interaction vertex in the Nambu-Gorkov basis. We study three 
interactions, the quark-gluon vertex 



Vu = ( "* ' m ) , (B.8 



l,X A /2 



-(7mA A /2) t 



the quark-magnetic gluon vertex 



7iX A /2 

y, A = A T , b.9 

-( 7t A A /2) r 



and the quark-instanton vertex, for which 



where 



■«G | -0^(v^{p)V^ + V&SgtoV^zSS (B.10) 



Sg] = -e ik ei l -(36%6] - fifty (B.ll) 



and 



^H ;i+75> r ), and ^C 1 ^ ° J. ,B, 2 , 
l + Tfef/ V 1- 7 5 T 
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In the case of the ipC^ijj condensate for the full gluon interaction we obtain the gap 
equation, which after rotation to Euclidean space becomes 



1 = 16G 



dp d 3 p 4 (A 2 + /j 2 + p 2 + p 2 
(2tt) 4 W 



(B.13) 



where 



W = A 4 + // + (p 2 + p 2 ) 2 + 2A 2 (/i 2 + p 2 + p 2 ) - 2 l 2 2 (-p 2 + p 2 ) . (B.14) 



The p integral can be explicitly evaluated, 



7T J 



P 



P 



'a 2 + (p + /i) 2 v/a 2 + (p - 



(B.15) 



The momentum integral can be performed analytically, giving 



A = 2^/A 2 ^ 



A 2 - 3/i 2 
1 — — ' 



7T 



4/i 2 G 



(B.16) 



for A < /i. 

B.2 Gap equation summary 

Here are the gap equations for the attractive channels. In the following, positive 
square roots are implied and we define p 2 = (p x ) 2 + ij>y) 2 - 



d\p\ = I d\p\, / dp r dp z = 

Jo J Jo 



dp r / dp z 



B.2.1 C75 and C gap equations 



1 = j d\ P \ 



IPI 



+ 



IpP 



A 2 + (H-/i) 2 v/A 2 + (|p| + / i) 2 



(B.17) 
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where N is a constant that differs for each interaction. 



Instanton N = A 

Magnetic + Electric Gluon N = 2 

Magnetic Gluon N — § 

The C channel produces an identical gap equation for both the full gluon and 
magnetic gluon interactions. The instanton interaction is not attractive in this chan- 
nel. 



B.2.2 Cao3 and CVT0375 gap equations 



with 



G 

1 = N— I dp r dp z 



El = A 2 + /1 2 + \p\ 2 ± 28 
where N is a constant that differs for each interaction. 



(B.l* 



Instanton 
Magnetic Gluon 



TV 
N 



The Co"o375 channel produces an identical gap equation for magnetic gluon in- 
teraction. The instanton and the full gluon interactions are not attractive in this 
channel. 



B.2.3 C(cr i ± icr 2) gap equation 

d Pr *. f *° (B.19) 

7T 3 7 i-00 VP 
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where 



W = // + {-{ P y + b| 2 ) 2 + 2A V - (p ) 2 - p\ - 2p%) - 2^{{ P y + | P n 



and N is a constant that differs for each interaction. 

Instanton 
Magnetic Gluon 



N = 2 
N = \ 



B.2.4 C73 gap equation 



G 



1 = iV — / dp r dp 



Pr{£ +P 2 Z ) , Pr{£-P 2 z 



£E_ 



with 



E 2 = A 2 p 2 z + ^ 2 (\p\ 2 + m 2 ) 
El = A 2 + /j 2 + m 2 + \p\ 2 ±2E 

where N is a constant that differs for each interaction. 



(B.20) 



Magnetic + Electric Gluon 
Magnetic Gluon 



N 
N 



B.2.5 C7375 gap equation 

This channel is not attractive for instantons and its gap equation with electric or 
magnetic gluon interaction is the same as for the massless C73 channel, i.e. 



1 = N 



G 



7T 



dp r dp z 



p r (£ +p 2 z ) P r (£-P 



EE A 



EE- 



(B.2i; 



with 



£ 2 = A 2 p 2 z + ^ 2 \p\ 2 
El = A 2 + /i 2 + \p\ 2 ±2S 



163 



where 

Magnetic + Electric Gluon N = | 

Magnetic Gluon N = \ 

B.2.6 C(7i ± ^72) gap equation 

1 = ^ / ^ dft f dp* ^ El _E l+ + (m 2 + , 2 - (/) 2 + rf)^) 



7T 3 7 rz 7_oo ^ J Bi + (4A 4 + 4A 2 (m 2 + /i 2 -(p0) 2 +p2) + jEl _ jBl+ ) 

(B.22) 



with 



£ 1± = m 2 -(/i±p°) 2 +p 2 
£ 2± = m 2 -(fi±p°) 2 + \p\ 2 

where N is a constant that differs for each interaction. 

Magnetic + Electric Gluon N = 1 

Magnetic Gluon N = ~ 



B.2.7 C7075 gap equation 



l = N^jd\ P \ 



\p\ 2 (£ + \ P \ 2 ) + \ P \ 2 (£-\p\ 2 ) 



EE, EE^ 



with 



E 2 = A 2 \p\ 2 + f i 2 (\p\ 2 + m 2 ) 
El = A 2 + /x 2 + m 2 + \p\ 2 ±2E 

where N is a constant that differs for each interaction. 

Magnetic + Electric Gluon N = \ 



Magnetic Gluon N = | 

For m = this reduces to 



'B.23) 



1 = - ^ (B.24) 

3vr 2 V ^T7? 



164 



B.3 Orbital/spin content of the condensates 

In the non-relativistic limit it is meaningful to ask about the separate contributions of 
the orbital and spin angular momenta to the total angular momentum of the diquark 
condensates. We can identify these by expanding the field operators out of which the 
condensates are built in terms of creation and annihilation operators, 



m 



k,s 



u s (k)at ia e-^ + v s (k)b s k ] a e 



:b.25) 



Inserting the explicit momentum- dependent spinors in any basis allows the creation/annihilation 
operator expansions of the condesates to be calculated. 
In the Dirac basis, 



/ 



u?(k) = A 



1 



Bk 3 
y B(k 1 + ik 2 ) J 



( 



<(k) = A 



-Bk 3 

1 



<(k) = A 



\ 



<(k) = A 




1 

B(h - ik 2 ) 
-Bk 3 
( Bk 3 \ 
B(ki + ik 2 ) 
1 

V 



(B.26) 



A 



E + m\ 2 



B 



1 



(B.27) 



2m J E + m 

Eq. (B.28) shows the result of performing such a calculation for the ipCip conden- 



sate, 



E 
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+2(a 1 pia a 2 _ pjB + b^l pjP )pz] S% (B.28) 

where S l £ 3 is the color-flavor matrix which is symmetric under the interchange i =± j 
(flavor) and a =± j3 (color). A sum over momentum p should be performed on the 
right hand side. 

Once the operator expansions have been obtained it is a relatively simple proce- 
dure to obtain the angular momentum content by rearranging the terms and inserting 
the relevant spherical harmonics. It is important to include contributions from mo- 
menta k and —k together, since they involve the same creation/annihilation operators. 
For example, for the condensate in (B.28), 

p = k : ^[al ia a\ jB {k 1 -ik 2 )-a\ ia a\ kjl3 {k l + ik 2 ) + 2^ (B.29) 
p=-k : - [-al^aljpih - ik 2 ) + a}_ kiot a\ j8 (k x + ik 2 ) - 2ai kia a 2 kjg k 3 ]{B.30) 
^ l a lia a -kjp{ki ~ ik 2 ) - alia al -kjf3(ki + ik 2 ) + 2a 2 kia a 1 _ kjg k 3 ] 

where we have relabelled k — > —k, i <-> j, a <-> j3 in the last line. The final result is a 
sum over k, a, (3, i, j of 

2 

^ [alia a -kjp(ki ~ ik 2 ) - afcia a -fcj/j( fc i + + (<Aia a -kjp + a li a a -kji3)h] (B.31) 

Upon inserting the relevant spherical harmonics and using standard arrow notation 
for the spins we obtain 



E 



73 ln>yfl + 7! lu)Fl1 " 737! [ln)+ liT)]Fl ° 



'B.32) 



which has precisely the correct Clebsch-Gordan structure to be interpreted as a state 
with orbital angular momentum 1 = 1, which gives an antisymmetric spatial wave 
function, and spin s = 1, which gives a symmetric spin wavefunction, combined to 
give j = 0. We write this as \l = 1a, s = Is). Applying this to all the condensates 
we studied, we can make a table of the particle-particle (as opposed to particle-hole) 



166 



content of each of them, 

These results are summarized in the "BCS-enhanced" column of Table 4.1. We can 
see explicitly that the ipCjo'-f^ip condensate has no particle-particle component in the 
massless limit, which is why it cannot occur at high density for the up and down 
quarks. This reflects basic physics: the condensate has spin zero, so the two spins 
must be oppositely aligned. But it is an "LR" condensate (see Table 4.1), so in the 
massless limit the two quarks, having opposite momentum and opposite helicity, have 
parallel spins. 
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4V2^ 



E 



4V27T 



71 p 

3E 
V 3 E 



5 ,s = A > 
1a, a = Is) 
s ,s = A ) 

lA, 8 = 1 S ) 

1a, s = Oa) 



ir 



-2zV2tt 



3 E(E + m) 

(£ + m) 2 -±p 2 " 
E(E + m) 



\l = 2 a ,8=l S ) 



\l = Q s ,s = l s ) 



ir 



-2V2tt 



3 E(E + m) 

(E + m) 2 - \p 2 ~ 
E(E + m) 



\l = 2 s ,s=l s ) 



\l = Q S ,8=l S ) 
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